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Zusammenfassung 



Ein Hauptzicl dicscr Arbeit ist die Bestimmung dcr mittlcrcn Matcricdichtc im Universum. 
Die Matcricdichtc ist ein wesentlichcr kosmologischcr Parameter, dcr die Zukunft des Univer- 
sums als Ganzem mitbestimmt. Es wurde dazu eine rontgenselektierte und rontgenfluiJbe- 
grenzte Stichprobe der 63 rontgenhellsten Galaxienhaufen am Himmel (ohne das galaktische 
Band, genannt HIFLUGCS) zusammengestellt, basierend auf der ROSAT Himmelsdurch- 
musterung. Die FluBgrenzc betragt 2 x 10""^^ ergs s^"*^ em^^ im Energieband 0.1 — 2.4 keV. 
Anhand mehrerer Tests wurde gezeigt, dai3 eine hohe Vollstandigkeit erreicht wurde. Diese 
Stichprobe kann, aufgrund der hoch angesetzten Flui3grenze, fiir eine Vielzahl von Anwen- 
dungen benutzt werden, die eine statistische Galaxienhaufenstichprobe benotigen, ohne Ko- 
rrekturen an das cffektive Durchmusterungsvolumen anbringen zu miisscn. 

Zur Bestimmung von Fliissen und physikahschen Haufenparamctcrn wurden hauptsachhch 
tief belichtete pointierte Beobachtungen benutzt. Es wurde gezeigt, dafi zwischen der Rontgen- 
leuchtkraft und der gravitativen Masse eine enge Korrelation besteht, wobei HIFLUGCS 
und cine crwcitcrc Stichprobe von 106 Galaxienhaufen benutzt wurde. Die Relation und die 
Streuung wurden quantifiziert mit Hilfe verschiedener Anpassungsmethoden. Ein Vergleich 
mit einfachen und erweiterten theoretischen und numerischen Vorhersagen zeigt insgesamt 
Ubereinstimmung. In grofien Rontgenhaufendurchmusterungen oder Simulationen dunkler 
Materie kann diese Relation direkt fiir Konvertierungen zwischen der Rontgenleuchtkraft 
und der gravitativen Masse angewendet werden. 

Daten des Galaxienhaufens Abell 1835, aufgenommen wahrend der 'Performance Verifi- 
cation' Phase des kiirzlich gestarteten Rontgensatellitenobservatoriums XMM-Newton wur- 
den ausgewertet, um die in dicscr Arbeit benutzte Annahmc, dafi das Intrahaufengas in den 
aufieren Gcbieten des Haufens isotherm ist, zu testen. Es wurde gefundcn, dafi das gemessene 
aufiere Temper aturprofil konsistent mit einem isothermen Profil ist. In den inneren Regionen 
wurde ein klarer Abfall der Gastemperatur um einen Faktor zwei gefunden. 

Physikalische Eigenschaften der Galaxienhaufenstichprobe wurden untersucht, indem Re- 
lationen zwischen verschiedenen Haufenparamctcrn analysiert wurden. Die Gesamteigen- 
schaften sind gut verstanden, aber im Detail ergaben sich Abweichungen von einfachen Er- 
wartungen. Es wurde gefunden, dafi der Anteil der Gasmasse an der Gesamtmasse nicht 
als Funktion der Temperatur des Intrahaufengases variiert. Fiir Galaxiengruppen (fceTx < 
2keV) wurde jedoch ein steiler Abfall dieses Anteils gefunden. Keine klare Tendenz fiir 
eine Variation des Oberflachenhelligkeitsprofils, d.h. P, als Funktion der Temperatur wurde 
beobachtet. Es wurde gefunden, dafi die Relation zwischen der Rontgenleuchtkraft und der 
Temperatur steiler als von einfachen selbstahnlichen Modellen erwartet verlauft, wie bereits 
in friiheren Arbeiten festgestellt. Allcrdings wurden keine klaren Abweichungen von der 
Form eines Potenzgesetzes bis zu einer gemessenen Gastemperatur /cb^x = 0.7 keV gefun- 
den. Die hier gefundene Relation zwischen der Gesamtmasse und der Temperatur ist steiler 
als von selbstahnlichen Modellen erwartet und die Normierung ist niedriger im Vergleich zu 
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hydrodynamischen Simulationcn, in Ubereinstimmung mit friihcrcn Resultatcii. Vorgeschla- 
gene Szenarien, daruntcr Hciz- und Kiihlprozesse, zur Bcschrcibung dieser Abweichungcn 
und Schwierigkeiten bei dem Bcobachtungsprozei3 wurden dargestellt. Es scheint, dal3 eine 
Uberlagerung verschiedener Effektc, moglicherweise inklusive einer Veranderung der mittleren 
Entstehungsrotverschiebung als Funktion der Galaxienhaufenmasse, benotigt ist, um die hier 
vorgelegten Beobachtungcn zu beschreiben. 

Unter Benutzung von HIFLUGCS wurde die gravitative Massenfunktion in dem Mass- 
enintervall 3.5 x 10^^ < M200 < 5.2 x 10^^ h^^ Mq bestimmt. Vergleich mit Press-Schechter 
Massenfunktionen fiihrte dazu, dafi die mittlere Materiedichte im Universum und die Ampli- 
tude der Dichtefluktuationen eng eingegrenzt werden konnten. Das groBe iiberdeckte Massen- 
intervall erlaubte einc individuelle Eingrenzung der Parameter. Im einzelnen wurde gefunden, 
dai3 rim = 0.12+[j;[j^ und as = OMtoH (90% Konfidenzintervall statistische Unsicherheit). 
Dieses Resultat ist konsistent mit zwei weiteren unterschiedlichen Abschatzungen fiir in 
dieser Arbeit. Der mittlere Anteil des Intrahaufengases an der Gesamtmasse von Galax- 
ienhaufen, bestimmt mit Hilfe einer erweiterten Sticliprobe von 106 Haufen, kombiniert mit 
Vorhersagen der Theorie der Elemententstehung deutet an, daB Jim ^ 0.34. Das Masse zu 
Licht Verhaltnis in den Haufen multipliziert mit der mittleren Leuchtkraftdichte impliziert 
~ 0.15. Eine Anzahl von Tests auf systematische Unsicherheiten wurde duchgefiihrt, 
darunter ein Vergleich der Prcss-Schecliter Massenfunktion mit den neuesten Resultaten von 
grol3en Vielteilchenrechnungen. Diese Tests ergaben Abweichungen kleiner als die statistis- 
chen Unsicherheiten. Zum Vergleich wurden die Werte der besten Anpassung von fiir 
gegebenes dg bestimmt, was zu der Relation cjg = 0.43 Jl""'^^ fiihrte. 

Die Massenfunktion wurde integriert, um den Anteil an der gesamten gravitativen Masse 
im Universum zu bestimmen, der in Galaxienhaufen enthalten ist. Normiert auf die kriti- 
sche Dichte ergab sich ^cluster = 0-012to'oo4 Galaxienhaufenmassen groBer als 6.4^06 ^ 
10^^ h'^Q Mq. Dies impliziert mit dem hier bestimmten Wert fiir Jim, dafi sich ca. 90% 
der Gesamtmasse des Universums aufierhalb von virialisicrten Haufenregionen befindet. Auf 
ahnliche Weise wurde gefunden, dafi der Anteil des Intrahaufengases an der Gesamtmasse des 
Universums mit J^b.Ciuster = 0.0015l|];[^[]2^ h^^-^ fiir Gasmassen groBer als e.Ql^lgX lO^^ ^^-5/2 
sehr klein ist. 
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1 Introduction 



Observational cosmology is currently a very active field of astronomical research. The aim is 
the determination of the past, present, and future status of the universe. Optical observa- 
tions of the magnitude-redshift relat ion of distant su pernovae (SNe) indicate an accelerating 
universe, i.e. 2 I^a > (^e.g-. iPerlmutter et al.lll999l ^. where is the normalized mean mat- 
ter density of the universe and Q\ the normalized cosmological constant. In the microwave 
regime, measurements of the fluctuations in the Cosmic Microwave Background (CMB) from 
recent satellite and ballo n borne experiments indicate a flat geometry, i.e. Om + ^^A = 1 (e.g., 
de Bernardis et al. 200(t ). From the X-ray side, clusters of galaxies, as the most massive 



collapsed objects known in the universe, are the ideal and most commonly used cosmological 
probes and indicate a low density universe, i.e. < 1, almost independent of ^}\. 

The SNe and CMB measurements clearly are very important cosmological tests, never- 
theless both suffer from inherent problems. The SNe are used as standard candles but the 
detailed physical processes that take place during a SN explosion are not well understood 
even at redshift z ~ 0. Furthermore additional dimming of the apparent brightness caused 
by gas and dust in the host galaxy of a distant SN is difficult to quantify. The presence of 
outliers in the SN distribution is worrying and a large number of objects is needed to clearly 
identify such eve nts fearlv results based on a smaller number of objects actually seemed to 
indicate i}\ ~ 0. iPerlmutter et al.lll997^ . CMB measurements yield information on the state 
of the universe at z ~ 1000, whereas galaxy clusters yield information from relatively low 
redshifts (0 < z < 1). It is thus important to take advantage of an independent method using 
galaxy clusters as performed in this work. 

Galaxy clusters can be used in a variety of ways to constrain cosmological parameters. 
For instance one may determine the typical mass to light ratio in clusters and multiply 
it by the measured total luminosity density to determine the mean mass density (e.g.. 



Bahcall et al. 19951 . and Sect. 17. 3|) . The underlying assumption is that the mass to light 
ratio in clusters is a good approximation to the mass to light ratio of the universe. An- 
other method uses the amount of mass contained in the intracluster gas as compared to 
the total gravitational cluster mass to set an upper limit on Q,^ b y comparison to pr edic- 
tions for the baryon density from the theory of nucleosynthesis (e.g.. White et al. 1993bl . and 
Sect. I7.H|1 . This approach assumes that the gas fraction in clusters resembles the baryon 
fraction in the universe. Furthermore within the framework of hierarchical structure forma- 
tion ~ small objects, e.g. galaxies, form first and assemble to larger structures, e.g. galaxy 
clust ers, afterwards under the influence of gravity - the merger rate depends also on 0,^ 
(e.g., Lacev and Colel 1993h . A comparison of observed cluster substructure frequencies with 
predictions of sp ecific models therefore allows in principal to put constraints on important 
parameters (e.g., Schuecker et al.l2001c^ . Moreover within this framework analytical prescrip- 
tions have been developped, which allow statistical predictions of the cosmic mass distribution 
for given cosmological models (e.g.. Press and Schechtei 19741 : Bond et al. 199ll : Boweilll99ll : 
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Lacev and Cold 19931 : Kitavama and Siitol 19961 : Schuecker et al. 2001a[ ). These predictions 



have been tested against a number of A^-bod y simulations and in general good agreement 
has been found over wide mass ranges fe.g.. lEfstathiou et al.lFl988l : iLacev and Colelll994l : 
Cover nato et al. 19991 : Ijenkins et a l. 2001)^. In this work the mean density of the universe 
is estimated utilizing three different methods. Most weight is given to the method, where 
the local cluster mass distribution is determined and compared to predictions, because the 
mass distribution is the most fundamental predicted quantity from cosmological models and 
simulations of structure formation. 

The observed galaxy cluster mass function, i.e. the cluster number density as a func- 
tion of mass, is particularly sensitive to t he fundamental cosmological parameters ari d 
amplitude of the density fluctuations (e.g., Henrv and ArnaudI 199l| *Bahcall an ' 
Previous local galaxy cluster mass functio ns ha ve been derived bvlBahcall and Ce: 



Biviano et al.Nl993h. Cirardi et al. ( 1998al ) , and Cirardi and CiuricinI (|2nnfll . for galaxv groups) . 



Bahcall and CenI ((1993) used the galaxy richness (measure of the cluster galaxy content. 



Sect. 



to relate to cluster masses for optical observati ons and an X-rav tempera ture 



mass relation t o convert the tempe r ature function given by iHenrv and ArnaudI (Il99lh to a 



mass function. Biviano et al. (jl99,'j ). Cirardi et al. ( 1998al ). and lCirardi and Ciuricin ip 




used velocity dispersions for optically selected samples to determine the mass function. Here 
a different strategy is used for the determination of the mass function. A statistical cluster 
sample is constructed taking advantage of the availability of an all-sky X-ray imaging survey. 
Furthermore the large number of archival cluster observations is exploited allowing detailed 
gravitational mass determinations through X-ray imaging and X-ray spectroscopy for the 
clusters included in the sample. 

The overall physical processes determining the main properties of clusters and their ap- 
pearance in X-rays are well understood. Intracluster gas (ICC) is trapped and heated to 
lO'^-lO^ K in the cluster gravitational potential. Thermal bremsstrahlung emission from this 
ICC makes clusters luminous X-ray sources. Only bright quasars exceed the typical cluster 
luminosities of ~ 10^^ ergss"^. 

The X-ray luminosity is well correlated with cluster mass (as will be shown in this work) 
as opposed to the measured galaxy richness, which has often been employed as selection 
criterium for optically selected cluster samples. Therefore X-ray selection effectively selects 
galaxy clusters by their mass. This property is vital for the construction of the mass function. 

The ROSAT All-Sky Survey (RASS), being the only all-sky su rvey carried out by an 
X-ray satellite with imaging capabilit ies to date (e.g.. Tn"impei]ll99,'j ) , has yielded a wealth of 
newly discovered X-ray sources fe.g.. lVoges et al. 1999h . A variety of galaxy cluster catalogs, 
together covering the whole sky, have been homogeneously built from the RASS (see refer- 
ences in Chap. El). These catalogs are utilized in this work for the construction of an X-ray 
selected and X-ray flux-limited sample of the brightest galaxy clusters in the sky (excluding a 
strip of ib20deg from the galactic plane, in order to ensure a high completeness). To perform 
a detailed characterization of the thereby selected galaxy clusters, including the determina- 
tion of the gravitational cluster mass, high quality (deep exposure) pointed observations of 
the PSPC detector onboard the ROSAT observatory are then analyzed and intracluster gas 
temperatures, mainly determined with the ASCA satellite because of its superior spectral 
resolution and larger sensitive energy range compared to ROSAT, are compiled from the 



more detailed discussion is presented in Sect. 12.2.3.21 
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literature. 

Within the framework of hierarchical structure formation the properties of galaxy clus- 
ters are also expected to follow certain scali ng relations and sim ulations have shown that 
scaled dark matter profiles look similar (e.g., Navarro et al. 1995h . Despite the well estab- 



lished overall understanding of galaxy clusters, in detail deviations from this picture have 
been found from observationally accessible quan tities, e.g., by the relation between X-ray 
luminosity and intracluster gas tempera ture fe.g., David et al.lll99^ ). and the gas properties 
in the center of groups of galaxies (e.g., IPonman et al.lll999l ). A variety of models has been 
suggested to explain these deviations (Sect.jT^)- Tests of detailed predictions of these models 
unfortunately are still compromised by observational difficulties. For instance the observed 
gas mass fraction has been found in the recent literature to either stay constant, decrease, 
or increase as a function of cluster temperature (Sect. [77I|l . The homogeneously selected and 
analyzed cluster sample presented here, comprising more than 100 galaxy groups and clusters, 
is therefore used to determine physical quantities like the X-ray luminosity, intracluster gas 
density distribution, temperature, and mass, as well as the gravitational mass over a wide 
temperature range from 0.7 to ISkeV. The relations between these quantities are analyzed 
and compared to predicted relations. 

The structure of this work is as follows. In Chap. |2l the different components of galaxy 
clusters are introduced with emphasis on the gas and gravitational mass determination. Fur- 
thermore some relevant cosmological background is given, including the calculation of model 
mass functions. Last not least the observing instruments are briefly described with most of 
the weight assigned to ROSAT and the RASS, according to their importance for this work. 
The sample construction is described in Chap. El The details of the data reduction and anal- 
ysis, and the determination of relevant quantities are given in Chap. 21 The gas temperature 
structure is very important for the X-ray mass determination of clusters. Therefore due to 
its importance for the present investigation the temperature profile for an example cluster is 
determined using brandnew data from the X-ray satellite mission XMM-Newton, which are 
in many respects superior to ROSAT and ASCA data (Chap.lS]). Before the mass function 
is determined in Chap. El - being the first galaxy cluster mass function constructed from an 
X-ray selected and X-ray flux-limited sample based on the RASS - the physical properties 
of the cluster sample are examined. Especially the correlation between X-ray luminosity and 
gravitational mass is of major importance here. In Chap. [3 tests of the sample complete- 
ness are performed and the cluster masses determined here are compared to independent 
determinations. The relations found between physical cluster properties are discussed. The 
mass function is compared to previous determinations and the cosmological implications of 
this mass function are presented, including a fit to model mass functions. Tight constraints 
on are derived. Previous work indicated that the mass fraction contained in galaxy 
clusters mav compri se already a fairly large fraction of the total mass in the universe (e.g., 
Fukugita et al. 1998l l. The well determined mass function given in this work is therefore used 



to determine the mass fraction in bound objects above a minimum mass to test these results. 
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2 Theoretical Background 



2.1 Galaxy Clusters 



Clusters of galaxies are believed to consist of four main components. As indicated by the 
name galaxy clusters have been discovered as conglomerates of galaxies. The space between 
these galaxies is not empty but contains huge amounts of intracluster gas (ICG). The largest 
portion of the total gravitating mass in clusters, however, exists in the form of dark matter. A 
possible forth component is a population of highly relativistic electrons, i.e. electrons having 
velocities close to the speed of light. Some charac teristics of the se components are briefly 



SarazinI 1986l l. Since this work mainly 



summarized in this Section (for a review see, e.g., 
deals with the intracluster gas and its implications for the dark matter content, these two 
components are awarded more attention. 



2.1.1 Cluster Galaxies 

How many galaxie s make a clu ster? An assembly of more than 4-5 galaxies is called a 
galaxy group (e.g., iHicksonI 11983 ) ■ ~ 100 galaxies make a cluster, and ~ 1 000 galaxies a rich 
cluster. These rough numbers exclude 'dwarf gal axies, which are difficult to count due to 
their faintness, except in the most nearby clusters. lAbelll (|1958') introduced the richness as a 
measure for clusters. The richness is determined by the number of galaxies above background 
fulfilling certain criteria. The two main criteria are that only galaxies be counted that a) are 
not more than two magnitudes fainter than the third brightest member galaxy, and b) have 
a projected distance from the center not larger than the Abell radius^ ta = 3 h'^^ Mpc. 

The galaxy population in clusters differs from the field population, i.e. galaxies not con- 
tained in clusters, especially in the following properties. 

• Morphology. The relative number of elliptical (E) and lens shaped (SO) galaxies in 
clust ers is larger and the relative num ber of spiral galaxies is smaller than in the field 



le.g. 



Dressleilll983 : lOemler ■TrJll992h . 



• Color. Spirals and irr egular galaxies iii clusters are redder on average than the same 
types in the field fe.g.. lOemler Jr. 19921 ). 



• Gas content. Especially spirals close to the cluster center con tain less amounts of neutral 
hydrogen than spirals in the field (e.g., Cavatte et al. 199nh . 



cD galaxies. These giant elliptical galaxies are found in the center of most groups and 
clusters. The most strik ing propertv of these c D galaxies is a very extended halo of low 
surface brightness (e.g., Matthews et al. 19641 ). 



^hso is defined in Sect. l2T2l Ipc = 3.085678 x 10^* cm. 
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2.1.2 Intracluster Gas 

The intracluster gas is the most massive visible component of galaxy clusters. Its mass ex- 
ceeds the (gravitating) mass contained in the cluster galaxies by a factor of ~ 2-5. The 
temperature, Tgas, of the ICG is in the range 1 < k^Tgas ^ 15keV (here IkeV corre- 
sponds to 1.16045 X 10^ K). The central gas number density, ngas(O), is in the range 10~^- 
10~^ particles cm~^. The collisionally ionized plasma is optically thin and emits thermal 
radiation in X-rays. For Tgas ^ 2keV the main component is bremsstrahlung (free-free 
transitions), for lower temperatures recombination (free-bound transitions) and line emission 
(bound-bound transitions) become more important. The emissivity depends on the density'^. 
A parameterized radial gas density distribution can be determined analytically from the ob- 
served surface bri ghtness distributio n. Numerical deprojections using onion shell models are 
also applied fe.g.. iFabian et aL 198ll ). The procedure of the analytic deprojection is outlined 



in Sect. 12.1.2.11 The observational determination of Tgas is described in Sect. 12.1.2.21 



2.1.2.1 Gas Density 

Assuming King 's (|l9fi2l 'l approximation to an isothermal sphere for the galaxy density dis- 
tribution, pgai, leads to an a nalytical representation of the radial gas density distribution 



(the 'standard /3 model', e.g., Cavaliere and Fusco-Femiano "197^; Sarazin and Bahcalll 1977 



Gorenstein et al. 19781 : Jones and Formaniiil984i : iSaraziniil986i) 



Pgas(O) 1 + -2 



(2.1) 



by using pgas oc p^^^, as implied by assuming the gas to be ideal, isothermal, and in hydrostatic 
equilibrium, and the galaxies to have an isotropic velocity dispersion, where /? denotes the 
ratio of the specific kinetic energies of the galaxies and the gas. The shape of the gas density 
distribution is therefore determined by the core radius, rc, and the shape parameter, p. The 
asumptions leading to the /? model may be violated in detail. The justification for its wide 
spread usage comes from the fact that the surface brightness profile derived from it (see 
below) represents the measured profile well in the relevant radial ranges. The gas mass. 



Mgs.s{< r) = 47r / Pgs^{r)r'^dr , 
Jo 

may for illustrative purposes be approximated for large radii and small (3 values by 

47r/,gas(0)r3 / .X -3/3+3 



(2.2) 



Mgas(< r) 



-3/3 + 3 



— > 1 A/3 < 1 



(2.3) 



The main constituents of the ICG are Hydrogen and Helium, where a good approximation for 
the number densities is nne = nn/lO. Due to the high temperature the gas can be considered 
completely ionized and the mean molecular weight including the electrons 



z>i 

2+ ^ wziZ + 1) 
z>i 



(2.4) 
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where Z is the atomic number and wz the relative weight (e.g., here W2 = 0.1 and wz = 
for Z > 2). Therefore one has fi ~ 0.61 and 

Pgas ~ l-17nemp. (2.5) 

Because of this proportionahty between electron number density, rig, and gas density it follows 
from ((TTJ) 

ne(r) = ne(0) 1 + ^ . (2.6) 



Before the connection between the observable surface brightness and the gas density is made 
a few more important quantities are introduced. The luminosity, L,^, i.e. the energy radiated 
per unit time at the frequency v is given by 



L^= / e^^dV, (2.7) 
Jv 

where the emissivity 

eu = nennA^{Tgi^,A) . (2.8) 

The emission coefficient, A^, mainly depends on the gas temperature and metalli city, A. 
Howe ver, it varies only weakly in the energy range where ROSAT is sensitive (e.g., B ohringei 
I995I ). for the relevant cluster gas temperature range (2~10keV). The emission measure is 



defined as 

E^= I nldV . (2.9) 



For the X-ray surface brightness, i.e. the number of photons detected in a defined energy 
range per unit time and per unit solid angle, one has 

/oo 
nidi, (2.10) 
-00 

where the integration is along the line of sight (/ = at the cluster center). With (|2.6j) it 
follows 

/■oo / 2 \ -3/3 

5xoc / 1 + - dl. (2.11) 

J —00 V J 

This integral can be reduced to a form solved in, e.g., Bronstein and Semendiaiewl (|l980l . 
Integral No. 39) and one finds 

Sx{R) = 5x(0) (^1 + > (2-12) 

where R denotes the projected distance from the cluster center. S'x(O) depends on ne(0), Tc, 
j3, Ajy(Tgas, ^), and redshift, z. Equation 1)2. 12() is used as a fitting formula to fit the observed 
surface brightness profile. With the obtained fit parameter values for S'x(O), r^, and (3 the 
gas density profile can be determined with 1)2. 1|1 . where Pgas(O) is obtained from (|2.5|) . The 
important step for the determination of the gas density distribution from p. 12)1 is the emission 
mechanism 1)2. 8p . which is well understood. The (3 model has been applied successfully already 
for many years, but also other models have been used, e.g., gas density distributions (e.g.. 
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Making et al. 19981 ) based on the Navarro-Prenk- White profile ( Navarro et al. 19961 . 199?! ). 
which is a fitting formula that represents well the cluster dark matter distribution found in 
iV-body simulations for varying cosmological models (but see Sect. 

Some clusters exhibit a central excess emission not well approximated by 1)2. 12() . To get a 
more accurate decription of the gas density profile in such double P model has been 

used by different authors (e.g., Ikebe et al. Il996l : Mohr et "all 19991 ) to fit the data, where 
the surface brightness takes the form Sx = Sx,i + Sx,2- The motivation is to have one 
component accounting for the central excess emission and the other component accounting 
for the overall cluster emission. It follows from the proportionality 1)2.10(1 that the gas density 
can then be determined from rig = [?^e,i'^ + 'T'e,2^]^^^- It has been shown, however, that the 
gas mass determination is not biased by the presence of central excess emission for instance 
byllleiprich (1998), who compared gas masses determined using single and double /3 models. 

It is worth noting that a new method to determine t he gas mass in clusters is becoming 
more and more important (e.g., ICarlstrom et al. l ll99(il V which uses the distortion of the 



CMB photon sp ectrum caused by inverse Compton scattering on the h ot ICG, the Sunyaev- 

Eaei; Isunva^ev and Ze1dovichlF97a) . 



Zeldovich effect ((Zeldovich and Sunvaev 



2.1.2.2 Gas Temperature 

When clusters of galaxies had been discovered as strong X-ray emitte rs more than 30 years 
ago (for references of the first detections and interpretations see, e.g.. ISarazinlll98(^ ) several 
possible emission mechanisms were di scussed. The d etection of li ne emission due to highly 
ionized iron in the X-ray spectra (e.g., Mitchell et al.l[l^7 6: Serlemit sos et al. 1977h . however, 
made clear that the major contribution is thermal emission. The main mechanism to heat 
the intracluster gas to the high temperatures observed is expected to be shocks. These shocks 
are caused by the gravitational assembling of the cluster from subunits. 

The electron temperature can be determined by fitting model spectra (folded with the 
instrument response) to the observed X-ray spectra (Chap. IS)). Assuming electrons and ions to 
be in ther mal equilibriuin this X -rav temperature corresponds to the gas temperature. Within 
2 h'^Q Mdc IFox and Loebl ( 1991^ ) have shown that this assumption should be satisfied. Since X- 
ray temperatures are seldom available for radii larger than 2 /i^q Mpc they should generally 
be good indicators of the gas te mperatures. Several spectral codes for hot, opticallv thin 
plasm as have been published (e.g. 

mi- 



3.1 SP 

Ravmond and Smithlll977l : iMewe et al.lll995l : ISmith et al 



The general dependence of the gas temperature on the distance from the clu ster center 
has been discussed controversely recently utilizing data from various satellites (e.g. , Rikazawal 
19971 : iMarkevitch et al.lll998l : llrwin et al.lll999l : Iwhitj bood : llrwin and BregmanI I200r h. In- 
cluding the latest findings from XMM-Newton (M. Arnaud, private communication; Chap.EJ 
the gas seems to be isothermal out to at least h alf the virial radius. In the very central part, 
where process es related to cooling flows (e.g.. 



1 11 994 and references therein) or cD 

galaxies (e.g., Mulchaev 2000l : Makishima et alJboOll . and references therein) may become 



important, a temperature drop is often found. 



2.1.3 Dark Matter 

The sum of the mass of all visible galaxies does by far not provide enough gravitational 
attraction to hold these galaxies in a cluster (e.g.. IZwickvlE933 ). Now, after the detection of 
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the large amounts of gas present in clusters, does the gas mass suffice to retain the galaxies 
and the gas? The answer is no. Assuming the laws of gravitation to be the same at the 
distance and at the scale of clusters still about 3/4 of the mass is 'missing'. 

Several candidates for this 'dark' matter have been and are being discussed. While, 
for instance, observations of the large scale clustering of galaxies rule out neutrinos (candi- 
dates for Hot Dar k Matter, HDM) as forming the only component of the dark matter (e.g., 
IWhite et al, the recent strong evidence that neutrinos with finite rest mass do exist 

(e.g.. IFukuda et al.. .l998i) leaves the possibility that at least part of the missing mass is pro- 
vided by neutrino s. One of the frequently cite d possible Cold Dark Matter (CDM) particles 
is the axion (e.g.. lOverduin and Wesson 199,'^: a, lso the heavier neutralino and gravitino are 
often discussed fe.g.. lOverduin and Wesson il997l ^ . 

Clusters of galaxies form a natural laboratory - obviously quite a bit larger than any 
experiment that could be built on Earth - filled abundantly with dark matter particles 
and may therefore be utili zed to actually place constrain ts on the nature of dark matter 



candidates. Recently, e.g., iSpergel and Steinhardtl (120001 ) suggested that elastic collisions 



of weakly self interacting particles may provide an explanation for the discrepancy between 
simulated CDM halos and observations of galaxies and clusters of galaxies. The discrepancy 
arise s when radial dark m atter profiles from simulations of coUisionless dark matter particles 



(e.g., Navarro et al.lll99^'l are compa red to d ark matter profi les indicated by rotation curves 
of dwarf galaxies (e.g.. iBurkert 19951: but se e lKravtsov et alJ 1998il and by radial gas density 
profiles of clusters (e.g., Makino et al. 19981 '!^. 



A more empirical dark matter density profile suggested bv iBurkertI (|l995h better repro- 
duces the data ori dwarf galaxies and also on the gas density distribution in clusters (e.g., 
Wu and Xujl2000l ^. 



This work mainly concentrates on the observational determination of the amount of grav- 
itational mass. Therefore in this paragraph a widely used method for this determination, 
which has also been used here, is decribed. The basic assumption is that the ICG is in 
hydrostatic equilibrium, i.e. 

dPge^jr) ^ -pgas(r)GMtot(< r) 
dr ' 
where Pgas represents the gas pressure, G the gravitational constant, and Mtot the cluster's 
gravitational mass. With the ideal gas equation, 

-fgas ~ Pgas^gas ; (^•1'^) 

/i?Tip 

this leads to 

Mto,« r) = (^_dp^ ^ I dT,^{r) \ ^^.15) 

firUpG \Pgas{r) dr Tgs^{r) dr J 

Inserting ()2.1|) and, based on the recent findings of XMM-Newton (Sect. I?!l.2.2() . assuming 
the cluster gas to be isothermal yields 

3A;BTgasr3/3 / 1 



Mtot(< r) = ° " ^— ^ . (2.16) 



^Note, however, that lYoshida et al] ll2000ll have shown that simulations, placed in a cosmological context, 
do not allow a simple model of dark matter particles with a finite cross section for elastic collisions to account 
for the discrepancy in dwarf galaxies and clusters simultaneously. 
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This equation is used throughout for the determination of the gravitational mass. The influ- 
ence of a possible non isothermality of the cluster gas on the results is discussed in Chap. [7| 
The assumption of hydrostatic equilibrium can be motivated by considering that the sound 
speed in the ICG Cg ~ 1000 km s~^ ~ IMpcGyr"^. The time a sound wave needs to cross 
the cluster therefore is small c ompared to th e cooling time and the time scales for sug- 
gested heating mechanisms (e.g., SarazinI 1986h . More quantitatively A^-body/hydrodynamic 
simulations have shown that as long as extreme merger situations are excluded, where the 
assumption of hydrostatic equilibrium probably breaks down, the mass estimates based on 
the hydrostatic equation give unbia. sed, accurate results with an uncertainty of 14-29 % (e.g., 
Schindleilll996al : Evrard et al. 199fih . The influence of bulk flows of the gas within the cluster 
and magnetic fields has been neglected. Cluster wide magnetic fields in a range of reasonable 
strengths about 1 /iC have been shown in ma gneto- hvdrodvnamic simula tions to provide only 
a minor non thermal pressure support (e.g., Dolae and Schindlei 2000l l . which is negligible 
for an X-ray mass determination based on (|2.15j) . 

Two other basic independent methods exist currently to estimate cluster gravitational 
masses. The oldest one utilizes the velocity dispersion of the cluster galaxies, and the youngest 
takes advantage of the alteration of images of background galaxies due to the gravitational 
field excerted by the cluster. The latter method can be subdivided into weak and strong 
gravitational lensing, the latter utilizing multiple and/or strongly distorted images of one or 
a few background galaxies and the former using a statistical approach on weak distortions 
of many background galaxies. Comparisons between the three mass determinations on low 
to medium redshift clusters seem to generally find good agreement between the velocity 
dispersion. X-ray, and weak lensing methods, whereas the strong lensing methods - probin, 
the ve rv center of clusters - yield factors of 2-4 higher masses (e.g., lAllen. ,1998. ,: .Wu et al 
1998lj ). The good agreement for the cluster masses at large radii is encouraging, since for 



the current work the overall cluster mass is important. In Sect. 17.21 masses for clusters 
determined here are compared to independent optical and X-ray estimates. The influence 
of weak systematic differences on the estimation of cosmological parameters is tested in 

Sect, mni 



2.1.4 Relativistic Electrons 



Extended diffuse synchrotron emission has been detected in radio imag es of gala xv clusters in- 
dicat ing the presence of highly relativistic electrons (e.g., Willsoniri97(]l . see, e.g., Govoni et al 



I2OO1I for a rece nt comparison of r adio and X-ray cluster properties) and large scale magnetic 
fields (see also lClarke et al.ll200ll ). Another indication is the detection of emission in excess 
of the bremsstrahlung and line emission expected for the intracluster gas temperature and 
metallicity on the soft (< 0.5 keV) and hard (> lOkeV) side, if inverse Compton scattering of 
CMB photons is the emission mechanism (e.g., [Sarazi nl ll999h . However, th e significance and 
especi allv the abundan ce of the soft exce ss is still debated vigorously (e.g.. lLieu et al.iri999l : 
bowver et al.l ll999: B ergh5fer et al.lll999l 'l. There seems to be less confusion about the hard 
excess being present in a few clusters (e.g.. iFusco-Femianol 1999h . Cluster mergers have been 
discussed as being responsible for the production of relativistic elec trons, since part of the 
released energy may go into particle acceleration (e.g., Observationally some 

indications for a correlation bet ween the presence of dif fuse radio emission and substructure 
in clusters have been found (e.g., Schuecker et aklbnOld ). The total energy contained in rela- 
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tivistic particles, however, is likely to be small compared to the total thermal energy conten t 
(~ 5 X 10^^ergss~^x ~ 10^'^ yr 2 x 10^^ ergs) of a typical cluster (e.g., |Sarazi 
Therefore the pressure supplied by these particles is negligible for the mass determination. 



2.2 Cosmology 

One of the main goals of the present work is to constrain important cosmological parameters. 
Sects. 12.2.11 and [2.2.21 giv e a brief overview o f the relevan t theoretical ba ckgroun d. See, for 
instan c e, .Mattig (1958). Wei^berd (Il972l. Sandagj Jioii) . ICarroll et al.l il992). Raschewskil 



r , ^ ' li^.^S). Weinberd (^9721. Sa^dl^ (BJ( 

|l995l ) , jPeacockl (|l999|) , and iGiulini and StraumannI ||2000|) for more extensive information 
and details. In Sect. l2.231 some observational and theoretical tools are given to turn knowledge 
about cluster gravitational masses into tests of cosmological models. 

2.2.1 Some Basics 

In an attempt to motivate the formulae used throughout this work the Einstein field equations, 
the Robertson- Walker Metric, the Friedmann-Lemaitre equation, and a few other important 
relations are introduced. 

The law of conservation of energy and momentum in special relativity (SR) takes the 
form 

^T^^ = a.T^''^ = , (2.17) 

where = {0,1,2,3} and, e.g., = ct, = x, x'^ = y, x^ = z are contravariant 
coordinates, where the vacuum speed of light c = 299792.458 km s"^ T^'•' is the energy- 
momentum tensor, describing the distribution and displacement of energy and momentum. 
In general relativity (GR) SR can be approximated locally, e.g., a free falling elevator of 
small dimensions (system of inertia) . If 1)2. 17() is to be written in GR for a general coordinate 
system then the partial derivatives have to be replaced by absolute (covariant) derivatives 

V^T^"' = 0. (2.18) 

Note that (|2.18j) is an invariant relation but in general cannot be interpreted as conservation 
of energy and momentum anymore since the energy and momentum induced by gravitation 
are not included in T^'^ . Now the hypothesis that the distribution and displacement of (non 
gravitational) energy and momentum is related to the geometry of space-time is introduced. 
Let the geometry be described by the Riemann curvature tensor, B7^y„, and the metric tensor, 
g^^'^ , then the simplest way to express this hypothesis and fulfil ()2.18|) is 

rp^^U ^ _ ^J^g^^U ^ Qt^U ^ (2.19) 

where the Ricci tensor R^j^y = R^^a ^'^^ the curvature scalar R = g^^R^u and G^" is the 
Einstein tensor. The simplest extension of (|2.19|) to fulfil 1)2. 18(1 is obtained by adding a 
constant 

T^^ oc G^''' + Ag^"" , (2.20) 

where A is the so called cosmological constant. The constant of proportionality is obtained 
by considering the limit of weak gravitational fields where Einstein's theory must go over to 
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Newton's. Einstein's field equations then read 

SttG 



c4 



T^"' = G^"" + Ag^"" . (2.21) 



A A > counteracts gravity and therefore theoretically allows a static universe, which was 
the main reason why Einstein introduced it into the field equations, before it was discovered 
that the universe actually is expanding. 

The next building block needed is the metric. Assuming that the universe is isotropic for 
every observer comoving with matter, the line element of the metric is given by 

ds^ = c^dt^ - Ritf + r'^dO'^ + sin^ Odct?] , (2.22) 

\ 1 — Kr^ J 

where R{t) is the scale factor at cosmic time t (not to be confused with the curvature scalar), 
r,6,(j) are unitless comoving coordinates, and the curvature index k = {0, ±1}. The metric 
defined by 1)2. 22() is often referred to as Roberts on- Walker metric. 

Inserting the metric tensor given by (|2.22j) in Einstein's field equations (|2.21j) yields the 
Priedmann-Lemaitre equation. This equation can, however, illustratively also be found within 
Newton's theory by considering a particle on a spherical shell with radius R encompassing a 
mass M{< R) 

]-R^ - G^ = const. , (2.23) 
2 R 

with 

M = 4/3Trp^^R^ , (2.24) 
where pm is the mean matter density within R, 

R^ —pmR^ = const. (2.25) 

follows. Including A und the constant one finds the Friedmann-Lemaitre equation 



i?2 SvrG A kc 



,2 



In the same picture another important relation is (illustratively) obtained by considering 

M 

R = -G^ , (2.27) 

with (fT^ 

47rG 

R = —pmR (2.28) 

follows. Including the pressure, and A one finds 



i2=-^(Pm + %)i?+^i?. (2.29) 
3 3 



Note that 1)2. 29(1 in general cannot be obtained directly by differentiating ()2.26|1 because pm 



may depend on t. Instead utilizing (|2.26j) and (|2.29j) . and setting = it is found that 
PraR^ = const., and 

^^m + ^^A + S^k = 1 , (2.30) 
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where the normahzed cosmic matter density Oin(t) = SuGp^it) / (2>H{ty'), the normahzed 
cosmological constant ^\{t) = A/{3H{t)'^), and the normahzed curvature index Ok(i) = 
-kc^ /{H{tfR{tf). The Hubble parameter H{t) = R{t)/R{t). For a flat universe {k = 0) 
one has Jim + f^A = 1- For the deceleration parameter go ^ —R{to)/[R{to)HQ), where 
Hq = H{to), one finds in a similar way 

Qo = - ^A,o , (2.31) 

where Om,o = ^m(io) and r2A,o = f^A(io)- For qo > {qo < 0) the expansion of the universe 
is decelerating (accelerating). The critical cosmic matter density is defined as 

Without subscripting from now on present day values, i.e. t = to = today, will be assumed 
for Qra, ^A, and fik- 

2.2.2 Practical Formulae 

For the physical distance, d, of two objects on the hypersphere of constant cosmic time, t, 
follows from 



R{t) arcsin r 
R{t)r 

R{t) arcsinh r 

For k = one therefore finds that the physical volume of a sphere with radius d is given by 



d = R{t)l ^J^ = { Rit)r 



k = +l, 

k = 0, (2.33) 
k = -l. 



V{r) = ^ TTR^{t) : = . (2.34) 

The redshift, z, compares a photon's wavelength measured in the observer rest frame, Aq, 
with the wavelength emitted in the source rest frame, Ai, by 

z^^V^. (2.35) 
Ai 

The redshift is caused by a change in the scale factor between the time of emission, ti, and 
absorption in the detector, to- Therefore one has 



The luminosity distance is defined as 



Dl^J^, (2.37) 
y 47r/Bol 

where /boI is the observed bolometric energy flux and LroI is the energy per unit time emitted 
by the source, the bolometric luminosity. Furthermore with R{to) = Rq 

Di^ = Ror{l + z). (2.38) 
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The angular diameter distance and the proper motion distance are related to Z^l by 

Da = Di^ {I + z)-^ and Dm = Di^ {I + z)'\ (2.39) 
respectively. The dependence of the Hubble parameter on redshift is given by 

H{zf = Hi E{zf , where E{z) = [17^(1 + zf + rJk(l + zf + n^f'^ . (2.40) 
Setting A = in 1)2. 26() the important Mattig equation can be derived for non vanishing p-^ 



c zqo + (go - l){y/2qo z + 1 - 1) 
Ho qi{l + z) 



Ror = — 2 /-I , -A ■ "A = 0Ai2m>0. (2.41) 



For JIa = and = 1 (|2.31j) yields qo = 1/2 and therefore 

2c / 1 



2c / 1 

i?or = — 1 - ) : 17a = A Jim = 1 ■ (2.42) 



Using (|2.38|1 one finds 



z = :^L»L-^ + y^I)L + ^ : f7A = 0Af)m = l. (2.43) 

The more general (valid also for ^ 0) relation between distance measure and redshift 
is given by 

r \n^\-^/^sm[\n^\y^Fiz;n^,nA)] k = +i, 
Dm = -^1 F{z;nrn,nA) 

m n-V2 -^1 roV2 



17//^sinh[17;^/^F(z;17m,17A)] 



fc = , (2.44) 
A; = -1, 



where 



F(z; ^a) = Hil + z'f{l + 17„,z') - z'(2 + z')nAr^^^dz' . (2.45) 
And comoving volumes for k ^ can be calculated by 

1, 

(2.46) 



^ 2//3l7k I GB^(l + aG'2 )^/2-J7i:'/'arcsinh(G^^f7^/') : k = -1 , 



where 



- ^ . (2.47) 



Remember that for k = one simply has V{Dm) = 4/3-7rD^ from 1)2.34(1 . 

Throughout Hq = 50 h^o kms^-*^ Mpc^-*^, h^Q = 1, = 1, f^A = 0, and = is assumed 
if not stated otherwise. Note that the determination of physical cluster parameters has a 
negligible dependence on and Oa for the small low redshift range under consideration 
here, as will be shown later. Therefore it is justified to determine the parameters for this 
specific model but to discuss the results also in the context of other models. 
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2.2.3 Mass Function 

A major motivation of this study is to determine the mean density of the universe, which 
is one of the key quantities that determines the fate of the universe. In this work this aim 
is achieved by determining an unprecedentedly accurate observational galaxy cluster mass 
function and comparing it to predicted mass functions. In this Section the basics for an 
observational determination of the mass function are outlined and it is shown how mass 
functions are predicted from cosmological models. 



2.2.3.1 Observational Determination 



The commonly used definition of the g alaxy cluster m ass function is analogous to the def- 
inition of the luminosity function (e.g., Schechtei 197d ): the mass function, (/>(M), denotes 



the number of clusters, A^, per unit comoving volume, dV ., per unit mass in the interval 
[M,M + dM], i.e. 

(/)(M) = N{M) /{dV dM) = dn{M) /dM . (2.48) 
Assuming con stant density the classical V^ax estimator (e.g., ISchmidtiFmfi^ : iFeltenI IToT^ : 



Binggeli et al.iil988, l can be used for the estimation of luminosity functions, i.e. 

1 ^ 1 

to = X7El7— • (2.49) 



T4nax is the maximum comoving volume within which a cluster with given luminosity for a 
given survey flux limit and sky coverage could have been detected. Combining ()2.34|) and 
(|2.38|) . and replacing 47r by the actual solid angle covered by the survey (this 'sky coverage' 
may in general be a function of flux), to, one finds that the surveyed volume enclosed by a 
cluster sitting at redshift z is given by 

where Dl is calculated by combining 1)2. 38(1 and 1)2. 42() . The maximum surveyed volume that 
a cluster with luminosity L could possibly enclose is given by 

Wi) = ^ (2.51) 

where -DL,max is determined by ()2.37|1 . replacing the measured flux with the flux limit, /um- 
-Zmax = -z(^L,max) IS given by (|2.43|) . Note, however, that in general only a finite energy band 
is available for flux measurements. Due to the different redshifts of detector [z = 0) and 
source the emitted spectrum for a given source rest frame (SF) energy range differs from the 
measured spectrum in the observer rest frame energy range (OF) and a correction factor has 
to be applied leading (for the ROSAT band) to 



rSF 

^[0.i-2.4kcv] ^ ^2.52) 



47ri<'(rgas,2:max)/i°^[0.i_2.4kcV] 



where K{Tgi^,z) = -^^|fi_2 4keV]/-^[o^-2 4keV] ^ cluster at redshift z. Unless noted other- 
wise in this work fluxes are quoted in the OF energy range and luminosities in the SF energy 
range. 
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2.2.3.2 Theoretical Determination 

In this Section the currently most important formahsm for the prediction of the abundance 
of massive objects is introduced. It is one of the basic ingredients to describe the growth 
of structure in the universe. This formahsm is used in Sect. 17.6.21 to obtain constraints on 
cosmological parameters. 

The following prescrip tion for the mass function (Eq. I2.53|) is based on the work of 



Press and Schechten (|1974f ). Vari ous assumptions of this ' PS' mass function have been relaxed 



by a number of works. See, e.g.. ISchuecker et al. for a compilation. The underlying 



idea of the extended PS formalism teond et al.lll99lh is to filter the initial mass density field 



on successively smaller scales. If a filtered region around the comoving location r exceeds a 
given density threshold it is assumed to end up as a collapsed object of the enclosed mass. 
Note that in this way the counting of objects contained in larger objects is avoided. This 
'cloud in cloud' problem was present in the original PS recipe and required the ad hoc intro- 
duction of a factor of 2. Under the assumption of a random initial density fluctuation field 
this picture is analogous to a random walk with an absorbing barrier, where the filter scale is 
interpreted as the time variable (in this picture time increases from large to small scales, i.e. 
from small to large fluctuations) and the density contrast, 6{r) = p[r)/p — 1, filtered on the 
corresponding filter scale, where p is the mean comoving density, is interpreted as the spatial 
coordinate. The mass fraction in collapsed objects above a minimum mass is associated with 
the fraction of volume elements which have a density contrast above the threshold. 

Three of the basic assumptions of this approach are the following. First, the initial density 
fluctuations are Gaussian and are uniquely described by their power spectrum, P{k), where 
k is the comoving wave number. Second, in order to get analytic results a top hat filter in 
Fourier space (sharp k space filter) is used. A third assumption, entering through the density 
contrast threshold, is that all objects form by a spherical collapse. 

As yet there see ms to be no comp elling observational evidence against the assumption 



mp 

of Gaussianity (e.g., IWu et al.ll2001al l. The second condition, wh ich requires una. t tractiv e 



oscillating filters in configuration space, has been relaxed recently bv lSchue^ker et alJ », 



who derived mass functions for more realistic filter functions. These new mass functions, 
however, are rather complicated to apply and the standard mass function has been used in 
this work to allow direct comparison to previous results. Note that these two assumptions 
imply that each point of each trajectory of a random walk has no memory of the past if 
followed from small to large scales^ (in this picture the time coordinate is reversed). This 
means a galaxy sized object at an early time before most clusters formed, s ay z = 3, has no 
information if it will end up in a cluster or in a void at z = (e.g.. White il996h . However, 
there is evidence that the formation of objects depends on their environment. For i nstance 
iV-bo dv simulations indicate that objects often align with large scale filaments (e.g., Whitel 



Il996h and observations show differences in the galaxy populations in the field and in clusters 
(Sect. I2.1.1j) . If this property of the extended PS formalism were realistic then possibly the 
infiuence of the environment would have to be only effective between z = 3 and z = 0. A 
model constructed completely free of this assumption might allow to better understand the 
discrepancies found if individual volumes (actually mass particles) are followed up in time in 



^Note that the use of, e.g., a Gaussian instead of the sharp k space filter does introd uce some correlati ons 
between different scales - with the dra wback of either having only numerical solutions iBond et alJll99ll) or 
a non negligible increase of complexity (|Schuecker et alJl2^^1a^ . 
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simulations and the final masses are compared to the predictions (e.g., Whitel 1996h . i.e. a 
comparison on the halo by halo basis, and possibly also the violation of the PS prediction that 
objects grow monotonically. As to the assumption of spherical collapse, simulations based 
on the hierarchical scenario rat her indicate that objects grow by accretioii of matter along 
filaments and/or merging (e.g., IColberg et Zl llQQfll : lOottlober etaP 1200 j). Many cluster 
syste ms haye also been observed which seem to be in a state of merging fe.g.- lMarkevitch et al. 
This and the violated predictions mentioned above have l ed to the su ggestion of 



replacing the spherical collapse model with ellipsoidal collapse (e.g.. lSheth et al.ll20 01). 

For the purpose of applying the analytical mass function as a fitting formula to observa- 
tional mass functions, these discrepancies are of minor importance as long as good agreement 
to simulated mass functions is shown, i.e. if there is good statistical agreement for simulated 
and predicted distributions . And really, good agreement has been found for many years 
re.g.. lEfstathiou et ablFmS^: fwhite et ajjlToQ^P: lliacev and ColellTflfli: Imo et a].lF99fiV Just 
recently large simulations covering very large mass ranges have conv i ncingly shown sligh t 
deviations at the high and low mass end (e.g., Governato et al. 199fll : Jenkins et al. 200 ll ^. 
However, the importance of these deviations for the current investigation is shown to be 
small in Sect. 17.6.21 In summary the justification for the usage of the standard formalism 
to be outlined below for the present work comes from the sufficient agreement with A^-body 
simulations. 

The PS formalism to predic t cluster mass funct ions for given cosmological models may 
be summarized as follows (e.g. 



Borgani et al. 1999h . To allow easier comparison with the 
theoretical literature on this subject in this Section /iioo = ^5o/2 is used. The mass function 
is then given by 



dn{M) 
dM 



2 pq 5c{z) 



vr M (T{Mf 



da{M) 



dM 



exp 



'2a{Mf 



Here M represents the object (cluster) virial mass and 



po = 2.7755 X lO" /i^oo ^0 Mpc" 



(2.53) 



(2.54) 



is the present mean matter density. The linear overdensity (density contrast threshold) com- 
puted at present 5c{z) = (^^(z) -D(O) D(z)^^, where the linear overdensity at the time of virial- 
ization , 5^ (-^)) is computed using the spherical collapse model summarized in lKitavama and Sutol 
(!l99fil ). 

3(12 7r)2/3 



and 



where 



3(12 7r)2/3 



20 



1.686 



20 



{1 0.0123 log [0/ (2)]} 



Om < 1 A rik = , 



an (1 + zf + (1 - - r?A) (1 + zf + 



(2.55) 
(2.56) 
(2.57) 



where z is the redshift of cluster formation. In this work the recent formation approximation 
is adopted and the observed cluster redshift is assumed to be the formation redshift. With the 
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assumptions ()2.56|) the second term in the denominator always vanishes. The hnear growth 
factor 



D{z) = 2.5 n^E{z) / {1 + z') Eiz')-'"^ dz' . 

J z 

The initial linear variance of the cosmic mass density fluctuations, 

1 



a{Ry 



k'P{k) \W{kR)ydk, 



where spherical symmetry has been assumed and the power spectrum is taken as 

P{k) oc k'^T{kf, 

may be expressed as 

joo r^Q^Y \W{k R{M))\^ dk 



a{My 



4 



f^k'^+nT{kY\W{kSh^^QMvc)\^dk ' 



(2.58) 



(2.59) 



(2.60) 



(2.61) 



where ag represents the amplitude of density fluctuations within a radius of 8 h^^Q Mpc. 
Recent measurements of the CMB an isotropics indicate that the primordial power spectral 
index, n, has a v a lue close to 1 (e.g., iBalhi et al. 20001 : .Taffe et al. 2001 : Prvke et al. 2001 



Wang et al 



2001 



de Bernardis et al.l 1200 ll ^ and is therefore set to 1 throughout. For the 



transfer function the fitting formula for CDM power spectra provided bv lBardeen et al. 
is used 



T{k) = T{q{k)) = ln(l + 2Mq)/{2Mq) 
X [1 + 3.89g + (IG.lqf + (5.46^)^ + (6.71^)^]-^/"^ 



(2.62) 



for 



q{k) = k/{ThiooMpc-'^), (2.63) 

where the shap e parameter is g iven by (modified to account for a small normalized baryon 
density fib > 0, ISugivam al ll99.5l i 

2 / 



m 'ilOO 



2.7K 



exp 



^100 
. 5 ^rn 



(2.64) 



The temperature of the CMB Tq = 2.726 K dMather et al.l Il994h and Qh /t?nn = 0.0193 
( Buries and Tvtlerlll998h . for the latter equation and (ITT)!!) /iioo = 0.71 (|Mould et alJl2000h 
will be used in the comparison to observed mass functions. The comoving filter radius 

R(M) = [3M/(47rpo)]^/^ (2.65) 

for the top hat filter function in configuration space, which is given in k space by 

Ty(x) = 3 (sinx — xcosx)/x^ : x = kR , (2.66) 

which is adopted in this analysis, because the cluster masses will be determined with a top 
hat filter, too (Sect. [2121) • I* is customary to use this filter despite the fact that the extended 
PS formalism, which predicts the correct normalization, has been derived using the sharp k 
space filter. Again the justification comes from the comparison to A^-body simulations. 

The two important parameters here, which will be constrained by a comparison to the 
observed mass function later, are the mean density of the universe, po, i.e. f^m when nor- 
malized by the critical density, and the amplitude of density fluctuations within a radius of 
8/i]"oqMpc, as- 
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Table 2.1. Instrument parameters 



Satellite 


Energy resolution 
FWHM [eV] 


Position resolution 
FWHM [arcsec] 


Field of view (diameter) 
[arcmin] 


ROSAT PSPC 


410 (J 


il.OkeV 


25@ l.OkeV 


120 


ROSAT HRI 






5 


40 


RASS 


410 « 


il.OkeV 


40@ l.OkeV 


unlimited 


ASCA GIS 


210 (c 


il.SkeV 


180=^ 


50 


XMM-Newton EPIC-pn 


110(c 


il.SkeV 


7@1.5keV 


30 



'^Half power diameter. 



2.3 Instruments 



This work is based on observations per formed by the X-ray satelli tes ROSAT (e.g., Triimper 



1993), ASCA fe.g.. lTanaka et al.lll994l 'l. and XMM-Newton fe.g.. I.Tansen et al.ltoll 'l. Since 



ROSAT is the most important instrument for this study, it is decribed briefly below. Some 
details on XMM-Newton are given in Chap. IHl 

A Delta rocket put ROSAT into an orbit 580 km above the Earth on June 1st 1990. After 
the calibration m easurements the first (and up to now last) All-Sky Survey (RASS, e.g., 
Voges et al.lll999h with an imaging X-ray telescope was performed with an average exposure 



of about 500 s. From 1991 till 1999 many much deeper exposures were taken of special fields 
(pointed observations) upon request by guest observers. 

The main component of ROSAT is the Wolter I telescope, which maps the X-ray photons 
onto one out of three focal plane detectors. It consists of four gold coated nested mirrors. 
The front part of the mirrors is slightly parabolically and the back part slightly hyperbolically 
shaped. Either one out of two proportional counters (Position Sensitive Proportional Counter, 
PSPC) or a channel plate detector (High Resolution Imager, HRI) can be placed in the 
focal plane. Additionally ROSAT carries a UV Camera (Wide Field Camera, WFC). The 
instrument parameters energy resolution, position resolution and field of view are given in 
Tab. 12. ll The effective area as a function of energy is compared to other missions in Fig. 12.11 
The high (particle induced) background rejection efficiency (> 99 %) and the large field of 
view of the PSPC allow to trace the source emission out to large apparent radii, an important 
feature for the study of extended objects like galaxy clusters. 

The average positional resolution of the RASS is lower compared to the on-axis resolution 
of the pointed observations performed with the PSPC (Tab.E^. This is due to a decreasing 
resolution with increasing off-axis angle and the fact that RASS photons have been collected 
at various different detector positions. Despite this lower resolution the nearby clusters 
which are relevant for this work still clearly appear as extended sources in the RASS which 
discriminates them from other X-ray sources like stars and AGN. 

Figure ITT] also shows that the major X-ray missions are sensitive exactly in the energy 
range where clusters of galaxies have their emission maximum. Because of this and their 
high luminosity the nearby galaxy clusters clearly stick out of the background. This can be 
appreciated in Fig. l2.2[ where the region of the Shapley supercluster of galaxies is shown. Most 



20 2. THEORETICAL BACKGROUND 



- ROSAT PSPC 

XMM EPIC PN+SxMOS 

. BeppoSAX SxMECS 

. Chandra ACIS-! 




Energy (keV) 



Figure 2.1 Comparison of the effective area of various instruments (Figure provided by V. 
Burwitz). Notice the log scale. 

of the dominant sources in this high cluster density region are galaxy clusters. On average it 
is estimated that about 10 % of all detected sources in the RASS are galaxy clusters. 
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Figure 2.2 Raw RASS image of the Shapley supercluster with names of the main individual 
galaxy clusters overlaid. The box size is 13.8 deg. The apparent large scale variations of the 
background are caused by differences in the exposure time (find a better version of this image 
at www.reiprich.net). 
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3 Sample 



The mass function measures the cluster number density as a function of mass. Therefore 
any cluster fulfilling the selection criteria and not included in the sample distorts the result 
systematically. It is then obvious that for the construction of the mass function it is vital to 
use a homogeneously selected and highly complete sample of objects, and additionally the 
selection must be closely related to cluster mass. In this work the RASS, where one single in- 
strument has surveyed the whole sky, has been chosen as the basis for the sample construction. 
Using the X-ray emission from the hot intracluster medium for cluster selection minimizes 
projection effects and the good correlation between X-ray luminosity and gravitational mass 
convincingly demonstrates that X-ray cluster surveys have the important property of being 
mass selective (Sect. I6T|) . 

Another aim of this work is the characterization of physical cluster properties. For such a 
statistical investigation it is important not to introduce any bias caused by selection effects. 
It is therefore desirable to select clusters by objective criteria, as performed here. 

Several cluster catalogs have already been constructed from the RASS with high com- 
pleteness down to low flux limits (see refs below). These have been utilized for the selection 
of candidates. Low thresholds have been set for the initial selection in order not to miss any 
cluster due to measurement uncertainties. These candidates have been homogeneously reana- 
lyzed, using higher quality ROSAT PSPC pointed observations whenever possible (Chap.Q. 
A flux limit well above the limit for candidate selection has then been applied to define the 
new flux-limited sample of the brightest clusters in the sky. 

In detail the candidates emerged from the following input catalogs. Table 13.11 lists the 
selection criteria and the number of clusters selected from each of the catalogs, that are 
contained in the final sample. 



1) The R EFLEX (ROSAT-ESO Flux-Limited X-ray) galaxy cluster survey (|Bohringer et al 
boflllj ) covers the southern hemisphere (declination 5 < +2.5 deg; galactic latitude 
|6ii| > 20.0 deg) with a flux limit /x,iim(0.1 - 2.4 keV) = 3.0 x lO'^^ g^gg g-i (,jn-2_ 



2) The NORAS (Northern ROSAT All-Sky) galaxy cluster survey (jRnhringer et a].ll200nl l 



contains clusters showing extended emission in the RASS in the northern hemisphere 
{S > O.Odeg; |6ii| > 20.0deg) with count rates Cx(0.1 - 2.4keV) > 0.06 ctss'^ 

3) NORAS II (J. Retzlaff et al., in prep.) is the continuation of the NORAS survey project. 
It includes point like sources and aims for a flux limit /x,iim(0.1 — 2.4 keV) = 2.0 x 
10~^^ ergss"^ cm~^. 



4) The BCS (ROSAT Brightest Cluster Sample) (jEbeling et al.lll998l ) covers the northern 



hemisphere {6 > O.Odeg; galactic latitude |6n| > 20.0 deg) with /x,iim(0.1 — 2.4keV) 
4.4 X 10^"*^^ ergss^""^ cm^^ and redshifts z < 0.3. 
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Table 3.1. Selection of candidates 



Sample 




Cx [ctss-i] 


fx [10-" erf 


;ss 1 cm ^] 




Ref. 




0.1 - 2.4 keV 


0.5-2.0 kcV 0.64 - 2.36 keV 


0.1 - 2.4 keV 


0.5 - 2.0 kcV 




REFLEX 




0.9 


1.7 




33 


1 


NORAS 




0.7 


1.7 




25 


2 


NORAS II 




0.7 


1.7 




4 


3 


BCS 


1.0 




1.7 




1 


4 


RASS 1 








1.0 





5 


XBACs 


1.0 




1.7 







6 


Abell/ACO 




0.7 









7 


early type 




0.7 









8 


previous saf* 













9 



^AU clusters from this catalog have been flagged as candidates. 



Note. — Only one of the criteria, count rate or flux, has to be met for a cluster to be selected as candidate. The 
catalogs are listed in search sequence, therefore Nqi gives the number of candidates additionally selected from the current 
catalog and contained in the final flux-limited sample. So in the case of NORAS a cluster is selected as candidate if it 
fulfils Cx(0.5 - 2.0 keV) > 0.7ctss-l or /x(0.1 - 2.4keV) > 1.7 x IQ-" ergss-^ cm-^ and has not already been selected 
from REFLEX. This candidate is counted under Nqi if it meets the selection criteria for HIFLUGCS. 



Re ferences, 
| 1998t (5) 
Il999t (9) 



)lBohringer et alj|2001bt (21 



De Grandi et a l. 1999: (6) Ebeling et al 
Lahav et al Higsg^-Edge et al...l99a . 



Bohringer et alj|2000l: (3) J. Retzlaff et al 
Il996t (7) H. Bohringer 1999 



m prep., 

private communication; ('Sl lBeuin^^^aT 



f4) lEbeling et al 



5) The RASS 1 Bright Sample of Clusters of Galaxies dPe Grandi et al.lll99fll 'l covers the 
south galactic cap region in the southern hemisphere {5 < +2.5 deg; bu < — 20.0deg) 
with an effective flux limit /x,iim(0.5 — 2.0 keV) between ~ 3 and 4 x 10~^^ ergss""^ cm~^. 



6) XBACs (X-ray Brightest Ab ell-type Clusters of gala- xies) ( Ebehng et al. 1996h i 



IS an 



all-sky sample of lAbelll (Il958h /ACQ dAbell et al.lll989l 'l clusters hmited to high galactic 
latitudes |6n| > 20.0 deg with nominal AGO redshifts z < 0.2 and X-ray fluxes /x(0.1 — 
2.4keV) > 5.0 x lO^^^ g^gss"^ cm'^. 

7) An all-sky list of Abell/ACO/ACO-supplementary clusters (H. Bohringer 1999, private 
communication) with count rates Cx(0.5 — 2.0 keV) > 0.6ctss~^. 

8) E a,rlv type galax i es wit h measured RASS count rates from a magnitude limited sample 
of lBeuing et!^ (Il999l 'l have been been checked in order not to miss any X-ray faint 



groups. 



All clusters from the sample of Lahav et al. ( 1989h and Edge et al. (|l990h , where clus- 
ters had been compiled from various X-ray missions, have been checked. 



The main criterion for candidate selection, a flux threshold 1.7 x 10~^^ ergs s~^ cm~^, 
i been chosen to allow for measurement uncertainties. E.g., for REFLEX clusters with 

X 10~^^ ergss"^ cm~^ < /x < 2.5 x 10~^^ ergss"^ cm~^ the mean statistical flux error is 
5 than 8%. With an additional mean systematic error of 6%, caused by underestimation 
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of fluxes due to the comparatively low RASS exposure times^, the flux threshold of 1.7 x 
10~^^ ergs s~^ cm~^ for candidate selection then ensures that no clusters are missed for a final 
flux limit /x,iim = 2.0 x 10~^^ ergss"^ cm~^. 

Almost none of the fluxes given in the input catalogs have been calculated using a mea- 
sured X-ray temperature, but mostly using gas temperatures estimated from an Lx-Tx rela- 
tion. In order to be independent of this additional uncertainty clusters have also be selected 
as candidates if they exceed a count rate threshold which corresponds to /x(0.1 — 2.4 keV) = 
2.0 X 10~^^ ergs s~^ cm~^ for a typical cluster temperature, Tgas = 4keV, and redshift, 
z = 0.05, and for an exceptionally high neutral hydrogen column density, e.g., in the NORAS 
case nh = 1.6 x lO^^cm"^. 

Most of the samples mentioned above excluded the area on the sky close to the galactic 
plane as well as the area of the Magellanic Clouds. In order to construct a highly complete 
sample from the candidate list the following selection criteria that successful clusters must 
fulfil have been applied: 

1) redetermined fiux /x(0.1 — 2.4 keV) > 2.0 x 10^-*^-^ ergss^^ cm^^, 

2) galactic latitude |6n| > 20.0 deg, 

3) projected position outside the excluded 324 deg^ area of the Magellanic Clouds (see 
Tab.ES), 

4) projected position outside the excluded 98 deg^ region of the Virgo galaxy cluster (see 
Tab. 1312)2. 

These selection criteria are fulfilled by 63 candidates. The advantages of the redetermined 
fiuxes over the fiuxes from the input catalogs are summarized in Chap. IH In Tab. 13. l1 one 
notes that 98% of all clusters in HIFLUGCS have been fiagged as candidates in REFLEX, 
NORAS, or in the candidate list for NORAS II; these surveys are not only all based on the 
RASS but all use the same algorithm for the count rate determination, further substantiating 
the homogeneous candidate selection for HIFLUGCS. 

The fraction of available ROSAT PSPC pointed observations for clusters included in HI- 
FLUGCS equals 86%. The actually used fraction is slightly reduced to 75% because some 
clusters appear extended beyond the PSPC field of view and therefore RASS data had to be 
used. The fraction of clusters with published ASCA temperatures equals 87%. If a lower 
fiux limit had been chosen the fraction of available PSPC pointed observations and published 
ASCA temperatures would have been decreased thereby increasing the uncertainties in the 
derived cluster parameters. Furthermore this value for the fiux limit ensures that no correc- 
tions, due to low exposure in the RASS or high galactic hydrogen column density, need to be 
applied for the effective area covered. This can be seen by the effective sky coverage in the 
REFLEX survey area: for a flux limit /x,iim(0.1 — 2.4keV) = 2.0 x 10~^^ ergss~^cm~2 and 

^This has been measured by comparing the count rates determined using pointed observations of clusters 
in this work to count rates for the same clusters determined in REFLEX and NORAS. If count rates are 
compared also for fainter clusters, not relevant for the present work, the mean systematic error increases to 
about 9% (BohriuEcr ct al. 2000). 

^The large scale X-ray background of the irregular and very extended X-ray emission of the Virgo cluster 
makes the undiscriminating detection/selection of clusters in this area diffcult. Candidates excluded due to 
this criterium are Virgo, M86, and M49. 
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Table 3.2. Regions of the sky not sampled in HIFLUGCS 



Region RA Range DEC Range Area (sr) 



LMC 1 


58 ^ 103° 


-63 - 


-77° 


0.0655 


LMC 2 


81 -> 89° 


-58 - 


-63° 


0.0060 


LMC 3 


103 108° 


-68 - 


-74° 


0.0030 


SMC 1 


358.5 20° 


-67.5 


^ -77° 


0.0189 


SMC 2 


356.5 ^ 358.5° 


-73 - 


-77° 


0.0006 


SMC 3 


20 30° 


-67.5 


^ -72° 


0.0047 


Virgo 


182.7 ^ 192.7° 


7.4 - 


^ 17.4° 


0.0297 


Milky Way'' 


^ 360° (111) 


-20 ^ 


20° 


4.2980 



^Galactic coordinates. 



Note. — Excised areas f or the Magellanic Clouds are the 
same as in iBohringer et al ] J2001bl) . because REFLEX forms 
the basic input catalog in the southern hemisphere. 



a minimum of 30 source counts the sky coverage amounts to 99%. The clear advantage is 
that the HIFLUGCS catalog can be used in a straightforward manner in statistical analyses, 
because the effective area is the same for all clusters and simply equals the covered solid angle 
on the sky. 

The distribution of clusters included in HIFLUGCS projected onto the sky is shown in 
Fig. 13.11 The sky coverage for the cluster sample equals 26 721.8 deg^ (8.13994 sr), about two 
thirds of the sky. The cluster names, coordinates and redshifts are listed in Tab. 14.11 Further 
properties of the cluster sample and completeness tests are discussed in Sect. |3 

For later analyses which do not necessarily require a complete sample, e.g., correlations 
between physical parameters, 43 clusters (not included in HIFLUGCS) from the candidate 
list have been combined with HIFLUGCS to form an 'extended sample' of 106 clusters. This 
sample is not a purely flux-limited sample with clearly defined selection criteria. Nevertheless 
this extended sample is not dominated by subjective selection and therefore one may take 
advantage of the increased statistics. The difference between the results of relations using 
HIFLUGCS and the extended sample is discussed in Sect. 16. ll 
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Figure 3.1 AitofF projection of the 63 HIFLUGCS galaxy clusters in galactic coordinates 
(filled circles). Additionally shown are 11 clusters above the flux limit but with |6n| < 20.0 deg 
(open triangles) and 32 clusters with fluxes below the flux limit (open squares). 



28 3. SAMPLE 



4 Reduction and Analysis 



This Chapter describes the derivation of the basic quantities in this work, e.g., count rates, 
fluxes, luminosities, and mass estimates for the galaxy clusters. These and other relevant 
cluster parameters are tabulated along with their uncertainties. 



4.1 Flux Determination 

Measuring the count rate of galaxy clusters is an important step in constructing a flux-limited 
cluster sample. The cou nt rate determi nation performed here is based on the growth curve 
analysis method (Bohrin ger et al.1 1200(1 ). with modifications adapted to the higher photon 
statistics available here. The main features of the method as well as the modifications are 
outlined below. 



The instrument used is the ROSAT PSPC ()Pfeffermann et al.Ml987^ . with a low inter- 
nal background ideally suited for this study which needs good signal to noise of the outer, 
low surface brightness regions of the clusters. Mainly pointed observations from the pub- 
lic archive at MPE have been used. If the cluster is extended beyond the PSPC field of 
view making a proper background determination difficult or if there is no pointed PSPC 
observation available, RASS data have been used. The ROSAT hard energy band (channels 
52 — 201 ~ 0.5 — 2.0keV) has been used for all count rate measurements because of the higher 
background in the soft band. The reduction of the soft X-ray background is about a factor of 



four i n the hard band, while still 60-100 % of the cluster emission is detected (jBohrineer et al 
l200(l ). Therefore the signal to noise ratio is multiplied by a factor 0.92-1.25 if the hard band 
is used. 

Two X-ray cluster centers are determined by finding the two-dimensional 'center of mass' 
of the photon distribution iteratively for an aperture radius of 3 and 7.5 arcmin around 
the starting position. The small aperture yields the center representing the position of the 
cluster's peak emission and therefore probably indicates the position where the cluster's 
potential well is deepest. This center is used for the regional selection, e.g. |6n| ^ 20.0. The 
more globally defined center with the larger aperture is used for the subsequent analysis tasks 
since for the mass determination it is most important to have a good estimate of the slope 
of the surface brightness profile in the outer parts of the cluster. 

The background surface brightness is determined in a ring outside the cluster emission. To 
minimize the infiuence of discrete sources the ring is subdivided into twelve parts of equal area 
and a sigma clipping is performed. To determine the count rate the area around the global 
center is divided into concentric rings. For pointed observations 200 rings with a width of 15 
arcsec each are used. Due to the lower photon statistics a width of 30 arcsec is used for RASS 
data and the number of rings depends on the field size extracted (100 — 300 rings for field sizes 
of 2 X 2deg^ — 8 X Sdeg^). Each photon is divided by the vignetting and deadtime corrected 
exposure time of the skypixel where it has been detected and these ratios are summed up 
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in each ring yielding the ring count rate. From this value the background count rate for the 
respective ring area is subtracted yielding a source ring count rate. These individual source 
ring count rates are integrated with increasing radius yielding the (cumulative) source count 
rate for a given radius (Fig. 14. I'll . Obvious contaminating point sources have been excluded 
manually. The cut-out regions have then been assigned the average surface brightness of the 
ring. If a cluster has been found to be clearly made up of two major components, for instance 
A3395n/s, these components have been treated separately. This procedure ensures that 
double clusters are not treated as a single entity for which spherical symmetry is assumed. For 
the same reason strong subtructure has been excluded in the same manner as contaminating 
point sources. In this work the aim is to characterize all cluster properties consistently and 
homogeneously. Therefore if strong substructure is identified then it is excluded for the 
flux/luminosity and mass determination. 



20 

rad 1 us_ 



30 

( arcm i n ) 



Figure 4.1 Cumulative source count rate as a function of radius (solid line) for the cluster 
A2029 (pointed observation). The vertical dashed line indicates the outer significance radius, 
rx. The dashed lines just above and below the source count rate indicate the l-cr Poissonian 
errors. 



An outer significance radius of the cluster, rx, is determined at the position from where 
on the Poissonian l-cr error rises faster than the source count rate. Usually the source count 
rate settles into a nearly horizontal line for radii larger than rx- In some cases, however, the 
source count rate seems to increase or decrease roughly quadratically for radii larger than rx 
indicating a possibly under- or overestimated background fFig. 14.2]) . Therefore a parabola of 
the form y = mx^ + b has been fitted to the source count rate for radii larger than rx and 
the measured background has been corrected. An example for a corrected source count rate 
profile is shown in Fig. 14.31 

Figure 1^31 shows for the extended sample (106 clusters) that the difference between mea- 
sured and corrected source count rate is generally very small. Nevertheless an inspection 
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Figure 4.2 Cumulative source count rate as a function of radius for the cluster EXO0422, 
shown as an extreme example (RASS data). The parabolic dashed line indicates the best fit 
parabola for count rates larger than rx- 




Figure 4.3 Corrected cumulative source count rate as a function of radius for the cluster 
EXO0422. The count rate correction is less than 5 %. 
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of each count rate profile has been performed, to decide whether the measured or corrected 
count rate is adopted as the final count rate, to avoid artificial corrections due to large scale 
variations of the background (especially in the large RASS fields). The count rates are given 
in Tab. O 




Figure 4.4 Comparison of measured and corrected source count rates for the extended sample 
of 106 galaxy clusters. The solid line indicates equality. 



The conversion factor for the count rate to flux conversion depends on the hydrogen 
column density, rih, on the cluster gas temperature, Tgas, on the cluster gas metallicity. A, on 
the cluster redshift, z, and on the respective detector responses for the two different PSPCs 
used. The rih value is taken as the val ue inferred from 21 cin radi o measurements for our 
galaxy at the pro jected cluster position dPickev and LockmanI 1990l : included in the EXSAS 
soft ware pa ckage. Zimmermann^ et al.| 19981: photoelectric absorption cross sections are taken 

" Il983h . 



from Morrison and McCammon 



19831 ). Gas temperatures have been estimated by compiling 
X-ray temperatures, Tx, from the literature, giving preference to temperatures measured with 
the ASCA satellite. For clusters where no ASCA measured temperature has been available, 
Tx measured with previous X-ray satellites have been used. The X-ray temperatures and 
corresponding references are given in Tab. 14.21 For two clusters included in HIFLUGCS 
no measur ed temperature ha s been found in the literature and the Lx{< 2/i^q^ Mpc) — Tx 
relation of iMarkevitrhl (jl 99HI 'I has been used. The relation for non cooling flow corrected 
luminosities and cooling flow corrected/emission weighted temperatures has been chosen, 
because the luminosities determined here have not been corrected for cooling flows and cooling 
flow corrected temperatures should be a better estimate of the cluster virial temperature. For 
/i50 = 1 this relation reads 



Tx 



Lx(< 2 Mpc) 
6.28 X 10^4 ergs s- 



1/2.09 



X 6 keV . 



(4.1) 



Since the conversion from count rate to flux depends only weakly on Tgas in the ROSAT 
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energy band for the relevant temperature range a cluster temperature k^Tg^s = 4keV has 
been assumed in a first step to determine Lx(< 2 /igQ^ Mpc) for the clusters where no gas 
temperature has been found in the literature. With this luminosity the gas temperature 
has b een estimated. Th e metallicity is set to 0.35 times the solar value for all clusters 
(e.g., Arnaud et al. 19921 . see also Chap. The redshifts have been compiled from the 



literature and are given in Tab. 14. l1 together with the corresponding references. With these 
quantities and the count rates given in Tab. 14.1 1 fluxes in the observer rest frame energy range 
0.1 — 2.4 keV have been calculate d applying a modern version of a Raymond-Smith spectral 
code ( Ravmond and SmithI 1977 ). T he results are listed in Tab. 14.11 The flux calculation 



has also been checked using XSPEC ()Arna,udlll99^ by folding the model spectrum created 
with the parameters given above with the detector response and adjusting the normalization 
to reproduce the observed count rate. It is found that for 90 % of the clusters the deviation 
between the two results for the flux measurement is less than 1 %. Luminosities in the source 
rest frame energy range 0.1 — 2.4 keV have then been calculated within XSPEC by adjusting 
the normalization to reproduce the initial flux measurements. 

The improvements of the flux determination performed here compared to the input cat- 
alogs in general are now summarized. 

1) Due to the use of a high fraction of pointed observations the photon statistics is on 
average much better, e.g, for the 33 clusters contained in REFLEX and HIFLUGCS 
one finds a mean of 841 and 19580 source photons, respectively. Consequently the 
cluster emission has been traced out to larger radii for HIFLUGCS. 

2) The higher photon statistics has allowed a proper exclusion of contaminating point 
sources (stars, active galactic nuclei (AGN), etc.) and substructure, and the separation 
of double clusters. 

3) An iterative background correction has been performed. 

4) A measured X-ray temperature is used for the flux determination in most cases. 

Simulations have shown that even for the HIFLUGCS clusters with the lowest number of 
photo ns the determined flux shows no significant trend with redshift in the relevant redshift 



range ()lkebe et al.ll200ll ). 

The parameters of the table columns of Tab. 14.11 are described as follows. Column (1) 
lists the cluster name. Names have been truncated to at most eight characters. Columns (2) 
and (3) give the equatorial coordinates of the cluster center used for the regional selection 
for the epoch J2000 in decimal degrees. Column (4) gives the heliocentric cluster redshift. 
Column (5) lists the column density of neutral galactic hydrogen in units of 10^*^ atoms cm~^. 
Column (6) gives the count rate in the channel range 52-201 which corresponds to about (the 
energy resolution of the PSPC is limited) the energy range 0.5 — 2.0 keV in units of ctss~^. 
Column (7) lists the relative (1 — o" Poissonian) error of the count rate, the flux, and the 
luminosity in percent. Column (8) gives the significance radius in h'^Q Mpc. Column (9) lists 
the flux in the energy range 0.1 — 2.4 keV in units of 10~^^ ergss"^ cm~^. Column (10) gives 
the luminosity in the energy range 0.1 — 2.4keV in units of h'^Q 10'^^ ergs s^^. Column (11) 
gives the bolometric luminosity (energy range 0.01 — 40keV) in units of /i^q^ 10^^ ergs s~^. 
Column (12) indicates whether a RASS (R) or a pointed (P) ROSAT PSPC observation has 
been used. Column (13) lists the code for the redshift reference decoded at the end of the 
table. 
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4. REDUCTION AND ANALYSIS 



Table 4.1. Cluster properties 
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Table 4.1 — Continued 
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Table 4.1 — Continued 
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Note. — The columns are described in detail at the end of Section f4.1l 
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4.2 Mass Determination 

The gravitational mass has been determined as outlined in Sect. 12.1.31 utilizing the surface 
brightness profile (Sect. I2.1.2.1|) and the gas temperature (Sects. 12.1.2.21 14.1() . Since the 
general cluster temperature structure for radii larger than about half the virial radius has 
not been studied well observationally until now, isothermality is assumed. The influence of 
the isothermal assumption on the results is discussed in Chap. [Tj 

Having determined the integrated mass as a function of radius using (|2.16|) . a physically 
meaningful fiducial radius for the mass measurement has to be defined. The radii commonly 
used are either the Abell radius, r200 5 oi' '^500 • The Abell radius is fixed at ta = 3/i5q^ Mpc. 
The radius r2oo ('^soo) is the radius within which the mean gravitational mass density is equal 
to a multiple of the critical density of the universe, (/Otot) = 200 (500) p^. It has been shown 
that a correction for red shift is not necessary for the nearby clusters included in HIFLUGCS 
( Finoguenov et al. 200 ll ) and the zero redshift value for pc (see Eqs. 12.321 12.40|) is used for 
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all calculations, i.e. pc = 4.6975 x lO^^'^g cm~^, unless noted otherwise. In Sect. IT^F^ it is 
shown that the results are not affected if instead the cluster redshift is incorporated and a 
strong evolution model is assumed. 

In order to treat clusters of different size in a homogeneous way the cluster mass is 
determined at a characteristic density but also the mass determined formally at a fixed radius 
is given for co mparison. Spherical collap se models predict a virial density (pvir) ~ 178 pc for 

= 1 (e.g., Kitavama and Sutd 1996l l. so a pragmatic approxima tion to the virial mass 
is to use r2oo as the outer boundary. Simulations performed by Evrard et al.l (|l996 ) have 
shown, however, that isothermal X-ray mass measurements may be biased towards high 
masses for r > rsoo- Furthermore for most of the clusters in HIFLUGCS (86%) up to 
rsoo no extrapolation outside the significantly detected cluster emission is necessary, i.e. 
^^500 < '''X) whereas the fraction is lower for r2oo (25 %) and ta (17%). In summary the most 
accurate results are expected for Mtot(< J^soo) = -^500i but for a comparison to predicted mass 
functions M200 is the more appropriate value (see Sect. [2?2.3.2|) . Results for all determined 
masses and their corresponding radii are given in Tab. 14.21 




350 - 



0.6 0.65 0.7 

Position /9 

Figure 4.5 Confidence contours of the fit parameters Tc and /? for the cluster A3532 as a typical 
example. Shown are the 68, 90, and 95% confidence levels (for two interesting parameters). 
The points on the 68 % error ellipse mark the value pairs used for the calculation of the 
uncertainty of the mass determination. 

A major source of uncertainty comes from the temperature measurements. However, this 
(statistical) error is less than 5 % for one third of the clusters, therefore also other sources 
of error have to be taken into account, in particular one cannot neglect the uncertainties 
of the fit parameter values when assessing the statistical errors of the mass measurements. 
Therefore mass errors have been calculated by varying the fit parameter values, /? and rc, 
along their 68 % confidence level error ellipse (Fig. 14. 5() and using the upper and lower bound 
of the quoted temperature ranges. The statistical mass error range has then been defined 



38 



4. REDUCTION AND ANALYSIS 



between the maximum and minimum mass. Note that a simple error propagation apphed 
to (|2.16|) would underestimate the uncertainty on M200 and M500, since r2oo and rsoo also 
depend on Tgas, /?, and (weakly) r^- The individual mass errors have been used in subsequent 
calculations, unless noted otherwise. In log space errors in general have been transformed as 
Alog(a;) = log(e) {x^ — x~)/{2x), where x~^ and x~ denote the upper and lower boundary 
of the quantity's error range, respectively. A mean statistical error in the gravitational mass 
estimate of 23% for clusters included in HIFLUGCS and a mean error of 27% for the 
extended sample has been found. 

The parameters of the table columns of Tab. 14.21 are described as follows. Column (1) 
lists the cluster name. Column (2) gives the (3 parameter value and the corresponding 68 % 
c.l. statistical uncertainty for two interesting parameters. Column (3) gives the core radius 
in h^Q kpc and the corresponding uncertainty. Column (4) lists the X-ray temperature along 
with its error. Columns (5) and (7) give M500 and M200 and their uncertainties in units 
of 10^^ h'^Q Mq. Columns (6) and (8) list rsoo and r2oo and their uncertainties in /i^q^ Mpc. 
Column (9) gives Ma ^ Mtot(< '^a) in units of 10^^ h'^Q Mq. Column (10) lists the code 
for the temperature reference decoded at the end of the table. Temperatures for codes 1-7 
have been determined with ASCA, code 8 with ROSAT, code 9 with EXOSAT, code 10 with 
Einstein, and code 11 with a ROSAT-ASCA Lx-Tx relation. Temperatures for code 11 are 
enclosed in parentheses and the corresponding errors have been calculated using the scatter 
in the Lx-?x relation. 

In the following the parameters of the columns of Tab. 14.31 are described. Column (1) 
lists the cluster name. Column (2) gives the central gas density in units of 10~^^gcm^'^. In 
column (3) and (4) the gas masses are given in units of 10^=^ h^Q^^ Mq. The errors include 
the uncertainty of rsoo and r2oo- The respective gas mass fractions (/gas = -^gas-M^"^) are 
listed in columns (5) and (6). The errors are determined from the (symmetrized) uncertainty 
of Mtot. 
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1 
i 


U'VT^'D A A 

ti 1 JJKA-A 


n cyo+o.oos 

U.0(0_o.o03 


OU_i 


A 

4, 


Q/-,+0.40 
•^"-0.40 


q 7C+0.58 

3-7bIo.55 


-1 

i, 


qo+O.OT 
■'5°-0.07 


r Q4+0.91 
O-y^-0.84 


z. 


■17+0.11 
■I'-O.IO 


01 

o.zl 


-1 
i 


AlUbU 


n cn7+0.040 
0.607_o 034 


94_i3 


Q 


94+0.06 
•^^-0.06 


9 cc+0.34 

2.66„o.28 


-1 

1. 


9q+0.05 
■^•J-0.04 


A 9/1+0.55 
4.24_o47 


-1 

1, 


QP-+0.08 
■yO_o.08 


fi C /I 
D.04 


z 


A 1 OCi'7 

AioD < 


n (;oc;+0.035 
0.b95_(,o32 


QQQ + 24 

■JO'J_22 


Q 


rr+0.08 
•00-0.08 


q q/i+0.36 
'J-'J^-0.32 


1 

1, 


q9+0.05 
■•5^-0.04 


c (50+0.63 
0-09_o.56 




z, 


1/1+0.08 
■14-0.07 


Q no 




z 


IVlrv VV 4 


U.44U_ooo5 


11+1 
■'■■'■-1 


1 

1, 


71 +0.09 
•'1-0.09 


0-D4_o 06 


n 
U. 


7C+O.O2 
■'"-0.03 


1 nn+010 

l-UU_o 09 


1, 


90+0.04 
■^"-0.03 


Z.ol 


Z 




n aiQ+0.038 

0.819_oo34 


/1 01 +28 
431_25 


(b. 


c-o+0.89\ 
•OO-0.78/' 


7Q+3.00 
o- ' »-2.29 


1 

±. 


Qq+0.19 
■°^-0.18 


1/1 c:9+4.92 
14.52_3 g7 


Z, 


Qq + 0.30 

■^•3-0.27 


^ A Q^ 
14. yi 


1 1 
11 


IN 01^4000 


n 4Q1+0 032 
U.'±31_o 027 


c+3 
D-2 


n 

U. 


70+0.01 
•'O-O.Ol 


n 99+0.03 
"•^^-0.02 


U. 


irq+0.02 
■O'J-0.02 


n qc:+0.04 
"-■^0-0.04 


n 
U. 


oc+0.03 
■°O-0.03 


1 OA 


A 
4 


A QCCOfi 

AoozD 


0.495_ooio 


07+5 


Q 


(.0+0.06 
•OO-O.O6 


9 qQ + 0.15 

Z.oy_o.i3 


1, 


1 0+0.02 
■l°-0.02 


Q -70+0.23 

0- '8-0.18 


1 

1, 


07+0.04 
■»'-0.03 


fi n? 
D.U / 



Z 


A 1044 


n cyn+O.lll 
u.O(a_Q 074 


Qnn+128 


/I 

4, 


7^+0.90 
• '"-0.70 


^ m+2.64 
4.iU_i 41 


1, 


49+0.26 
■^^-0.18 


(, 7q+4.54 
b.73_2.38 


Z, 


97+0.43 
■^' -0.31 


o.yo 


1 n 
lU 


A 1 ticn 
AlDOU 


n 7n/i+o.i3i 
0.704_oo8i 


001+104 


D, 


70+0.80 
• '"-0.80 


Q R9+4.91 
^^•"^-2.92 


1, 


QQ+0.28 
.»»_0 21 


15.b0^4 g5 


Q 


01 +0.45 
■01-0.35 


lo.Db 


1 
i 


A 1 1 
Algol 


n «/iq+o.oi4 
O.b43_o.oi3 


181+10 

181_io 


a 
b 


10+0.40 
•l"-0.40 


7 /(c + l.OO 
'•45_o.95 


1 

1, 


7q+0.07 
■ ''J-0.08 


1 1 ai +1.60 

11-91-1.52 


z, 


7r+0.12 
■ 'O-0.13 


1 Q m 
lo.Ul 


1 
i 


CUMA 


0.b54_oo2i 


A A +22 

344_2i 


Q 

0, 


qO + 0.34 

•OO-0.34 


11-99-1.29 


z, 


nq+0.07 
■"■J-0.08 


1 O QO + 2.08 

19.38_i 97 


Q 
0. 


99+0.11 
■^^-0.11 


1 Q ni 
lo.Ul 


z 


IN VjL>0U44 


n c;9/l+0.002 
'-'•^^^-0.003 


11+1 


1, 


07+0.01 

•07_o.oi 


r, /I 1+0.01 

0-41-0.01 


U. 


(,(,+0.01 
•DD-0.01 


cc+O.Ol 

0-o5_o.oi 


1 

1, 


o/i+O.Ol 

■04_o.oi 


1 QT 
1.0 



z 


A1736 


0.54210,1- 


"^'^-130 


3 


c;n+0.40 
•OU_o 40 


9 iq+1.23 
^•1^-0.74 


1. 


1 C-+0.18 
•10-O.lB 


3-78lftl 


1, 


07+0.34 

■**'-0.25 


6.22 


1 


A3558 


o.58oi;!:™- 


224+^ 


5, 


t-o+0.40 
•O"-0.40 


c q7+0.70 
O-'J'-0.64 


1, 


i;i;+0.07 
•OO-0.O6 


8-64tl:- 


2, 


4g+o.io 

■^"-0.10 


10.56 


1 


A3562 


U.KZ_Q 006 


99l5 


5, 


if.+0.16 
•1O-0.16 


3.68t°:^2l 


1, 


q7+0.03 
••J' -0.03 


5-83l0j« 


2, 


i(,+0.05 
■1O-0.04 


8.10 


3 


A3571 


fiis+0 010 

U.010_o 010 


1811? 


6, 


qn+0.20 
•^^"-0.20 


QQ+0.56 


1. 


7Q+0.04 
•'^-0.04 


iq qi+0.90 
±0.01_o 85 


2, 


oc+O.OO 
■°O-0.06 


14.04 


1 


A1795 


n c;Qf;+o.oo3 

U.09D_o 002 


78ti 


7, 


•°U-1.00 


9-75t?;0J 


1. 


OQ+0.12 
.09-0.14 


15.39_2 92 


2, 


qg+0.19 

■99-0.20 


15.46 


1 


A3581 


n c:4o+0.024 
U.O'±0_o.o22 


35tl 


1, 


oq+0.04 
•"''-0.04 


0.96lO:S^ 


0. 


07+0.02 
•°'-0.03 


1 (-9+0.16 
J--O^-0.13 


1. 


qo+O.OS 
■''°-0.04 


3.30 


5 


MKW8 


A ri 1+0.098 


1071™ 


3 


9Q+0.23 
•^■^-0.22 


^-l"-0.52 


1, 


14+0.13 
•1^-0.10 


q qq+1.45 
O-'J -3-0.83 


1, 


70+0.24 
■ '^-0.17 


5.60 


5 


A2029 


n c;S9+0.004 
U.OOZ_o 004 




9, 


in+1.00 
•1"-1.00 




2, 


ni+011 

■^1-0.12 


18.79lit-« 


3, 


90+0.18 
■^"-0.19 


17.62 


1 


A2052 


U.O^U_o 005 


37t^ 


3, 


nC!+0.04 
•"O-0.04 




1, 


in+o.o2 
■i"-o.oi 


3.iol»:'« 


1, 


7r+0.01 
■ 'O-0.02 


5.30 


3 


MKW3S 


n KOI +0.008 
U.Oi51_Q 007 


661^ 


3, 


70+0.20 

•'"-0.20 


O.UD_Q 3Q 


1. 


9Q+0.05 
■^^-0.04 


4.84lOo:^i 


2. 


nq+0.07 
■"'3-0.07 


7.16 


1 



40 



4. REDUCTION AND ANALYSIS 



Table 4.2— Continued 



Cluster 


/3 




Tx 


M500 




T'SOO 


M200 




T'200 


Ma 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 




(8) 


(9) 


(10) 


A 91 AO 


U.O»l_Q ggg 


1 c:4+6 
J-04_g 


y./u.i 10 


1 9Q+3-45 
±o.z,»_2 41 


9 

z. 


if,+0.17 
■J^U_o.i4 


91 n4+"-46 


Q 

0. 


01 +0.26 
■•^1-0.20 


1 Q nc; 


1 
1 


A 91 /17 




900+103 


4 Q1+0.28 
4-yi_o.28 


9 AA+0.92 
z.»»_0 g3 


1 

i 


90+0.12 

■^°-o.io 


^•°4_i.03 


9 

z. 


Aq+0.21 
■"■5-0.15 


7 91 
( .Zl 


9 
z 


A 91 fi^ 


n 7Qfi+0.030 
u. 1 »^J_o.028 


oiy_29 


1 90+0.64 
lO.Z3_0 g4 


rfi.»t)_3 74 


9 

z. 


01+0.12 
■^-L-0.11 


c;i QQ+6.96 

t)i.yy_g 13 


/I 


40+0.19 
■4y-0.18 




Q 


A 91 QQ 


U.OOO_Q 


1OA+10 

ioy_io 


4 lA+0.08 
4.iU_o.o8 


4 91+0.33 
4-^J^-0.29 


1 

1 


/1 0+0. 04 
■^•J-0.03 


a 7Q+0.52 
"•'>'-0.51 


9 

z. 


97+0.06 
■^'-0.06 


S 09 

o.yz 


9 
z 


A 99nA 


n c,q7+0 008 


b' -3 


y 91+0.25 
' -^^-0.25 


f!7+0.67 
°-"'-0.57 


1 

1 


09+0.05 
•"^-0.04 


1 q 7Q+0.96 
id./9_i 00 


9 

z. 


oQ+o.oe 

■°y-0.08 




Q 




U.DU/_o 015 


i^D_io 


7 1A+5.00 
' -^^"-2.20 


8 fi5+"-'*^ 
0-bC)_3 89 


1 

1 


09+0.59 
■°^-0.33 


1 Q 70 + 18.02 


9 

z. 


oq+0.92 
■"'^-0.52 


1 A "^"^ 




A2256 


0.914+°°t7 

— U.U4/ 


587+*? 


6.60+?,10 

— U.4U 


1 9 00+2. 38 

— J.UU 


2 


.07+°!? 

— U. il 


21.81+tS 

— 0.54 


3 


oR+o.ig 
— u..^u 


19.34 


1 


A2255 


„ +0.033 

u.(a(_o.o3o 


593li 


f. 07+0.20 

D-°'-0.20 


10.90li;i^B 


1 


Qfi+0.07 


18.65l?:67 


3 


18+0-11 
■-^°-0.09 


17.54 


3 


A3667 


"-"•^^-■--o.oos 


27911° 


(.UU_o go 


6.88ll:°« 


1 


f;o+0.08 
•O°-0.09 


111 Q+1.76 

11.13_1 gg 


2, 


■691o:m 


12.50 


1 


SllOl 


U.DOy_Q 007 


561^ 


AA+1.20 
O.UU_o 70 


9 co+1.76 
^•>JO_0.88 


1 


99+0.23 
■^^-0.16 


4^08Ii 39 


1, 


Q9+0.36 
■^^-0.25 


6.38 


9 


A2589 


0.596_o.oi2 


1181? 


3.70li:io 


3.141?:^^ 


1 


9Q+0.37 
■^y-0.22 


K A1+5.41 
5-0112.15 


2, 


nc+0.56 
■06_o.35 


7.33 


9 


A2597 


U.DJO_Q 008 


581^ 


'±•^1^-0.70 


4.52l'll? 


1 


/1 7+0.07 
-0.14 


1 /I+1.14 

7.14Ii72 


2, 


qi +0.11 

■''^-0.20 


9.27 


1 


AzDo4 


n c/in+0 051 
0.640Ioo43 


3641^^ 


q 7A+0.28 
•^^ '"-"-0.28 


1 c:+0.78 
■J-J-J_0.60 


1 


9Q+0.10 
•^^-0.09 


r qr+1.34 
0.00_i Q4 




z, 


in+0.l7 
■IOI0.I4 


(.11 


z 


A2657 


U.OOD_o.o07 


lists 


7A+0.30 


9 00+0.43 
^•°>5-0.39 


1 


9C+O.O6 
■^>^-0.()6 


A rr, + 0.68 

4-52I„fi2 


1 


AA+0.10 

■99Io.o9 


6.84 


1 


A4U38 


A +0.009 

0.541_ooo8 


59tt 


1 t;+0.03 
■J-l'J-0.03 


9 ,f,+0.09 
^•l"-0.08 


1 


14+0.02 
■J^^-0.02 


Q A 1 +0.14 

3^41_o.ii 


1, 


oA+0.03 
■OU_o.01 


5.67 


3 


A /I ncn 

A4Uoy 


n coo+o.olo 

U.OOZ_Q 010 


90tl 


4 4A+O.3O 
^•™-0.30 


Ar+0.52 
O.»0_o 4g 


1 


4Q+0.06 
■^'-'-0.06 


6.30_o.76 




z, 


99+0.09 
■^^-0.09 


CO 
O.OZ 


1 






Clusters from the extended sample not included 


in HIFLUGCS. 










A 97^A 


n R9/t+0.034 
U.UZ'±_Q 029 


212126 


(3.85ir5D 


S 4q+i-25 

•J-^-^-0.89 


1 


04+0.15 
■'J^-0.12 


c C7+1.98 
O.D^_l 4g 


9 

z. 


14+0.22 
■^4-0.21 


7 Q7 

i .y / 


1 1 


A 9ft77 
AZo / / 


n c;fifi+0.029 
U.t)bb_o.o25 


19011/7 


3.50l?;fo 


9 fil +3.32 
Z.01_i 24 


1 


99+0.39 
•^^-0.23 


A 94+5.28 


1 

±, 


Ar+0.60 
■^'-'-0.38 


fi c;7 


1 n 




r, 799+0.034 
U.(ZZ_Q 030 


241^ 


n 79+0.03 
' ^-0.02 


0.361°:°! 





/?q+0.03 


n CO+O.O8 
U-O8_o.06 


1 

J., 


■UU_o.o4 


1 7^ 


A 


AWA/T7 

7-yVV iVi 1 


n (S71 +0.027 


173+1^ 


7t;+0.09 
•J- ' '-'-(). 09 


3.791°:- 


1 


qO+0.05 

•''°-0.()4 


f\ nR+0.62 

"•"O-0.52 


9 

z. 


1Q+0.08 

■-■-^-o.oe 


o.oO 


9 
Z 


r^JZjrtOJZjU o 


U.f)4U_„ (,(,4 


641^' 


6.79l«- 


6.84i;;:- 


1 


RS+0.02 
■68-0.02 


10 SO+0-46 
iU.»U_o.4i 


9 

z. 


fifi+O-04 
■""-0.04 


1 9 90 
IZ.ZU 


9 
z 




U.bb4_o 133 


459l?«i 


(4.2ii:^:S) 


3.911?:^? 


1 


■^U-0.22 


f! 7c:+6.80 
"•'0-2.81 


9 

Z. 


97+0.60 
■^'-0.37 


Q HQ 

y.uy 


11 


9Q 


u.bUi_o.i3i 


qio+178 


O.DU_o.go 


5.68l|:l^ 


1 


t:o+0.40 
•^°-0.25 


q qr,+9.61 
y-OU_3 85 


9 

z. 


cq+0.67 
■0^-0.42 


1 1 ns 
11. Uo 


Q 

y 


A n^i^Q 


A rci +0.020 
U.0bi_o.oi8 


14811^ 


3.24lH^ 


9 00+0. 21 
Z.OO_o 19 


1 


1 0+0. 03 
■18-0.03 


7^+0.35 
■^•''*_o.34 


1 

1 . 


07+0.05 
■O'-0.06 




9 
z 


o04U 


A C41 +0.073 


130l';^f, 


(2.40l«i«) 


i.83l°:g 


1 


r,o+0.14 
■U8_o.l2 


9 Q'?+l-34 
z.»o_o 87 


1 

1, 


79+0.23 
■ ' ^-0.19 


Pi 1 Q 


1 1 
11 


A n^A^inr 
r\U04c5 W 


U.bbb_o.iii 


1981- 


(1.20lg:l?) 


U.OO_o 23 





7C+O.I6 
■'"-0.11 


1 Ofi+0-84 
^•U0_o.4i 


1 


90+0.26 
■^"5-0.19 


9 fid 


1 1 

J. J. 




U.40U_Q 013 


11811? 


1A+0.10 

■J.iU_o.io 


i.74l°:ii 


1 


^7+0.03 
■"' -0.04 


9 77+0.27 
^- ' ' -0.23 


1 

J., 


(,0+0.06 
■"°-0.05 


A 


Q 




A 081+0.619 
U.y»i_o 244 


6721^*^ 


5.001°:^° 


8.701^:?^ 


1 


09+0.51 
■"^-0.26 


1 47+18.82 
-'-^•4'-6.07 


9 

z, 


00+0-92 

■yy-0.46 


1 C 

10.00 


1 
1 




U-'4U_o 086 


142l« 


(2.58l«:;li) 


9 c-r+l.SO 
^•'-'1-0.83 


1 


90+0.21 
■^"-0.15 


4 09+2-41 
^•U^-1.33 


1 


01+0.33 
■^^1-0.23 


U-OO 


1 1 




U.bU»_o 006 


711^ 


7 91 +0.11 

' -^1^-0.11 


Q QQ + 0.35 

^•°°-0.28 


1 


00+0. 03 
•OJ-O.Ol 


1 4 1 9+0.56 
14.i2_0.53 


9 

z. 


A1+0.04 
■^^-0.04 


1 /I f^K 






u.^uu_Q oil 


203I7 


7.90l°1° 


12.501^?? 


2 


(,(,+0.12 
•UD_o.i2 


1 Q Qq+3.79 
i9.00_3 23 


Q 

0. 


94+0.21 
■^4-0.17 




1 
1 




n. ^(-9+0.217 

U. /OZ_o 123 


3441"° 


(1.181^:1?) 


0.611?,:^^ 





7r+0.15 
■ '>^-0.10 


1 1R+0-90 
i.iD_o.44 


1 
1 


9R+0.27 
■^"-0.18 


9 

z.yo 


1 1 
1 1 


A 1 A1 ^ 


u.bbu_o 015 


I79II1 


7 09+0.26 

' •■'^-0.24 


io.2ol°:g 


1 


09+0.05 
■^^-0.05 


1(S 90+1-49 
lu.zy 1 31 


Q 


Ac;+0.()9 

■"■J-o.os 




Q 




A cq9+o.007 
u.o»/;_o.oo7 


nt\ 


0.95l«:°? 


0.411°:°? 





(7(7+0.01 
■OO-O.Ol 


U.00_o.o2 


1 


04+0-02 

■"4-0.01 


1 R7 
1.0 ( 


A 

■4: 




"-"•"^^--0.030 




Q 4n+l-66 


2.89l?:«t 


1 


9-,+0.32 


^■DO_i 48 


9 

z, 


(,(,+0.61 
■""-0.23 


7 on 


Q 




A3528S 


A 4(50+0.013 
U.^DO_o 012 


lOll^ 


1 t;+0.89 
'J-l'J-0.59 


1 69+°-*'^ 


1 


0c:+0.16 
■UO-0.12 


9 7r,+ 1.39 
'"-0.80 


1 


(57+0.25 
■"' -0.19 


4.86 


8 


A3530 


„ „.,+0.114 
' ' 'J-0.085 


4211^5 


3.891°:- 




1 


/1 7+0. 15 
■^' -0.13 


' ^"^-1.80 


2, 


0(7+0.26 
■•^"-0.20 


9.82 


7 


A3532 


A ci:o+0.034 
U.UOO_o.o29 


28212I 


4.58l«:i? 


4 77+0.70 

'-0.52 


1 


40+0.07 
•4y_o.o5 


7.79lJ:i? 


2, 


00+0.12 
■"5°-0.10 


9.88 


7 


A1689 


0.6901H11 


1631^ 


Q 00+0.28 
y.zo_o_28 


10.'±»_1 00 


2 


9r.+0.06 
■^"-0.05 


24.6811:^° 


3 


rA+0.07 

■''"-o.io 


21.13 


3 


A3560 


0.566iHi 


2561^0 


(3.i6l[!:^D 




1 


14+0.12 

■-■-^-0.11 


50+1.30 
'J-'jy-0.95 


1, 


04+0.20 
■84_o.i8 


5.92 


11 


A 1775 




26011^ 


3.69l«:?? 


3.611°:™ 


1 


qc+0.06 
■■^"-0.05 


r Q1+0.83 
O.yi-0.56 


2, 


1 7+0.09 
■l'-0.07 


8.21 


3 


A1800 


0.766l°o;?^« 


3921^23 


(4.02l°;«^) 


4.75lt:«^ 


1 


4Q+0.38 
■^y-0.23 


7.97l«:?^ 


2, 


qA+0.65 
■"J^-0.37 


10.08 


11 


A1914 


A yri +0.018 
U.(01_Q 017 


23111J 


10.531^5^ 


21.4311:?^ 


2 


4fi+0.09 
•'±0-0.08 


33.99l3^^ 


3 


■°«-0.13 


26.20 


3 


NGC5813 


"-"•'""-0.103 


251^ 


(0.521H?) 


0.241°:;^ 





cc+O.ll 
■>^>^-0.07 


0.381°:- 





07+0. 16 
■8 '^-0.12 


1.32 


11 
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Table 4.2 — Continued 



Cluster 




/3 




Tx 


M500 




''500 


M200 




^200 


Ma 


Ref 


(1) 




(2) 


(3) 


(4) 


(5) 




(6) 


(7) 




(8) 


(9) 


(10) 


A2151w 





564tni^ 


681^ 


2.40«:°« 


1 r9+0.12 


1 


09+0. 03 
'-'^-0.02 


2.42l»;f3 


1 


C1 +0.03 
"^-0.04 


4.51 


3 


A3627 





rrr+0.056 
^^^-0.044 


299tlt 


fi 02+°-''** 

O-U^-O.OS 


r f;o+0.95 


1 


1:7+0.09 
^' -0.06 


9.2011:?^ 


2 


c;i+0.14 

>^J^-o.io 


11.03 


3 


TRIANGUL 





(51Q+O.OIO 
"^'-'-0.010 


279tiJ 


q an+O-ao 


13.42il:™ 


2 


in+o.09 

^^'-'-0.09 




3 


04+0.14 
■J^-0.11 


19.35 


1 


OPHIUCHU 





7/17+0.035 
-0.032 


279tg 




20.25l^:5J 


2 


41+0.10 

^-■--0.08 




3 


00+0. 16 

"■^-o.is 


25.32 


2 


ZwC11742 





71 7+0.073 
'^'-0.053 


oqo+46 
^•J^-38 


(5.23«:r3) 


6.88t;:«^ 


1 


O<5-0.18 




2 


C7+O.35 
-0.28 


12.42 


11 


A2319 





rqi +0.013 


285i;i 


8.80«-J» 


11.1611:^0 


1 


Q7+0.08 
^' -0.07 




3 


1(5+0.11 


17.17 


1 


A3695 





0*^-0.117 




(5.29«:?!) 


+5.29 
0.01 _2.i6 


1 


r7+0.39 
J' -0.24 


y.oz_3 74 


2 


rq+0.67 
>^'J-0.40 


11.12 


11 


IIZwlOS 





""^-0.097 


3651:111 


(3.44«:f8) 




1 


97+0.24 
^' -0.16 


5.04t?:«o 


2 


Q(5+0.40 
^"-0.28 


7.47 


11 


A3822 





(.oq+0.150 

"•J ="-0.093 




(4.9o«:r9) 


A q7+3.30 
-1.75 


1 


1:1+0.28 
^-■--0.20 


8.26l|:« 


2 


40+0.46 
^•J-0.34 


10.29 


11 


A3827 





QOQ+0.410 
y»y-0.192 


5931?!^ 




16.35t^^,f 


2 


9r+0.60 
^"-•-0.37 


97 4.4+29-53 
^'■**-11.46 


3 


f;9+0.99 
"^-0.60 


22.44 


11 


A3888 





Q9O+0.084 
■'^"-0.066 


40ltfo 


(8.84tl;«) 


22.00i36i» 


2 


40+0.31 
4O-0.26 


Qc; 74+15.07 
■J "J- '^-10.38 


3 


qfi+0.49 
y"J-0.44 


26.85 


11 


A3921 





7f;9+0.036 
'"^-0.030 


QOQ+26 
■J^O-23 


r 70+0.24 
^- ''J-0.23 


8.46tJi^ 


1 


on+0.08 


13.80tl:«^ 


2 


oq+0.12 

"■='-0.12 


14.37 


3 


HCG94 





r 14+0.007 

"-•J^^-o.ooe 




3.45t°:^3^ 


9 90+0. 36 
^•^"-0.34 


1 


1 7+0.06 
J^'-0.06 




1 


84 +0-09 
<54_o.09 


5.90 


6 


RXJ2344 





007+0.033 
-0.030 




(4.73to-i«) 




1 


(70+0.17 
"°-0.14 


n.27+_lit 


2 


(5q+0.27 
"''-0.25 


12.58 


11 



Note. — The columns are described in detail at the end of Section f4.2l 



References. -- fl ) iMarkevitch et al.|ll998l . ('21 iFukazawa et al.lll998l. f3l IWhitdl2000l. (4) lMatsushitalll997l (5) 
Ikebe et alJ[200ll (6) iFinoeuenov et alJl200ll (7) This work. (8) ISchindleilll996bl (9) lEdee and StewardliggiaT (10) 
David et al. I I1993I (11) Estimated from the Lx-Tx relation given bv lMarkevitchlll998l 
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4. REDUCTION AND ANALYSIS 



Tabic 4.3. Cluster properties 



Cluster 


Pe 


as(0) 


A/gas, 600 


■^gas,200 




/gas, BOO 




u 


18,200 


(1) 


(2) 




("31 








(5) 






(6) 




A0085 





.337 


16. 


,83411:?^, 


32 


q79+2,041 
•""^-2,222 


0, 


,246 


±0, 


.024 


0, 


,300 


±0 


.030 


A0119 


0, 


.026 


11. 


f.qq+0.984 
■""5 ''-(). 888 


22, 


070+1 ,373 
•°'°-l,533 


0, 


,187 


±0, 


.025 


0, 


.213 


±0, 


.029 


A0133 


0, 


.421 


5. 


' -1.089 


10, 


4+4,055 
•*J^*-2.080 


0. 


,194 


±0, 


.121 


0, 


.236 


±0, 


.148 


NGC507 


0, 


.226 


0. 


•O40_o.057 


1, 


OQr)+0.118 


0. 


,159 


±0, 


.016 


0, 


.216 


±0, 


.022 


A0262 


0, 


.158 


1, 


qoc+0.370 


4, 


1 sn+0,823 


0, 


,215 


±0, 


.022 


0, 


,293 


±0, 


.029 


A0400 


0, 


.039 


2, 


. ..+0.189 
■-'-^-'--0.177 


4, 


q4t;+0,355 


0, 


,167 


±0, 


.021 


0, 


,210 


±0, 


.028 


A0399 


0, 


.042 


18. 


rvi rr+3.889 

^J^' -3.344 


32 


•"l^*-4.984 


0, 


,180 


±0, 


.056 


0, 


,192 


±0, 


.063 


A0401 


0, 


.111 


24. 


42fi+l-716 
■^''"-1.428 


44, 


Q1 9+2-628 
■='-'■''-2.828 


0. 


,238 


±0, 


.023 


0, 


.271 


±0, 


.026 


A3112 





.544 


10, 


CI t; + 1.313 
■'J^'J-1.613 


19 


206+^'^^'' 

■^^^"-3,049 


0, 


,203 


±0 


.052 


0, 


,234 


±0 


.058 


FORNAX 





.018 


0, 


"J-" -0.032 





^90 + 0, 097 
■"^'^-0,091 


0, 


,046 


±0 


.010 


0, 


,044 


±0 


.010 


2A0335 


1, 


.066 


5, 


QOl +0,182 
■'-'^^-0.161 


9, 


Ofjc;+0,267 


0, 


,230 


±0, 


.010 


0, 


,264 


±0, 


.011 


IIIZw54 


0, 


.039 


2. 


c;qfi+0.840 

•^^'-'-0.624 


3, 


QCCC+1.199 
■="^^-0.958 


0. 


,110 


±0, 


.073 


0, 


.101 


±0, 


.069 


A3158 


0, 


.076 


13, 


/^nri+0,79,5 


23, 


00/I + 1.264 
■°O*-1.056 


0, 


,194 


±0, 


.013 


0, 


,211 


±0 


.015 


A0478 


0, 


.502 


23, 


.507l-«- 


40 


701+3,147 
•'°1-4,781 


0, 


,208 


±0, 


.042 


0, 


,228 


±0 


.047 


NGC1550 


0, 


.145 


0, 


,806irif7 


1, 


r-i r,+0,323 
•^l^"-0.249 


0, 


,117 


±0, 


.018 


0, 


,139 


±0, 


.022 


EXO0422 





.128 


3, 


•°"-0.997 


6 


990+2-216 
■^^"-1.535 


0, 


,135 


±0, 


.082 


0, 


.134 


±0, 


.082 


A3266 





.045 


26, 


■■5^0-1,780 


44 


■048l|^M 


0, 


,186 


±0 


.025 


0, 


,185 


±0 


.024 


A0496 





.651 


6, 


77c; + 0,244 
■ ' ' ■-'-(),:-iO,''j 


13 


77(1 + 0,639 
■ ' ' 'J-0,,S92 


0, 


,246 


±0 


.010 


0, 


,317 


±0 


.013 


A3376 


0, 


.020 


8, 


,873^i;?5« 


14, 


t- 40 + 1,075 
■>31O-1.109 


0, 


,140 


±0, 


.041 


0, 


,122 


±0, 


.035 


A3391 


0, 


.050 


8, 


c;77+1.306 
'-1.1B8 


16 


^99+2. 134 
■"^^-2.104 


0. 


,165 


±0, 


.038 


0, 


.198 


±0, 


.046 


A3395S 


0, 


.025 


9, 


finq+1-588 

.DUy_j 43g 


14, 


or7+2,202 
■O0'-2,118 


0, 


,109 


±0, 


.046 


0, 


.097 


±0, 


.043 


A0576 


0, 


.031 


5, 


Qo-i +1.377 
■9oJ^-1.154 


9 


f.74+1,916 
■"'*-1.962 


0, 


,112 


±0, 


.063 


0, 


,108 


±0 


.066 


A0754 


0, 


.086 


14, 


Q04+I.I88 
■^"J^-0.960 


24, 


qQ9 + 2.019 

■^»z_1.664 


0, 


,091 


±0, 


.013 


0, 


,093 


±0, 


.013 


HYDRA-A 


0, 


.634 


7, 


Q7/I +0.633 
•°' ^-0.653 


14, 


901+1.190 
■^"-'--1.085 


0. 


,209 


±0, 


.031 


0, 


.240 


±0, 


.035 


A1060 


0, 


.092 


2, 


97fl+0,257 
•^"-'-0,225 


4, 


„c:4+0.440 
•"0^-0.405 


0. 


,085 


±0, 


.010 


0, 


.096 


±0, 


.011 


A1367 


0, 


.025 


5, 


107+0,294 
-0,272 


9, 


Q-r 7+0,375 
■°1' -0,367 


0, 


,156 


±0 


.016 


0, 


,173 


±0 


.018 


MKW4 


0, 


.570 


1, 


nqq+0-097 
■'-'•^•^-0,087 


2, 


99r+0.209 
■ZZJ_0.184 


0, 


,163 


±0, 


.015 


0, 


,222 


±0, 


.021 


ZwC11215 


0, 


.052 


14. 


cor + 1-795 
, 000 775 


23 


9C-7+2.66O 
.ZO' -2.232 


0. 


,166 


±0, 


.050 


0, 


.160 


±0, 


.047 


NGC4636 


0, 


.328 


0, 


noe+0,025 


0, 


+0.0.50 
■1 ' "^-0,037 


0. 


,039 


±0, 


.004 


0, 


.050 


±0, 


.006 


A3526 


0, 


.286 


3, 


qo9+0,350 
■•J"^- 0,332 


6 


7oq+0,683 
■'°''-0.575 


0, 


,141 


±0, 


.008 


0, 


,180 


±0, 


.010 


A1644 


0, 


.043 


9, 


4(70+4,242 
■^'J0_2.921 


18 


■90711:°?! 


0, 


,231 


±0, 


.114 


0, 


.281 


±0, 


.145 


A1650 


0, 


.084 


15. 


77^+4,817 
■ ' '"-3,728 


26, 


r9q+6,834 
■'J^'5-5,777 


0. 


,164 


±0, 


.067 


0, 


.170 


±0, 


.070 


A1651 


0, 


.153 


15. 


147+1-156 
-1,209 


26, 


qi q + 1,877 

■'^l'^-2,036 


0. 


,203 


±0, 


.027 


0, 


.221 


±0, 


.029 


COMA 


0, 


.061 


21, 


/IS9+1-482 
■4°^-l,640 


38 


9Q1 +2,613 
■^°J--2,359 


0, 


,179 


±0 


.019 


0, 


,198 


±0 


.021 


NGC5044 





.672 


0, 


^9-1 +0,007 
•^^1-0,010 


0, 


or,Q+0,()21 
■O"y-0,019 


0, 


,103 


±0, 


.002 


0, 


.125 


±0, 


.003 


A1736 


0, 


.025 


5, 


Q9r+2,436 
•O^O_1.805 


13 


1 qq + 4,808 

■l'J'3_3.715 


0. 


,266 


±0, 


.119 


0, 


.347 


±0, 


.172 


A3558 


0, 


.085 


13. 


Qq9+0.988 

■OOZ_Q gg7 


26 


49J+1.840 
■^^-1-1. 720 


0. 


,258 


±0, 


.032 


0, 


.306 


±0, 


.038 


A3562 


0, 


.110 


7, 


■">^>^-0,448 


16, 


0-14+0,964 
•"l*-0,860 


0, 


,208 


±0, 


.013 


0, 


,275 


±0, 


.018 


A3571 


0, 


.144 


17, 


977+0,775 
'-0.681 


31 


r,99+l,244 
.UZZ_i 203 


0, 


,208 


±0, 


.014 


0, 


.233 


±0, 


.015 


A1795 


0, 


.499 


16. 


7gg+1.514 
■'""-1.566 


29 


.6031^:?^° 


0. 


,172 


±0, 


.035 


0, 


.192 


±0, 


.038 


A3581 


0, 


.314 


1, 


47q+0.155 
•^'•^-0.146 


2, 


70C+0.334 

■'°^-0.284 


0. 


,153 


±0, 


.015 


0, 


.184 


±0, 


.018 


MKW8 


0, 


.051 


2, 


+ 1.019 
■■J' 1-0.721 


4, 


7ei +1-806 
• '"1-1.309 


0, 


,113 


±0, 


.037 


0, 


,143 


±0 


.049 


A2029 





.564 


25, 


1 1 q+2,332 
■-'■^"'-2,322 


45, 


t-q9+3,748 
•'-'"'^-4,166 


0. 


,212 


±0, 


.037 


0, 


.242 


±0, 


.042 


A2052 


0, 


.521 


4, 


1 VR+0-201 
•l"-'-0.166 


8, 


-1 q9+0,334 
•l'^^-0.373 


0. 


,214 


±0, 


.007 


0, 


.262 


±0, 


.009 


MKW3S 


0, 


.314 


5, 


qq7+0.348 

•'^'^'-0.362 


9, 


(.^1+0.639 
■641Io.589 


0. 


,174 


±0, 


.018 


0, 


.199 


±0, 


.020 
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Table 4.3— Continued 



Cluster 


Peas (0) 


Afgas 500 


Afgas 200 




fgas 500 


/gas, 200 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 




A 91 49 


U.ZOo 


dd.0dU_3 ggg 


fin c;qs+7-866 
ou.oy»_g gi3 


u. 


.zDo ^ U.UOD 


0.288 ± 0. 


,063 


A 91 47 




7/11 7+2.132 
'•4J^'_1.537 


1U.303_3 Q92 


n 

u. 


,z^o in u.uuo 


0.350 ±0, 


,096 


A 91 f<'^ 


n 1 09 


t{f\ ^ Q/t+4-558 


ins c;fic;+6-794 
iU».0O0_g gi3 


n 
u. 


017-1-0 097 
.zi / ^ U.UZ / 


0.209 ± 0, 


,026 


A 91 QQ 


n 1 fi9 


7 c;n7+'^-^-^-^ 
/.OU<_Q 


10 anfi+0-796 

i^.»UD_Q 749 


u. 


1 7ft -1- 01 ^ 
.1 / ^ U.UlO 


0.190 ±0. 


,015 


A oor\A 


n QS7 


^d.Ddy_i 287 


4i.924_2.575 


u. 


97Q -1- 090 
.Z /O ^ U.UzU 


0.304 ± 0. 


,022 


A 9944 


n 9'^4 


1 1S fiS1+'''-5*° 
iO.D8i_4 


97 S8Q+-'2'847 
2/.»»y_7 019 


n 

u. 


,±Ol in U.-LU-L 


0.202 ± 0, 


,178 


A2256 


0.051 


— l.o2o 


35.923+?lrr 


0. 


,182 ± 0.031 


0.165 ±0, 


,029 


A2255 


0.032 


18.433ll:J«i 


31.898ltl]^ 335 


0. 


.169 ±0.016 


0.171 ±0. 


,017 


A3667 


0.066 


19.756tM°« 


40.30lli:f7| 


0. 


.287 ± 0.044 


0.360 ± 0. 


,055 


SllOl 


0.554 


4.366l,V.?*? 


7 "^1 p: + 1-732 
' •'^^^-1.163 


0. 


,169 ±0.086 


0.179 ±0, 


,091 


A o con 


0.120 




9.060l^:5f, 


0. 


,158 ±0.121 


0.181 ±0, 


,136 


A ocn'7 


0.711 


7 yi /I +0.649 
14-0.980 


12.988ll:^^«5 


0. 


,171 ±0.035 


0.182 ±0. 


,036 


A ^a.'iA 

Azoo4 


0.021 


4.0U9_Q 5gQ 


O.y00_o.898 


0. 


,143 ±0.031 


0.167 ±0. 


,037 


AZDO / 


0.098 


'±.0»D_Q 432 


».0»0_Q go7 


0, 


,173 ±0.025 


0.208 ±0, 


,030 


A A no o 


0.259 


4 064+"-i*^" 


7 797+0.396 


0, 


,188 ±0.007 


0.227 ±0, 


,008 


A4059 


0.202 


6-362I„5i2 


1 1 fi5n+"-''7-"'' 

ii.DOU_ 0.968 


0, 


,161 ± 0.020 


0.185 ±0, 


,023 




Clusters from the extended 


sample not included in HIFLUGCS. 




A OI'iA 
Az / 04 


0.063 


fi QC,1 +1 259 


1 1 791 +1-980 
11. (Zl_i 7gg 


0. 


,182 ± 0.056 


0.207 ±0, 


,063 


A 9877 


0.028 


2.574lJ;?»| 


5.072l2:|i| 


0. 


,099 ± 0.086 


0.120 ±0. 


,103 




0.204 


0.086l°:r9 


Q 101+0.019 
'-'•-'-'-'-'--0.018 


0. 


,024 ± 0.003 


0.023 ± 0. 


,003 


i\M\ iVi ( 


0.088 


t- Q7Q+0.461 
d.O 1 »_o.424 


9 lfi3+°-7i4 

y. 100-0.676 


0, 


,142 ± 0.013 


0.151 ±0, 


,014 




0.632 


i9.605i;;-:]l? 


37.15711:41? 


0, 


,287 ±0.012 


0.344 ± 0, 


,014 




0.024 


7.6381--? 


15.243llf2§ 


0. 


,195 ±0.128 


0.226 ±0, 


,161 


9Q 


0.034 


Q nTS"*"-^-^^-^ 
3.U( o_2 280 


17(^11+4.391 
1/.D11_3 g27 


0. 


,160 ±0.111 


0.189 ±0. 


,136 




0.061 


Q 7(57+0.266 
o. I U( _o.269 


'•28b_o.456 


0. 


,161 ±0.014 


0.195 ±0, 


,018 




0.078 




4.283lJ:^?« 


0, 


,134 ±0.050 


0.146 ±0, 


,055 


A n'^4S^wr 


0.019 


0.697l°;^207 


1 000+0.402 
1.000_Q 34g 


0. 


,112 ±0.064 


0.126 ±0, 


,075 


An^4j^p 


0.054 


3.090l°i|? 


6.484l°:f7i 


0. 


,177 ±0.016 


0.234 ± 0. 


,021 




0.020 


8.85312-287 


14.014l2:f« 


0, 


,102 ±0.075 


0.091 ±0, 


,073 




0.093 


2.497l»:««« 


Q Q97+1.092 
0.»Z( -0.864 


0, 


,099 ± 0.046 


0.098 ±0, 


,046 


PT^Cin74Pi 
r XVOU 1 4tO 


0.970 


24.105lJ:?i? 


42.225ll:l«? 


0. 


,271 ±0.010 


0.299 ±0, 


,012 




0.152 


17.31611:413 


27.890l|g^^ 


0. 


,139 ±0.024 


0.141 ±0. 


,025 




0.014 


0.968l°;«* 


9 09^+0-645 
Z.U/;o_Q 504 


0. 


,158 ±0.084 


0.175 ±0, 


,101 


A 1 /ll 


0.193 


i9.436ii:^;5^ 


Qo QQO+2197 
rf2.»».i_2 132 


0, 


,191 ±0.016 


0.202 ±0, 


,017 


iVi4y 
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Table 4.3— Continued 



Cluster 


Pg 
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A'/gas,200 








iOO 




.u 


IS, 200 


(1) 


(2) 




(3) 






(4) 






(5) 








(6) 




A2151w 


0, 


.159 


2, 


.300 


+0, 
-0, 


153 
.133 


4. 


.3131°: 


254 
257 


0, 


.151 


± 


0, 


.011 


0. 


.178 


± 


0. 


.013 


A3627 


0, 


.037 


10. 


.468 


+1. 
-1 


,316 
.088 


21. 
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,937 


0, 
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± 


0, 
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0. 


.231 


± 


0. 


.035 
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Note. — 



The columns are described in detail at the end of Section 



5 Temperature Structure of Abell 1835 



5.1 Motivation 

As mentioned in Sects. 12.1.2.21 and 12.1.51 knowledge about the temperature structure of the 
intracluster gas is vital for the determination of the gravitational cluster mass using the 
hydrostatic equation. Since for the aim of this work the mass is the most fundamental 
cluster property here I take advantage o f the availabilitv of data from the recently launched 
satellite observatory XMM-Newton (e.g., Jansen et al.ll200ll ^ on the galaxy cluster Abell 1835 



to study the temperature structure. 

XMM-Newton is a cornerstone mission of the European Space Agency's (ESA) Horizon 
2000 program. The most important advantage of this largest scientific satellite that has 
been launched by ESA to date, compared to previous missions is the huge eff'ective area 
(Fig. 12.1(1 . Compared to ROSAT the sensitive energy range and the spectral resolution has 
also been significantly improved and compared to ASCA the spatial re solution is much bette r 
(Tab.H^. Compared to the recently launched satellite Chandra (e.g., Weisskopf et al.lll996l ^. 



whose main virtue is the high spatial resolution of < 1 arcsec FWHM, XMM-Newton has the 
- especially for extended objects crucial - advantages of a larger field of view, and a larger 
effective area. 

Data from previous satellites have been used extensively in order to answer the question 
whether or not clusters show a trend of de creasing gas temperatu r es with increasing radius 



wnetner or not clusters snow a trend or de creasing gas temperatu r es witn mcreasmg radius 
in the outer clust er parts in ge neral (e.g., iMarkevitch et al. I I199HI : llrwin et al .ill 9991 : Iwhitfj 



200nl : llrwin and B regma nl l200(t ). The findings are contradictory. A significant temperature 



decrease would lead to an overestimate of total cluster masses if the gas is assumed to be 
isothermal (Sect. 17.2.1(1 . With the availability of cluster data obtained by XMM-Newton 
it should now be possible in principal to determine the temperature profile accurately out 
to large physical radii. The large effective area, the good simultaneous spatial and spectral 
resolution combined with its field of view make XMM-Newton the ideal instrument to resolve 
the previous discrepancies. 

Abell 1835 is an especially interesting galaxy cluster. First of all because of its medium 
redshift, z = 0.252. Only a few detailed temperature profiles for clusters around this redshift 
have been published up to now and this distance also allows to trace the cluster emission 
far out in a single observation. Secondly A1835 is believed to host one of th e most massive 
cooling flows with mass accretion rates M ~ 2 000 Mq yr'^ (jAllen et a].lll99(ih . Following one 
of the major cluster specific XMM-Newton discoveries these mass accretion rates have been 
significantly reduced recently due to high resolution spectroscop y evidence obtained w ith the 
Refiection Grating Spectrometer (RGS) onboard XMM-Newton ((Peterson et al ]l200ll l. These 



facts make A1835 one of the crucial clusters to test the predictions of cooling fiow and alter 
native models. Furthermore strong gravitational lensing has been observed f e.g.. | Allen 199^'! 



and the Sunyaev-Zeldovich effect has been measured for this cluster ^e.g.. bauskoof et al 
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Table 5.1. Observation parameters 



ID 


Date 


R.A. 


Dec. 


Duration 


Exposure 


GTI 


Detector 


Filter 


0098010101 


2000-06-28 


210.2629 


2.8507 


56 ks 


42 ks 


25 ks 


EPIC-pn FF 


Thin 



The latter measurement in combination with X-ray data has been used to constrain 



Hq, wit h the largest uncertain ty being the uncertainty in the gas temperature available pre- 
viously (jMauskonf et a].l l2nnnV With a flux /x(0.1 - 2.4keV) = 1.2 x 10"^^ ergss"^ cm'^ 
A1835 is not included in HIFLUGCS or in the extended sample but may still be called a 
fairly bright cluster. 

In this Chapter the temperature profile of A1835 is deter mined from the centr al regions 



out to large radii utilizing data from the EPIC-pn camera ()Struder et al 
XMM-Newton. 



2nnil 'l onboard 



5.2 Data Reduction and Analysis 

A1835 was observed during the performance verification phase. The data are public and have 
been retrieved from the XMM-Newton archive at MPE. The basic observation parameters 
are listed in Tab. 15.11 

Data taken with the Chandra and XMM-Newton observatories both suffer from times of 
very high background composed of 'soft protons'. In Fig. 15.11 the lightcurve of the EPIC-pn 
observation of A1835 is shown^. To obtain good statistics all events with energies above 
0.3 keV have been selected. Most of the time the count rate, Cx, is stable at about 125 
counts per 10 s. But occasionally huge flares are visible, reaching peak count rates more 
than 50 times higher than normal. Additionally smaller flares and times when the count rate 
drops to zero (counting mode times) appear. To increase the signal to noise ratio count rate 
thresholds (100 < Cx < 150 counts per 10 s) for selection of Good Time Intervals (GTI) have 
been set after inspection of the count rate histogram (Fig. 15.2']) . The improvement is obvious 
in Fig. 15.31 Dozens of sources pop out of the background in the screened image. Since the four 
corners of the EPIC-pn camera are shielded against X-ray photons, the image on the right 
hand side in Fig. 15.31 also shows that most of the remaining background is Particle Induced 
Background (PIB; background induced directly or through instrumental lines by high energy 
particles). 

For the spectral analysis all obvious point sources have been excluded as well as hot pixels 
and columns. Furthermore a correction has been applied to account for events counted during 
read out (out of time events, most obviously visible as the 'jet' extending from the cluster 
center parallel to the hot column. Fig. 15.3]) . The correction is done by creating an events 
file within the XMM SAS (Science Analysis System), where events have been redistributed 
randomly along columns keeping all other info, e.g., pattern and time^. The CTI (Charge 
Transfer Inefficiency) correction is done after the redistribution. Then a spectrum (or image) 

^The events file has been created with SAS version xmmsas_20001220_1901. 
The time information in this events file is accurate up to one read out cycle. 
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lightcurve2_10.fits 

fflSTO 




7.855e+07 7.856e+07 7.857e+07 7.858e+07 7.859e+07 



TIME [s] 

Figure 5.1 Lightcurve binned in 10 s intervals. 



ctr_hist_10_l.fits 



Idhisto 




X [count] 



Figure 5.2 Central part of the count rate [counts per 10 s] distribution. 
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Figure 5.3 EPIC-pn images of A1835 in sky coordinates in the energy band 0.35-10 keV. 
The field of view is almost 30arcmin. Left: Full 42 ks observation. Right: 25 ks Good Time 
Intervals after application of the count rate thresholds. 

created from this events file is multiplied by a mode dependent factor (= 0.06 for Full Frame 
mode, FF) and subtracted from the original spectrum (or image). Neglecting the influence of 
out of time events may lead to wrong temperature estimates since these events are assigned 
an incorrect CTI correction. Moreover, assume the cluster center has a lower temperature 
than the outer parts (as will actually be found later), then a spectrum taken from a ring of 
the outer low surface brightness region is contaminated by photons having a lower energy on 
average and the temperature estimate may be biased low. 

The subtraction of the (remaining) background is performed using spectra created in the 
same way as outlined above. Various methods may be applied for background subtraction. 
The major drawbacks of two simple methods for the observation here are described as follows. 
One possibility is to use a 'blank sky' observation like the Lockman Hole observation and 
create a spectrum from the same detector region as the source spectrum, scaling by the source 
spectrum's exposure time. The disadvantage of this method is that it ignores variation 
in the Cosmic X-ray Background (CXB). The CXB is mainly composed of A) isotropic 
extragalactic emission (dominant for photon energies E >2 keV) and B) anisotropic emission 
from galactic diffuse sources and the 'hot bubble' around the solar system (important for 
E < 2keV). Galactic photoelectric absorption (important for E < IkeV) causes an additional 
anisotropy. The CXB therefore can b e considered as constant in time but varying with 
pointing direction (at soft energies) fe.g.. llshisakilll997t ). A further drawback is that a possible 
temporal variation of the normalization and/or spectrum of the PIB is ignored. Last not 
least if this method is used one must ensure that the flare filtering criteria are equivalent for 
science and background observation. This can be achieved, e.g., by assuming a constant PIB 
and selecting an energy band containing exclusively PIB and soft proton events (possibly 
E > 12 keV) with the drawback of a significant loss in statistics. Another simple method 
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for the background subtraction is to use a ring outside the cluster emission in the science 
observation, taking advantage of the medium redshift of A1835 and XMM-Newton's field of 
view. Here the drawbacks are that the structure on the detector of the PIB is ignored. Using 
this method also the question arises whether or not one should apply a vignetting correction 
for the background spectrum. If a correction is applied then the actually not vignetted PIB 
is overestimated. If no correction is applied then the vignetted CXB is underestimated. 

More complicated met hods ar e possib le circumventing some of the above problems (e.g., 
Maierowicz and Neumann.2001 : Pratt et al. 2001. ') . Here the background from not vignetting 
corrected data in a ring where the contribution from cluster emission is negligible (7-11 arcmin 
from cluster center) has been taken because the PIB dominates the background. The main 
structural feature of the PIB is a fluorescent copper line at 8 keV, which is almost absent at 
the detector center but increases towards the outer parts of the EPIC-pn detector. Therefore 
events with 7.9 < E < 8.1 keV have been excluded from the fit. With this method the CXB 
and PIB are corrected for in one step, since the CXB is not expected to change drastically for 
this small angular distance and variations in the PIB spectrum with time have no effect since 
exactly equivalent time intervals are used. Note, however, that the absolute contribution of 
the CXB is underestimated due to the vignetting increase with increasing radius. In the outer 
parts the fitted value of the column density of neutral hydrogen, n^, drops below the galactic 
value, which indicates that leaving as a free parameter compensates for the underestimated 
CXB. However, the decrease of the nh value is certainly not an ideal correction for this effect, 
since the spectral changes invoked by both effects probably do not cancel each other exactly. 
On the other hand using a vignetting corrected background for this method, where photons 
have been corrected individually, has resulted in worse fits and underestimated temperatures 
in the outer parts because of the wrongly corrected PIB. This is due to the dominance of the 
PIB at high energies and to the steeper vignetting increase with radius for higher energies. 

Only singles, i.e. photons having deposited their energy in a single pixel on the CCD, 
have been selected for the spectral analysis. Most events are single events, but the just 
recently made available response matrix for double events will allow to increase the statistics 
at least for higher energies and test the effects of a possible dependence of the single to 
double fraction on detector position soon. A 5-cr binning for the spectra has been used to 
ensure a significant signal in each bin. Choosing a 3-0" binning changes the fit results only 
well within the statistical errors. The energy range used has nominally been set to 0.3- 
10.0 keV. In the outer low surface brightness parts the number of significant high energy bins 
decreases. The fit parameter values have therefore been determined partly within smaller 
energy ranges. If two or more different temperatures are present for a specific radial bin then 
the availability of only a weak signal at high energies may bias the estimated temperature in 
the outermost radial bin low. The presence of a lower signal to noise ratio at higher energies 
is independent of the binning method. The spec tra have been mo deled and fitted within 
XSPEC using a single temperature MeKaL model ( Mewe et al. 19951 ) including photoelectric 
absorption (wabs*mekal). The temperature, metal abundance, absorption column density, 
and normalization have been left as free parameters. Leaving the redshift as a free parameter 
yields values consistent with the optically determined redshift and with encouragingly small 
uncertainties. The 2001-05-07 release of EPIC-pn response matrices has been used. The 
spectra have not been deprojected, i.e. in the innermost radial bins also contributions from 
physically larger radii are included. However, the very steep radial surface brightness profile 
of A1835 ensures that these contaminating contributions are small. Changing the background 
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500 1000 1500 2000 

Radius [kpc] 



Figure 5.4 Temperature profiles as obtain ed with different detectors. Filled circles: EPIC-pn, 
this work; open triangles: EP IC-M0S2 (|Maierowicz and NeumannlEooH 'l: open diamonds: 
Chandra fechmidt et al. 200 ll ). Vertical error bars indicate the 90% confidence level statis- 
tical uncertainty for one interesting parameter for the first two works and the 1-a statistical 
uncertainty for the last work. Horizontal bars indicate the radial bin size. 



normalization by ±10 % has a negligible effect on the fit parameter values. Inclusion of arf 
files that correct the source spectra for telescope vignetting also affects the parameter values 
only insignificantly. 



5.3 Results and Discussion 

The best fit cluster temperatures for six radial bins are shown in Fig. 15.41 as filled circles. 
The intracluster gas in the center (r < 200 kpc) is sigi iificantly coo ler than in the outer 
parts, which may be caused either by a cooling flow (e.g., [Ribian 1994 ) or by the In terstellar 



Medium (ISM) of the central brightest cluster galaxy fe.g.. iMakishima et al.ll200ll ). Within 



the statistical uncertainty the temperature profile for r > 200 kpc can be considered isother- 
mal, in agreement with the assumption made for the mass determination throughout this 
work. Although at present a temperature decrease in the very outer part obviously cannot 
be excluded. The main limiting factor for a precise determination of the outer temperature 
is the relatively high PIB. 

Also shown in Fig. 15.41 is the temperature profile determined using data from the EPIC- 



5.3. RESULTS AND DISCUSSION 



51 



12 



> '0 

CD 

^ 8 



6 
4 







500 1000 1500 
Radius [kpc] 



2000 



Figure 5.5 Overall gas tei nperature as estiin ated with EPIC-pn (upper filled circle; this work) 
and ASCA (shaded area: lAllen et al.l[l99^ . 



M0S2 camera onboard XMM-Newton bv lMaierowicz and NeumannI tOQlh (open triangles). 
There is perfect agreement between the results from the two EPIC cameras out to ~ 600 kpc. 
For larger radii the temperature estimates are well consistent. 

The open diamonds in Fig. 15.41 additi onally indicate t he te mperature estimates obtained 



for A1835 with the Chandra satellite by ISchmidt et al.l (|200ll l. The agreement in the yery 



center (r < 100 kpc) between the Chandra and EPIC-pn results is yery good. For larger radii 
Chandra gives systematically higher temperatures, but note that 1-a errors are shown for 
Chandra. 

Using a two temperature model does not result in a significant improvement in any of the 
EPIC-pn fits, indicating the absence of evidence for gas with fcRTy below ~ 2 keV. This is in 



F ^1 — I 

agreement with the results obtained with the RGS ()Peterson et al .1120011 ^. 



To obtain a good estimate of the overall cluster temperature outside the central region, 
where the temperature is assumed to be constant - consistent with observations (Fig. 15.411 . 
a spectrum has been extracted within 250 < r < 1 500 kpc. The best fit temperature for this 
region k^Tx = 8.7 ± 0.6 keV (90% c .l.). In Fig. J_5. 51 th is result is compared to the previous 

The ASCA data suggest a higher 
. Note that the errors have been 



19961 '!. 



estimate from ASCA (shaded area; lAllen et al 
temperature (9.81^3 keV) but the results are consistent 
decreased significantly with EPIC-pn. 

Comparing this overall cluster temperature and Al835's X-ray luminosity as measured 
with a ROSAT PSPC pointed observation (rx = 2.85 Mpc) with the Lx-Tx relation of the 
part of the clusters in extended sample (Sect. that have measured temperatures one finds 
that within the scatter it falls nicely onto this relation (Fig. 15.6]) . Moreover the exceptionally 
strong central cool emission in A1835 may boost the total X-ray luminosity, which would 
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Figure 5.6 Luminosity and temperature of A1835 (open square) as compared to the Lx-Tx 
relation of 88 clusters compiled during construction of HIFLUGCS. 
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Figure 5.7 Metal abundance profile as determined with the EPIC-pn camera. 
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explain the slightly higher luminosity than expected from the relation. 

The metal abundance profile relative to solar abundances ( Anders and Grevessel 1989') is 
found to be flat within the statistical errors (Fig. 15.71) . The mean abundance is about 0.35 
solar, in agreement with the value adopted in the flux determination fSect. mT|) . 



5.4 Conclusions 

The intracluster gas in the central region is significantly cooler than in the outer part of A1835. 
The innermost circle with radius 74 kpc yields a temperature estimate Tx = 4.3ib0.2 keV. The 
gas in the outer part is consistent with being isothermal, in agreement with the assumption 
made for the mass determination throughout this work. No indications for gas temperatures 
below ~ 2keV have been found. The metal abundance profile is consistent with being flat 
with A ~ 0.35. The cluster temperature within 250 < r < 1 500 kpc has been determined as 
fce^x = 8.7 lb 0.6 keV. This temperature is almost exactly equal to the temperature of the 
innermost circle multiplied by a factor of two. However, it is possible that the temperature 
decreases further a bit for even smaller radii. The significantly reduced error of the overall 
cluster temperature will allow to put tighter constraints on Hq in combination with the SZ 
effect. Assuming the cluster temperature to stay isothermal beyond 1 500 kpc an estimate of 
the total gravitational mass has been obtained utilizing the gas density profile determined 
from a ROSAT PSPC pointed observation. Under the assumption of hydrostatic equilibrium 
it is found using (l2Tfi|l that Mtot = 1-3 ± 0.2 x 10^^ h^Q Mq within rsoo = 2.10 Mpc. 
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6 Results 



In this Chapter it is shown that a tight correlation exists between the gravitational cluster 
mass and the X-ray luminosity. This ensures that HIFLUGCS is essentially selected by 
cluster mass. Also relations between other physical properties of the clusters are shown. In 
the last part of this Chapter the cluster mass function is constructed, including the proper 
treatment of the scatter in the Lx-Mtot relation. 

6.1 Mass— Luminosity Relation 

Since one of the main aims of this work is the construction of a mass function from a flux- 
limited sample it is important to test for a correlation between X-ray luminosity and gravi- 
tational mass. In Fig. 16.11 Ly. given in the ROSAT energy band, is plotted as a function of 
M20O) showing clearly the existence of a tight (linear Pearson correlation coefficient = 0.92) 
correlation, as expected. 




Figure 6.1 Gravitational mass-X-ray luminosity relation (solid line) for the extended sample 
of 106 galaxy clusters. The dashed line gives the best fit relation for the 63 clusters included 
in HIFLUGCS (filled circles only). The bisector fit results are shown. One-o" statistical error 
bars are plotted for both axes, however, only the mass errors are larger than the symbol sizes. 

To quantify the mass-luminosity relation, a linear regression fit in log-log space has 
been performed. The method used allows for intrinsic scatter and errors in both variables 
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6. RESULTS 



Table 6.1. Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.496 


0.089 


-17.741 


1.320 


bootstrap 


1.462 


0.089 


-17.238 


1.327 


BCES(M 1 L) 


1.652 


0.085 


-20.055 


1.261 


bootstrap 


1.672 


0.086 


-20.357 


1.278 


BCES-Bisector 


1.571 


0.083 


-18.857 


1.237 


bootstrap 


1.562 


0.083 


-18.717 


1.239 


BCES-Orthogonal 


1.606 


0.086 


-19.375 


1.283 


bootstrap 


1.609 


0.088 


-19.419 


1.308 



Note. — Best fit parameter values and standard de- 
viations for the extended sample (106 clusters) for a fit 
of the form 16.111 . The rows denoted 'bootstrap' give the 
results obtained for 10 000 bootstrap resamplings. 



( Akritas and Bershadv 1996l l. Tables [^^^6 . 5 1 give the results for different fit methods, where 



minimization has been performed in vertical, horizontal, and orthogonal direction, and the bi- 
sector result is given, which bisects the best fit results of vertical and horizontal minimization. 
The fits have been performed using the form 

, / Lx(0.1-2.4keV) -N ^ / J^A ^ 



It is found, as noted in general by previous authors (e.g., Isobe et al. IQQOl l. that the chosen 



fitting method has a significant influence on the best fit parameter values^. For different 
applications different fltting methods may be required. In this work the appropriate relation 
for the application under consideration is always indicated. 

The difference between the fit results for 106 (Tab. 16. 1|) and 63 (Tab. 16. 2|) clusters may 
indicate a scale dependence of the Lx~Mtot relation, since the difference is slightly larger than 
the uncertainty evaluated with the bootstrap method. The small number of low luminosity 
clusters in HIFLUGCS compared to the extended sample may be responsible for the less 
steep relation obtained using the HIFLUGCS clusters only. Note that only two out of the 
six clusters with Lx < h'^Q 10^^ ergs s~^ are included in HIFLUGCS. To reliably detect any 
deviations from the power law shape of the Lx-Mtot relation, however, even more clusters 
with Mtot < 10^^ h^Q Mq (and possibly Mtot > 3 x 10^^ h^^ Mq) would need to be sampled. 
In Sect. l7.3l this will be further discussed. As will be seen later, in the procedure used here for 
the comparison of observed and predicted mass functions the precise shape of the Lx^Mtot 
relation is not important. 

The best fit relation for LboI (Tab. I6.3() will be discussed in Sect. 17.31 One notes that 
this relation is steeper than the one for Lx(0.1 — 2.4 keV). This is caused by the fact that 
the bolometric emission coefficient depends stronger on the gas temperature. As expected 



^This also implies that for a proper comparison of relations which have been quantified by many differe nt 
authors, e.g., the Lx-Tx relation, one and the same fitting statistic ought to be used ('e.g.. lWu et alJll99ft) . 
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Table 6.2. Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L ! M) 


1.310 


0.103 


-14.935 


1.526 


bootstrap 


1.256 


0.103 


-14.146 


1.531 


BCES(M 1 L) 


1.538 


0.105 


-18.320 


1.568 


bootstrap 


1.584 


0.113 


-18.995 


1.681 


BCES-Bisector 


1.418 


0.097 


-16.536 


1.434 


bootstrap 


1.407 


0.096 


-16.368 


1.427 


BCES-Orthogonal 


1.460 


0.105 


-17.157 


1.559 


bootstrap 


1.468 


0.110 


-17.274 


1.633 



Note. — Same as Tab. 16.11 but for the purely flux- 
limited sample (63 clusters). 



Table 6.3. Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.756 


0.091 


-21.304 


1.350 


bootstrap 


1.719 


0.090 


-20.746 


1.338 


BCES{M \ L) 


1.860 


0.084 


-22.836 


1.246 


bootstrap 


1.881 


0.084 


-23.144 


1.250 


BCES-Bisector 


1.807 


0.084 


-22.053 


1.251 


bootstrap 


1.797 


0.084 


-21.899 


1.241 


BCES-Orthogonal 


1.835 


0.085 


-22.473 


1.260 


bootstrap 


1.841 


0.085 


-22.563 


1.270 


Note. — Same as 


Tab. O but for L 


Bol- 




Table 6.4. 


Fit parameter values 




Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


1.504 


0.089 


-17.545 


1.298 


bootstrap 


1.469 


0.089 


-17.042 


1.300 


BCES(Af i L) 


1.652 


0.086 


-19.708 


1.254 


bootstrap 


1.671 


0.086 


-19.992 


1.260 


BCES-Bisector 


1.575 


0.084 


-18.590 


1.228 


bootstrap 


1.565 


0.083 


-18.445 


1.224 


BCES-Orthogonal 


1.609 


0.087 


-19.075 


1.274 


bootstrap 


1.611 


0.088 


-19.109 


1.290 



Note. — 



Same as Tab. l6.1l but for M500. 
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Table 6.5. Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(L 1 M) 


2.001 


0.124 


-25.484 


1.858 


bootstrap 


1.933 


0.126 


-24.474 


1.897 


BCES(M 1 L) 


2.488 


0.127 


-32.761 


1.902 


bootstrap 


2.519 


0.130 


-33.223 


1.940 


BCES-Bisector 


2.223 


0.118 


-28.793 


1.773 


bootstrap 


2.193 


0.119 


-28.357 


1.785 


BCES-Orthogonal 


2.404 


0.129 


-31.509 


1.938 


bootstrap 


2.415 


0.133 


-31.667 


1.984 



Note. — Same as Tab. l6.ll but for AIa for comparison 
(mass errors have not been talten into account here). 



Table 6.6. Measured scatter 



Scatter (L | M) {M \ L) Bisector 



106 clusters included in the extended sample. 



o-iog A/tot 0.21 0.21 0.21 

o-iogix 0.31 0.34 0.32 

cTiogi/M 0.17 0.18 0.17 



63 clusters included in HIFLUCCS. 



CTiogMtot 0.22 0.21 0.21 

CTiogix 0.29 0.32 0.30 

cTiogL/A/ 0.18 0.18 0.17 



Note. — Scatter measured for different rela- 
tions. 



6.2. INTRACLUSTER GAS FRACTION 



59 



the relation for M500 fTab. lH^ has a very similar slope and a higher normalization than the 
one for M200. It is given here for a possible comparison with subsequent works as well as 
the relation with the formally determined masses within an Abell radius (Tab. which is 
steeper because low luminosity clusters are assigned much higher masses. 

When constructing the mass function, the overall (measurement plus intrinsic) scatter in 
the Lx-Mtot relation may become important (Sect. IHTHl . After verifying that the scatter is 
approximately Gaussian in log space the scatter has been measured as given in Tab. 16.61 The 
scatter in Lx(0.1 — 2.4 keV), M20O) and orthogonal to the best fit line is given by cxiogLx' 
(JiogAftot) and (T\ogL/Mi respectively. The measured scatter does not vary strongly depending 
on the chosen fitting method. The slight differences are expected qualitatively. For instance 
o'logLx is smaller for (L | M) compared to (M | L) because the best fit is determined by 
minimizing the residuals in log Lx for the former. 



6.2 Intracluster Gas Fraction 

An important parameter for the understanding of cluster physics is the amount of intracluster 
gas relative to the total mass, since for instance heating and cooling processes are likely to 
affect this gas fraction. The mean gas mass fraction within r2oo for the extended sample of 

—3/2 

106 clusters has been found as (/gas) = 0.19 it 0.08 /i^q . This parameter is also important 
for cosmology as outlined in Sect. 17.6.21 In the next two Sections it is tested whether there 
are systematic variations of the gas fraction. 



6.2.1 Variation of Gas Fraction with Mass 

In the following relations between various quantities for the 106 galaxy clusters will be shown. 
Since for 18 clusters the gas temperature has been estimated using a luminosity-temperature 
relation, these clusters are omitted here and in all other Sections when the temperature is 
one of the quantities of interest. For a better illustration of possible deviations from pure 
powe r laws in some plots non parametric regression fits using a smoothed spline function 
(e.g., Simonofi 19961 . Sect. 5.6) are shown. These fits have been calculated with the standard 



statistics software package S-Plus using 5 degrees of freedom and assigning each point the 
same weight. 

In Fig. 16.21 the gas mass fraction is shown as a function of the gravitational mass. Note 
that the main uncertainty of Mgas,200 comes from the uncertainty of r2oo- Therefore the errors 
of Mgas,200 and M200 are strongly correlated and combining these errors via standard error 
propagation is not a useful option for assessing the uncertainty of /gas. Since for a given 
radius the error of the gas mass estimation is much smaller it is neglected here. Rather the 
error of /gas is determined from the error of Af2oo (including the uncertainty of ^200)- 

It is clear from Fig. 16.21 that /gas ^50^ ~ 10-30 % for most of the clusters. No clear trend 
with the gravitational mass is apparent, as long as Mtot ^ 2 x IO^^/ijq^Mq. For smaller 
masses, however, a drop of the gas fraction is indicated, even though caused by only 5-6 
clusters. A similar behavior is noted if /gas is plotted as a function of the gas mass (Fig. 16.3]) 
and of the gas temperature (Fig. 16. 4|) . Possible reasons for this sudden drop of /gas towards 
smaller systems are discussed in Sect. 17.41 
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M200 [hso"' Mq] 

Figure 6.2 Gas mass fraction as a function of gravitational mass. The solid line denotes the 
result of a smoothed spline fit, indicating a break around 2 x 10^^ h'^^ Mq. 




Figure 6.3 Gas mass fraction as a function of gas mass. The solid line denotes the result of 
a smoothed spline fit, indicating a break around 2 x 10^^ h^Q Mq. 



6.3. RELATIONS BETWEEN SHAPE PARAMETERS, TEMPERATURE, 
LUMINOSITY, AND MASS 
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Figure 6.4 Gas mass fraction as a function of gas temperature. The solid line denotes the 
result of a smoothed spline fit, indicating a break around 2keV. 

6.2.2 Variation of Gas Fraction with Radius 

To characterize the extent of the gas relative to the dark matter the ratio /gas,20o//gas,500 is 
introduced. This ratio compares the gas fraction at two characteristic radii, i.e. at two radii 
of constant overdensity. Obviously, since the overdensity decreases with increasing radius, 
/gas,20o//gas,500 > 1 indicates that the gas is more extended than the dark matter. 

In Fig. 16. 51 this ratio is plotted as a function of gravitational mass. Note that the errors of 
/gas, 200 and /gas,500 strongly correlated. Only the error in /gas,500 is therefore taken into 
account for the calculation of the uncertainty of /gas,20o//gas,500- The statistical errors still 
seem larger than the actual scatter of the points. This may be caused by the fact that the 
uncertainty in rsoo is not neglected. 

One notes that for most of the clusters the gas fraction is larger at larger radii. For 
some clusters, however, the dark matter seems to be more extended than the gas. Note that 
taking into account the error bars no cluster has /gas,20o//gas,500 significantly lower than 1. 
Interestingly, however, there appears to be a systematic trend for the larger systems that the 
gas extent becomes smaller with increasing mass relative to the dark matter extent. This 
apparent trend will be discussed further in Sect. 17.41 

6.3 Relations between Shape Parameters, Temperature, Lu- 
minosity, and Mass 

In this Section various relations between the parameters describing the surface brightness 
(and therefore gas density) profile, /? and Tc, the gas temperature, X-ray luminosity, gas and 
gravitating mass will be presented. The implications of the trends seen here will be discussed 
in Sect. 17.51 
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Figure 6.5 Gas fraction within r2oo divided by the gas fraction within rsoo as a function of 
M20Q- For most of the clusters the gas fraction increases with radius but for some of the 
high mass clusters there are indications for a decrease. The result of a smoothed spline fit is 
shown as solid line. 




Fi gure 6.6 versus core radius. Solid line: Best fit parabolic model (/3o = 0.55, Vg = 885 kpc) 
of Neumann and ArnaudI (jl999h . dashed line: same model with /3o = 0.57, rg = 775 kpc. 
Dot-dashed line: smoothed spline fit. 
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A correlation between the two fit parameters f3 and Tc as induced by the fitting process 
is indicated by the confidence contours (Fig. 14.51) . However, the typical uncertainty is much 
smaller than the range of P and rc values spanned by the sample. The origin of the corre- 
lation present in Fig. 1 6.61 may therefore be physical. Using a sample of 26 galaxy clusters 
Neumann and Arnaud (,1999. ) showed that a physical correlation exists and determined best 
fit parameters of the function 

P = Po[l + (r,/r,)2] (6.2) 

as given in the caption of Fig. 16.61 The smoothed spline fit shown in Fig. 16.61 indicates that 
this parameterization accounts well for the trend seen. However, for the 106 clusters used 
here a different set of parameter values appears to provide a better description of the data. 




Figure 6.7 Core radius versus temperature. The two parameters are clearly correlated. 

Figs. 16.71 and clear Iv show that the core radius is correlated with cluster temperature 
and mass. The temperature and core radius are determined completely independent from one 
another. Also the gravitational mass estimate is effectively independent of rc at large radii 
(r ^ Eg. 12.161) . ^ significant portion of the scatter present in these Figures is probably 
due to the presence (smaller Tc) or absence (larger Tc) of a central excess emission. First 
results of a cooling flow study of the HIFLUGCS clusters indicate a respective trend (Y. 
Chen et al., in preparation). 

In Fig. 16.91 B is plotted as a function of M20Q. Also shown is the result of a smoothed 
spline fit. For Af2oo < 3 x 10^^ h^^ Mq there is no indication of any trend of (3 with mass. 
For larger masses (3 appears to increase with increasing mass. Note that (3 is used for the 
calculation of the gravitational mass fEq. 12.1^ . The X-ray temperature determination does 
not depend on (3 (except possibly through the modeling of the surface brightness profile often 
performed to calculate the effective area for ASCA observations). Shown in Fig. 16.1(71 is (3 
versus Tx, indicating that the two parameters are almost independent of one another. For 
the hottest clusters there is, however, a weak suggestion of a trend in the same sense as in 
the previous plot. 
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Figure 6.9 The fit parameter /3 is plotted versus the gravitational mass. A smoothed spline 
fit is shown as solid line, indicating a possible trend at the high mass side. 
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1 10 
kTx [keV] 

Figure 6.10 /3 versus the intracluster gas temperature. A smoothed spUne fit is shown as solid 
Hne. No clear trends are indicated, except a very weak indication for the hottest clusters. 

The plot Lx(0.1 — 2.4keV) versus Tx has already been shown in Sect. 15.31 The correspond- 
ing best fit relations, determined using the linear regression method described in Sect. 16.11 
are given in Tab. 16.71 In Fig. 16. ill the bolometric X-ray luminosity is shown as a function of 
the gas temperature. The relation seems to be well represented by a power law of the form 
LboI tx r| (Tab. EH) over the entire temperature range. 

Fig. 16.121 shows that the X-ray luminosity and the gas mass are tightly correlated. The 
best fit relations are given in Tab. 16.91 Knowing already that there exist correlations between 
luminosity and temperature and between luminosity and gas mass it is not a big surprise 
that the gas mass correlates well with the temperature. This relation is graphically shown 
in Fig. 16.131 and tabulated in Tab. I6.1UI Note that in principal Mgas and Tx are determined 
almost independent from one another. However, since r2oo depends strongly on Tx clearly 
also the calculation of Mgas, 200 depends on Tx- 

Another relation of specific interest is the Mtot-7x relation. Its relevance derives not only 
from its importance in the study of physical processes but also from its role in the context 
of the interpretation of temperature functions (analogous to mass functions, see Sect. 16. 4j) . 
If one wants to use temperature functions to determine cosmological parameters, in some 
way the connection between temperature and mass has to be made. This may be done 
using theoretical predictions, simulations, or observations. It is therefore important to test 
the consistency of the different Mtot~Tx relations. In Fig. 16.141 the data points and the 
corresponding best fit relation (Tab. I6.11|) are shown. 

6.4 Mass Function 

In Sect. l6Tl it has been shown that the X-ray luminosity is closely correlated with cluster 
mass. Therefore the T^nax estimator (Sect. 12. 2. 3. If) can also be applied to estimate the mass 
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kTx [keV] 

Figure 6.11 Bolomctric X-ray luminosity versus temperature. The solid line represents the 
best fit result using the bisector method. 




Figure 6.12 X-ray luminosity versus gas mass in the energy band (0.1 — 2.4 keV). The solid 
line represents the best fit result using the bisector method. 
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Figure 6.13 Gas temperature versus gas mass. The solid line represents the best fit result 
using the bisector method. 




Figure 6.14 Gas temperature versus gravitational mass. The solid line represents the best fit 
result using the bisector method. 
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Table 6.7. Fit parameter values 



Fit 


a 


Aa 


A 




BCES(L 1 T) 


2.481 


0.130 


2.862 


0.093 


bootstrap 


2.367 


0.132 


2.934 


0.095 


BCES(r 1 L) 


2.724 


0.128 


2.708 


0.091 


bootstrap 


2.766 


0.131 


2.682 


0.093 


BCES-Bisector 


2.598 


0.125 


2.788 


0.089 


bootstrap 


2.552 


0.124 


2.817 


0.089 


BCES-Orthogonal 


2.692 


0.129 


2.728 


0.092 


bootstrap 


2.712 


0.132 


2.716 


0.094 



Note. — Best fit parameter values and standard 
deviations for the 88 clusters with measured tempera- 
tures for a fit of the form 

I hjTo 10"° crgss-1 / V ''^^ / 

The rows denoted 'bootstrap' give the results obtained 
for 10 000 bootstrap resamplings. 



Table 6.8. Fit parameter values 



Fit 


a 


Aa 


A 


A.4 


BCES(L 1 T) 


2.906 


0.126 


2.898 


0.090 


bootstrap 


2.776 


0.126 


2.981 


0.091 


BCES(T 1 L) 


3.067 


0.122 


2.797 


0.088 


bootstrap 


3.109 


0.125 


2.770 


0.090 


BCES-Bisector 


2.984 


0.121 


2.849 


0.087 


bootstrap 


2.933 


0.120 


2.881 


0.087 


BCES-Orthogonal 


3.050 


0.123 


2.808 


0.088 


bootstrap 


3.074 


0.125 


2.792 


0.090 



Note. — Best fit parameter values and standard 
deviations for tfic 88 clusters with measured tempera- 
tures for a fit of the form 
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Table 6.9. Fit parameter values 



Fit 


a 


Aa 


A 




BCES(I. 1 M) 


1.118 


0.034 


-11.291 


0.482 


bootstrap 


1.104 


0.033 


-11.096 


0.460 


BCES(M 1 L) 


1.162 


0.038 


-11.907 


0.533 


bootstrap 


1.184 


0.043 


-12.221 


0.606 


BCES-Bisector 


1.140 


0.035 


-11.596 


0.500 


bootstrap 


1.143 


0.036 


-11.647 


0.515 


BCES-Orthogonal 


1.143 


0.036 


-11.637 


0.513 


bootstrap 


1.149 


0.039 


-11.727 


0.546 



Note. — Best fit parameter values and standard de- 
viations for the 106 clusters of the extended sample for 
a fit of the form 

log f i2c^2.i^MiE2Vl'\ = A + a log f V 



Table 6.10. Fit parameter values 



Fit 


n 


An 


A 


A.4 


BCES(M 1 T) 


2.230 


0.106 


12.654 


0.075 


bootstrap 


2.126 


0.113 


12.721 


0.081 


BCES(T 1 M) 


2.246 


0.115 


12.644 


0.081 


bootstrap 


2.276 


0.115 


12.626 


0.081 


BCES-Bisector 


2.238 


0.108 


12.649 


0.077 


bootstrap 


2.199 


0.112 


12.675 


0.080 


BCES-Orthogonal 


2.243 


0.112 


12.645 


0.080 


bootstrap 


2.250 


0.116 


12.643 


0.082 



Note. — Best fit parameter values and standard de- 
viations for the 88 clusters with measured temperatures 

for a fit of the form 
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Table 6.11. Fit parameter values 



Fit 


a 


Aa 


A 


AA 


BCES(Af 1 T) 


1.710 


0.058 


13.735 


0.041 


bootstrap 


1.641 


0.052 


13.778 


0.037 


BCES(r 1 M) 


1.591 


0.052 


13.811 


0.036 


bootstrap 


1.633 


0.055 


13.784 


0.038 


BCES-Bisector 


1.649 


0.046 


13.774 


0.033 


bootstrap 


1.637 


0.046 


13.781 


0.033 


BCES-Orthogonal 


1.622 


0.046 


13.791 


0.032 


bootstrap 


1.635 


0.048 


13.782 


0.034 



Note. — Best fit parameter values and standard de- 
viations for tfie 88 clusters witfi measured temperatures 
for a fit of the form 



function; Imax then being a function of mass. Three different methods have been employed 
here to correct for the scatter present in the Lx-Mtot relation. If FmaxC-^x) is used instead 
of 14nax(-^tot) = ^max(-^(-^tot)), where L(Mtot) is the luminosity estimated from the Lx- 
Mtot relation using the determined cluster mass Mtot, the scatter is automatically taken 
into account. This method has be e n wid ely used in the construction of X-ray temperature 



functions, recently, e.g., bv iHenrvl (|2000f ). If Imax(-^tot) is used and the utilized Lx-Mtot 
relation is assumed to be the 'true' relation then the scatter in this relation has to be taken into 
ac count explici t ly. T heref ore following the m ethod employed for the temperature function 



by MarkevitchI (|l998l ^ and Ikebe et al. ( 2001 ) the mass function may also be estimated by 



determining V^i^^{Mtot), where the measured scatter in logLx is included. Specifically 

/oo 
y^ax(i')(2vr<^J-i/2 



I (logL^-(^ + 40)-alogMtot; I ^, 
X exp I dlogL (6.3) 

^'^logLx / 



has been used, where A and a are the best fit parameter values taken from the appropri- 
ate Lx~Mtot relation of the form ()6.H) and ciogLx i^ corresponding measured standard 
deviation in log Lx given in Tab. 16.61 

Now say that for a cluster at redshift zi to be included in the flux-limited sample it must 
have Lx > Li. Then consider the redshift shells zi — dz (shell A) and zi + dz (shell B). 
The volume enclosed by shell B is larger than the volume enclosed by shell A. This means 
that shell B contains more clusters with mass Mi which corresponds to Li. However due 
to measurement and intrinsic scatter clusters with Mi may be assigned a range of different 
luminosities Lq < Li < L2. As long as Li — Lq < L2 — Li, which is satisfied here since 
the scatter in L is roughly Gaussian in log space, this implies that more clusters with Mi 
from shell B will be included in the sample than clusters with Mi from shell A. This further 
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implies that the counted clusters with Mi have a mean L slightly larger than Li. Therefore 
the measured normalization of an L-M relation constructed with these clusters is expected 
to be slightly higher than the measured normalization of an L-M relation constructed with a 
volume-limited sample. This also implies that if the former L-M relation is used to calculate 
Vraax One gets on average a larger l^nax- This means the effect of the scatter in the L-M 
relation is already taken into account if the L-M relation determined using the flux-limited 
sample is used for the calculation of V^max(-^tot)- 

The drawback of using T^ax(-^x) is that a small number of clusters per mass bin possibly 
does not represent the true scatter well. To minimize this effect here at least ten clusters per 
mass bin are used. The drawback of using l^max(-^tot) or T4nax(-^tot), as noted bvtMarkev itchI 
(j 19981 ). is the reliance on the validity of the measured relation over the entire mass range. 
The first method and the method that accounts for the scatter explicitly (Eq. 16.3(1 have been 
tested by using Monte Carlo simulations for a precisely known Lx-Tx relation and scatter 
and have been shown to give accurate estimates of d>(T) for a large nu mber of clusters in the 
study of the HIFLUGCS temperature function bv llkebe et al. I (|2nnih . 

In Fig. 16.151 the HIFLUGCS mass functions determined by the two approaches using 
^ax(ix) and 14nax(-^tot) are shown. Each bin contains 10 clusters, apart from the highest 
mass bin, which contains 11 clusters. The highest and lowest mass clusters of the 63 clusters 
contained in the sample have been used to calculate the highest and lowest mass intervals. 
E.g., the high mass boundary of the highest mass interval is determined by adding half the 
difference between the largest and second largest mass to the second largest mass. 

As expected the method employing T^ax(-^x) prompts a mass function exhibiting a larger 
scatter, because in this case the scatter is accounted for by the actual scatter of the ten or 
eleven clusters in each mass bin. For comparison the two extreme mass functions calculated 
using T4nax(Aftot) are shown. Extreme is meant in the sense of using the steepest (A) and 
shallowest (B) Lx^Mtot relation for the HIFLUGCS sample, i.e. (M | L) with a = 1.538 and 
(L I M) with a = 1.310 (Tab.EHJ). At the low mass end (A) predicts a lower luminosity for a 
given mass than (B) resulting in a smaller T^nax and therefore a higher dn/dM. At the high 
mass side the effect is opposite resulting in a lower dn/dM for (A). The differences of these 
mass functions to the mass function calculated using Knax(-^^x) can be understood in a similar 
way and are caused partly by the indication of a deviation from a power law shape of the 
Lx-Mtot relation. Using V^iax(^^tot) results in similar mass functions as shown for the open 
symbols in Fig. 16.151 but the points lie systematically lower because the scatter is accounted 
for twice; first by using an Lx-Mtot relation derived from a flux-limited sample and secondly 
explicitly by weighting with a Gaussian function. For the comparison of the observational 
mass function to mass functions predicted by certain cosmological models T4nax(-^x) is used 
because it is independent of the precise shape of the Lx-Mtot relation and also because Lx 
has a much smaller measurement uncertainty than Mtot- The influence of the choice of the 
^max calculation on the estimation of cosmological parameters is investigated in Sect. 17.6?^ 

Fig. 16.151 shows that the mass function is decreasing fast with increasing mass, indicating 
that massive clusters are very rare objects. This fact shows an advantage of using a flux- 
limited sample for the construction of the mass function. The survey volume is largest for the 
most luminous - and therefore most massive - clusters. Thus the survey volume is largest 
for the rarest objects. Construction of a volume-limited cluster sample that samples the 
whole mass range at least equally well would require adding > 10000 clusters (for instance 
sampling 10 clusters with M200 ~ 2 x 10^^ h'^^ Mq implies sampling ~ 10 000 clusters with 
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X) 



1 1 

M200 [10" h,„-' Me] 

Figure 6.15 Gravitational mass functions for HIFLUGCS. The mass function plotted with 
filled circles has been determined using Vmax(-^x)) the ones with open symbols using 
Knaxl-^tot), where triangles correspond to the relation (L | M) and squares to (M | L) 
both for the flux-limited sample (see text Sect. I6.4() . Vertical error bars correspond to the 
formal l-cr Poisson errors, the horizontal bars indicate the mass intervals covered. 

M200 ~ 2 X 10^^ /igQ^ Mq in a volume- limited sample). 



7 Discussion 



A precise determination of distribution functions requires a high sample completeness. In 
Sect. 17.11 several completeness tests for HIFLUGCS are discussed, indicating a high com- 
pleteness. In Sect. 17.21 the cluster masses determined here are compared to independent 
determinations. The observed Lx-Mtot relation is compared to expectations in Sect. 17.31 and 
possible applications are indicated. Implications from the measured mean cluster mass to 
light ratio are outlined. Physical properties of the cluster sample are discussed in Sects. 17.41 
and l7.51 The cluster mass function is compared to previous determinations and to predictions 
of cosmological models in Sect. 17.61 Cosmological inferences from the mean cluster gas mass 
fraction are given. The total gravitational mass contained in galaxy clusters is compared to 
the total mass in the universe in Sect. 17.71 

7.1 Sample Completeness 

The sample completeness is important for the accuracy of the mass function. The selection 
criteria detailed in Sect.Elare met by 63 clusters with mean redshift {z) = 0.05 and with two 
clusters having z > 0.1. The sample is constructed from surveys with much deeper flux limits 
and high completenesses. A possible remaining incompleteness in these surveys is likely to be 
present at low fluxes close to their flux limits, which therefore would not effect HIFLUGCS. 
Nevertheless four completeness tests have been performed and are described in this Section; 
they all indicate a high completeness of HIFLUGCS. The log N — log S and Lx — z diagram 
are compared to expectations, the luminosity function is compared to luminosity functions 
of deeper surveys, and the V/Vma.x test is performed. 

Figure ITTTl shows the integral number counts as a function of X-ray flux ('log N — log 5'). 
The slope in the log — log 5 diagram is very close to the value —1.5 expected in a static 
Euclidean universe for uniformly distributed clusters. Due to the small number of clusters (4) 
the deviation is not significant for /x ^ 1 x 10~^^ ergss^^ cm^^. Since the average redshift is 
smallest for the highest fluxes large scale structure is not completely washed out at the high 
flux end, therefore the slight bump visible around /x 6 x 10^^^ ergss"^ cm~^ in Fig. 17.11 
suggests a deviation caused by cosmic variance. The effect of an expanding and finite universe 
on the log N — log S - flattening of the slope towards low fluxes - is small for the redshift 
range covered by the sample. The slope consistent with —1.5 towards the flux limit therefore 
indicates a high completeness of HIFLUGCS. 

In Fig. 17.21 the X-ray luminosity is plotted as a function of redshift. The flux limit is 
shown as a solid line^. One notes the increase in rare luminous systems with increasing 

^The correction K{Tgas, z) for converting observer rest frame luminosities to source rest frame luminosities 
(Eg. 12.521 depends on redshift and source spectrum (Tgas)- For source rest frame luminosities it is therefore 
not possible to plot the flux limit as one line in 2 dimensions (Lx, z), but rather as an area in 3 dimensions 
(Lx, z, Tgas). For consistency therefore in this 2d plot the observer rest frame luminosity is given (the correction 
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Figure 7.1 logA^(> fx) — log/x diagram. Fluxes are measured in the ROSAT energy band 
(0.1 — 2.4 keV). The dashed hne has a slope —1.5, expected for a uniform distribution of 

clusters in a static Euclidean universe ('thrcc-halves-law'), the line is normalized to produce 
the same cluster number at /x = 8 x 10~^^ ergss"^ cm~^. 




Figure 7.2 X-ray luminosity as a function of redshift. The flux limit is shown as a solid line. 
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redshift (volume). Because of the seeming underdensity of clusters in the redshift range 
0.10 < z < 0.15 a comparison with the expected number of clusters as derived from A^-body 
simulations has been performed. An Q CDM simulation, carrie d out for analysis of the power 
spectral densities of REFLEX clusters (jSchuecker et al. 1 12001 bl V adjusted to the HIFLUGCS 
survey volume in the southern hemisphere (roughly half of the total volume sampled by 
HIFLUGCS) has been used. The simulation yields 39 clusters while 33 H I FLU GCS clusters 
have been detected in this region. It is found that in fact not even one cluster with z > 0.1 is 
expected for this volume based on this simulation and the HIFLUGCS subsample also does 
not contain any cluster with a redshift larger than 0.1. This is a further piece of evidence for 
the high completeness of the sample. 

In Fig. 17.31 the HIFLUGCS X-ray luminosit y function is compared to luminosity func- 
tions of larger surve ys in the southern fREFLEX. lBohringer et al ] l2001ah and northern (BCS, 
Ebeling et al. 199?! ) hemisphere. These surveys have much deeper flux limits (Sect. and 
contain many more clusters. Very good agreement is found, which shows the high complete- 
ness and homogeneous selection of HIFLUGCS. 
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Figure 7.3 X-ray luminosity function for HIFLUGCS compared to surveys with deeper flux 
limits in the northern (N) and southern (S) hemisphere. Vertical error bars correspond to the 
formal 1-a Poisson errors (no cosmic variance), the horizontal bars indicate the luminosity 
intervals covered. 



The y/Vrr ,.^ test (e.g., iRowan-RobinsonI Il968l : ISchmidti Il96i : lAvni and Bahcalll Il98nl : 
jPeacockl Il999l . Sect. 14.5) can be used to asses a possible sample incompleteness. Assuming 
a uniform distribution of clusters a value 1/2 is expected on average. For HIFLUGCS 
{y/Vmax) = 0.47 lb 0.04, which is consistent with the expectation and this result is interpreted 
as a clear sign that HIFLUGCS covers a large enough volume for most of the Lx range to 
be representative of the local Universe with a high sample completeness. The local nature 



is less than 6 % for 90 % of the clusters anyway) . 
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of HIFLUGCS becomes obvious by noting that the result of the comoving V/Vmax test is 
almost identical to the result of the equivalent test assuming a Euclidean and non expanding 
space, i.e. ((/x//x,iim)"^/2) = 0.46. 

7.2 Comparison of Mass Determinations 

Next to a high sample completeness, reliable mass estimates are obviously important for the 
construction of the cluster mass function. In Sect. 12.131 it has been outlined that simulations 
show that X-ray mass estimates as performed in this work generally yield unbiased results 
with a relatively small scatter. In this Section masses for clusters determined using a similar 
method by different authors and using a completely independent method are directly com- 
pared to the results obtained here. Below a method for a proper comparison is described. In 
Sects. [7?2.1l and [7.2.2l the results of the comparisons with masses determined using X-ray and 
optical data, respectively, are given. 

Let us assume having two different measurement methods, A and B. For a sample of 
objects simply calculating the mean of the ratios of the results, M"^ and , of these two 
measurement methods, i.e. (-pB-/ (and clearly the ratio of the means, it7b\), is not the 



ideal way for a comparison. Consider for instance the measurement results for two objects 
(Xi,X2), = (a, l/o) and = (1,1). Both following means are always equal in this 
case and depending on a they can be arbitrarily different from 1, i.e. 



\ / M 



= 1 


a 


= 1 


> 1 


a 


> 1 


< 1 


a 


< 1 



(7.1) 



E.g. for a = 2 the former mean would imply that method A yields on average results which are 
a factor of 1.25 larger than method B. A calculation of the latter mean would imply exactly 
the opposite. Using log M* instead of Af * to calculate the above means avoids these apparently 
qualitatively contradicting results, but the results still differ quantitatively depending on 
which mean is calculated and they cannot be interpreted in a way that the values for the 
different measurement methods differ by a constant factor. Therefore here results of different 
measurements are compared using the mean and standard deviation of the differences of the 
logarithmic values, i.e. 

MD(M^, M^) = (log - log M^) , (7.2) 

which is straightforward to interpret (MD(M^,M^) = b ^ = l&M^ on average) and 
gives self consistent {MT){M^,M^) = -MD(M'^,M^)) results. 

7.2.1 X-Ray Masses 

Th e extended sample has 2 8 clusters in common with the sample analyzed independently 
by Finoguenov et al. ( 200ll ^. Independent is meant in the sense that different people have 



reduced and analyzed the data while the data themselves were obtained from the same 
satellite observatories, ROSAT and ASCA. Figure 17.41 shows a comparison of isothermal 
mass estimates within r2oo for these clusters. The differences of the logarithmic values are 
plotted as a function of the mean logarithmic mass in Fig. 17.51 There is no indication of 
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a correlation between these two quantities. Calculation of the mean difference according to 
TTIh yields MD(M^,M^) = 0.027 ± 0.113, indicating that = 1.06t|]|^M^ on average 
for masses determined within r2oo- The standard deviation is almost exactly equal to the size 
of the mean statistical mass measurement errors found in Sect. 14.21 This result shows that 
even comparing mass determinations for individual clusters not many significant differences 
are found. On average the results show very good agreement. 




[hso" Me] 



Figure 7.4 Mass determination for 28 grou ps and clusters. den otes masses determined 
in this work and masses determined bv iFinoeuenov et al.l (200l|). 



Finoguenov et al. ( 200 ll ) determined masses not only by employing the assumption of 
isothermality, M^''**°, but also using measured cluster gas temperature profiles, M^'^^'"^*^. A 
comparison for 38 clusters included in their sample yields MD{M^'^'^^'^, M^'^^°) = —0.097 ± 
0.099, indicating M^'^^'^'^ = 0.799lQ^g M^''^° on average. The influence of a possible overes- 
timation of masses determined here on the final determination of cosmological parameters is 
investigated in Sect. 17.6.21 



7.2.2 Optical Masses 

Girardi et all (Il998bl 'l determined virial masses for a sample of nearby galaxy clusters by com- 



piling optical velocity dispersions of cluster galaxies from the literature. Figure 17. HI shows a 
comparison of the mass determinations for 42 clusters common to their sample and the ex- 
tended sample presented here. No trend in the difference of the logarithmic values with mass 
is seen in Fig. 17.71 A calculation of the mean difference yields MD{M'^, M^) = 0.098 ±0.283, 
indicating that = 1.25^qq^ on average. This result shows that comparing mass 
determinations for individual clusters one may find significant differences. The average offset 
between the methods is of about the same size as the mean statistical errors found for the 
extended sample. The optical masses are about 25 % higher on average. Note, however. 
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Figure 7.5 Differences of individual mass estimates. 




Figure 7.6 Mass determination for 42 clusters . d enotes masses determined in tliis work 
and masses determined by Girardi et al. ( 1998bl l. 
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Figure 7.7 Differences of individual mass estimates. 



that since both measurements aim for determination of the virial mass, no rescahng of radii 
has been performed, i.e. masses are not determined at the same radius for each method, 
since this would underestimate the offset and standard deviation of the measurements as 
needed here. Moreover all clusters in common are compare d, i.e. critical case s like A548w, 
which is part of the larger compl ex structure of A54 8 (e.g., Davis et al. 1995l l. and A2319, 
known to host a subcluster (e.g., Oegerle et al. 1995h . have not been excluded. These two 
cases constitute the two upper extreme clusters in Fig. 17.71 The th ree lower extreme case s 
are A754, a famous cluster undergoing a major merger event (e.g., Henrv and Briell 19951 ). 
A2589, which shows no obvious indications of irregularity, except a velocity of the central cD 
galaxv offset of the c luster core velocity by about 50 % of the cluster velocity dispersion (e.g., 
Capelato et"aDll99ll ^. and A3921, which has been found to be in the state of merging (e.g.. 



Arnaud et al.lll99fih . Note also that the simple calculation of MD assigns each cluster the 



same weight, i.e. the uncertainty of the mass measurements is neglected. It is worth noting 
that (^Ifu-^ = 1-001 ± 0.787, misleadingly implying the complete absence of any systematic 
differences between the two measurement methods. The influence of a possible underesti- 
mation of masses as determined here on the determination of cosmological parameters is 
investigated in Sect. 17.6.21 Note that the isothermal mass measurements as performed here 
lie almost exactly in between the results using gas temperature gradients and using velocity 
dispersions on average. 



7.3 Mass— Luminosity Relation 

The close correlation between the X-ray luminosity and the gravitational mass found in 
Sect. I6.l1 is not surprising. If clusters form by collapse when exceeding a common overdensity 
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threshold (neglecting any redshift dependence for now), e.g., by spherical collapse fEas. 12.531 
and 12 .56^ but also any other kind of collapse from an overdensity independent of mass, then 
all clusters will share the same mean density. This picture implies that the cluster mass is 
proportional to its volume. Therefore 



Mtot oc R 



ch ' 



(7.3) 



where Rch is the virial radius, but for self similar clusters it could be any characteristic radius, 
e.g., the core radius. With the assumption of virial equilibrium one has 



Combining ()7.3() with bremsstrahlung emission. 



LboI oc />2 ri/2i23^ and Lx(0.1 - 2.4keV) oc i?; 



ch 



fSect. [Tl.2.1|) . where LboI is the bolometric luminosity, one has 

LboI oc /)2^,ri/2Mtot and Lx(0.1 - 2.4keV) oc p^^^Mt, 
Using /9gas oc Mgas-Rjof , (|7.3|1 . and (|7.4I) the relations 



LboI oc dsMtof and Lx(0.1 - 2.4 keV) oc M 



tot 



(7.4) 
(7.5) 
(7.6) 

(7.7) 



follow from (|7.(ij) . In this simple picture one therefore expects a correlation between luminosity 
and mass. For later use note that from ()7.3|) and (|7.4|) one finds 



Mtot oc r|/2 

and, combined with ()7.7|) . 

LboI oc f^^T^^ and Lx(0.1 - 2.4keV) oc f^^T^i,^ . 



(7i 



(7.9) 



O bservationally fro m the tight correlations between X-ray luminosity and te mperature 
(e.g.. lMarkevitchlll998h . and temperature and mass (e.g., Finoguenov et al. 2001 ') a correla- 
tion between luminosity and mass clearly is expe cted. Also correlations fou nd between X-ray 
luminosity and galaxy velocity dispersion (e.g., Edge and Stewari] 1991bl'l and X -rav lumi- 
nosity and mean shear strength from weak lensing studies fe.g.. lSmail et al■lll997^ indicate a 
correlation between Lx and Mtot- 

The X-rav luminosity has been comparecl direc tl y to gravitational mass e stimates bv 
Reiprich and Bohringei (ll999l.l200(]l l.lSchindlei< (ll999h.ljones and F ormanI (|l999l ^. lMiller et al, 
(|l999l V I Ettori and FabianI and Borgani and Guzzol (2001). where good correlations 

have been found in all these studies. 

In order to compare the empirical Lx-Mtot relation with predictions a quasi bolometric 
luminosity, LboIj has been calculated in the source rest frame energy range 0.01 — 40keV (for 
the relevant range of cluster gas temperatures at least 99 % of the flux is contained in this 
energy range). In Fig. 17.81 this LBonA'/2oo relation is compared to predicted relations. The 
solid line shows the best fit relation for the 106 clusters in the extended sample and the triple- 
dot-dashed line shows the best fit relation determined using HIFLUGCS. Here the bisector 
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Figure 7.8 Gravitational mass-bolometric X-ray luminosity relation compared to predicted 
relations. Shown are: best fit relation for the extended sample (solid line), best fit re- 
lation determined usi ng HIFL U G CS (trip le-dot-dashed line), self similar relation normal- 
ized by simulation s of iNavarro et al. I dSi) (dot-dashed line), pre-heated relation given by 
Evrard and Henrv (Il99lh . usin g a normalization taken from the simulations of pre-heated 
clusters bv lNavarro et alJ (|l99F 



(dashed line) . The dotted line gives the result of a smoothed 



spline fit. 



82 



7. DISCUSSION 



fit results have been used in order to trea t variables sy mmet 
method for a comparison to theory fe.g.. llsobe et ahllQ Qfll. 



/mmetrically, which is the appropriate 
The dot-dashed line shows the 



self si milar relation (LboI oc M^^^ , Eq. I7.7() normalized by the si mulations of 



Navarro et al 



and the dashed line shows the 'pre-heated' relation given bv lEvrard and Henrvl(jl99lh 
(LboI oc mAV^, see below ) , usin g a normalization taken from the simulations of pre-heated 
clusters bv lNavarro et alJ (|l995l ). The idea of pre- heating is that the intracluster gas is not 
cold initially, as in the self similar case, but is heated to a characteristic temperature by some 
form of non gravitational heat input, e.g., from SNe, early during cluster formation. Assuming 
the central regions of all clusters to have the same entropy yields the latter relationship. This 
can be seen by expressing the bolometric luminosity in terms of the central gas density, i.e. 



iBolOCPgas(0)2-^/^rl/2M( 



tot 



fe.g.. lEvrard and Henrvlll99lh . If the central specific entropy, 



3 , , Tgas 

s OC - /cr m ^- — r- 



is set to a constant value and taking 7 = 5/3 one finds 

Pgas(O) OC T|/2 . 



(7.10) 



(7.11) 



(7.12) 



Plugging this into (|7.1U|1 and using /3 = 2/3, which is close to the mean value in the extended 
sample, yields 

LBoiocrgi/,:^Mtot. (7.13) 

With H7.8() one therefore has 

(7.14) 



Fig. 17.81 shows that measured and predicted relations are in rough agreement, the differ- 
ence between the predicted relations being larger than the difference to the observed relations. 
Note, however, th at the X-rav l umino sitv is one of the most uncertain quantities to be derived 



from simulations. Frenk et al 



19991 ) recently showed in a comparison of 12 different cosmo- 
logical hydrodynamics codes that a factor of 2 uncertainty is a realistic es timate of the curren t 
accuracy. Including gas cooling in simulations worsens the situation (e.g., Balogh et al.ll200ll ). 
The slopes of the observed relations are closer to the pre-heated relation. Observationally the 
effect of pre-heating can also result in a decrease of the gas mass fraction for low temperature 
systems. This has actually been observed for the clusters in the sample (Sect. I7.4j) . The 
possibility that winds froi n SNe - originallv invoked to explain th e apparent low gas content 
of elliptical galaxies (e.g., Mathews and Baker 197l] : LarsonI 1974h - pre-heat and dilute the 
central gas a nd thereby break the self similarity has been pointed out by various authors (e.g., 
Kaiseilll98f)h. S uch a process would work most efficiently on th e least massive clusters (e.g.. 
White TTllll99l1 : iMetzier and Evrardlll 997l : IPonman et a].l ll999'V The spline fit (dotted line 
in Fig. 17. 8() indicates a weak bending, i.e. a slight deviation from a power law, over the entire 
M200 scale. There is some indication that the bend is strongest around 2 x 10^^ h7n Mg 



'0, 



suggesting a possible break. In the middle part, where most of the clusters are located, the 
slope of the dotted line is close to 4/3, while at the low mass end the slope appears even 
steeper than 11/6. 



7. 3. MASS-L UMINOSITY RELATION 



83 



The Abell cataloe; (lAbel]lll958l : lAbell et al.lll98fll ^ , as one of the first and largest systematic 
cluster catalogs, continues to be a widely used database. The Abell richness (Sect. in~T|) has 
been used frequently in the past as a select ion criterium for optical cluster samples (e.g., 



Colless and Hewett 1987 ; Katgert et al. 199f)h . The construction of the cluster mass function 



requires a selection closely related to the gravitational cluster mass. The X-ray luminosity 
has been shown above to exhibit a tight correlation with mass. In Fig . 17.91 the measured 
number of cluster member galaxies, A'gx, as taken from Abell et alJ ( 19891 ^ is compared to Lx 
as gravitational mass tracer. One notes that there exists a trend for higher mass clusters to 
have higher values of iVgx- However, there is a huge scatter. One of the main reasons for this 
poor correlation may be the much less homogeneous background compared to the X-ray case. 
It is clearly seen that a selection by X-ray luminosity is much more efficient than a selection 
by Abell richness in terms of mass. Even though only the X-ray surface brightness profile 
and neither its normalization nor the X-ray luminosity are directly used in the X-ray mass 
determination via the hydrostatic equation, it is nonetheless reassuring that a similar result 
is obtained when Ly and richness ar e compared to masses estimated from optical velocity 
dispersions ()Borgani and 
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Figure 7.9 Measured number of cluster member galaxies, A'gx, as taken from lAbell et al 
(|l989l ) and X-ray luminosity of the same (66) clusters as a function of the gravitational mass. 
Double clusters, whose components have been treated separately here, e.g., A3528n/s, are 
removed. Above the dashed line all clusters have an Abell richness R> 1. 



In Fig. l7.1Ul a plot is shown using a quantity somewhat related to Ng^. Girardi et al. ( 2000l l 
provide integrated blue band luminosities, Lsj, within the virial radius. Shown is the mass to 
light ratio, M2Q0/ Lsj, in units of /150 Mq/Lq (the mass to light ratio of the sun), as a function 
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Figure 7. 10 Mass to light rat io as a function of mass for 18 clusters in common with the 
sample of Ctirardi et al. (j2nflrih . where the integrated blue band luminosity values within the 



virial radius have been taken from. 



of M^nn for the 18 c lusters included in the extended sample used here and the C sample of 



Girardi et al. (j2nflrih . One notes that this ratio exhibits a large scatter. Furthermore a very 
w eak trend of a positiv e correlation with the mass is indicated in agreement with the findings 
of iGirardi et al. The median mass to light ratio found here equals WSh^o Mq/Lq, 

implying that the stars in clusters account only for a small fraction of the cluster mass. 
Taking into account the slight svstem atic differences in th e mass estimates (Sect. [7?2.2|) this 



value agrees with the value quoted by Girardi et al. ( 200(t ). 



One old and simple technique to estimate the mean density in the univers e is to assume 



this median mass to light ratio to be representative of the universe (see, e.g.. iBahcall et al 
Il999h . If = 1 one would expect M2oo/LBj ^ 675 hr^n M(d iLg^ based on estimates for the 



total blue band luminosity density (see refs in iGirardi et al.l I200Q 1 . The value found here 
therefore indicates ~ 0.15, independent of the Hubble constant. This value should be 
considered as a rough estimate because of the fairly small number of clusters and the scatter 
present, because of the difference in the determination of the virial radius for M200 and Lbj, 
because of the dependence on the accuracy and representativeness of the measured galaxy 
luminosity functions (which lead to the above mass to light ratio for VLj^ = 1), and because 
of the possibility that the mass to light ratio may increase or decrease on even larger scales. 
Nevertheless the estimate VL^ ~ 0.15 obtained here is in very good agreement with what will 
be found in Sect. ITTT^ when the mass function is compared to predictions. 

A wide range of possible applications becomes available with the quantification of the 
Lx-Mtot relation and its scatter. For large X-ray cluster surveys, where individual mass 
determinations are currently not feasible, luminosities can be directly converted to masses. 
No combination of observations, simulations, and theory is then needed, like the frequently 
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used approach of relating X-ray luminosities to X-ray temperatures by an observed rela- 
tion, and converting X-ray temperatures to masses using a relation where the slope is taken 
from theoretical argum ents and the normalization from hydrodynamical simulations (e.g., 
Moscardini et al. 200(1 ). The observational Lx-Mtot relation has fi rst been applied directl v 



in this sense in the power spectral analysis of REFLEX clusters ( Schuecker et al. boOlbh . 
An example of another direct application is given in Sect. 17.6.31 The Lx-Mtot relation may 
also be applied to convert theoretical or simulated mass functions to luminosity functions for 
comparison with observations, which is currently being performed in the interpretation of the 
REFLEX luminosity function. 

At this point it is important to note that even for the highest redshift cluster in the 
sample used here {z = 0.2) the dependence of the observational determination of Mtot and 
Lx on the chosen cosmological model is very weak. For instance at z = 0.2 the increase 
in the luminosity distance, D^, and the diameter distance. Da, is less than 5% going from 
(rim = 1.0, = 0) to {Qrn = 0.1, = 0). From (tOTHl one finds that Mtot(< r) oc r 
and therefore Mtot(< ''') oc D^, implying an increase of Mtot by less than 5% for the two 
models above. For Lx one has an increase of less than 10%. This means that the Lx- 
Mtot relation given here can be used unchanged for various cosmological applications (unless 
redshift ranges are probed where evolution becomes important, in this case a model dependent 
redshift correction has to be introduced). A similar calculation for T^ax shows that for the 
extreme case ^max = 0.2 the increase of T^nax is less than 14%, which is less than the size of 
the Poissonian error bars in Fig. 17.141 

More detailed investigations on the shape of the relation and the construction of a volume- 
limited sample, spanning a reasonably large range in luminosity /mass, to test how much the 
Lx-Mtot relation given here is affected by being estimated partly from a flux-limited sample 
are beyond the scope of this work but are envisaged for the future. 



7.4 Intracluster Gas Fraction 

Even though caused by only 5-6 clusters Figs. l6.2t46.4l clearlv suggest a drop in the gas fraction 
at low cluster masses/temperatures (~ 2keV). Note, however, that the X-ray emission of the 
low temperature clusters can in general be traced out to less large (compared to hotter 
clusters) radii relative to r2oo- Thus more extrapolation is performed rendering the (gas 
and total) mass estimates more uncertain. On the other hand if the gas mass fraction is 
indeed significantly lower on e would of course just expect to observe emission from the gas 
less far out. Furthermore, as lVikhlinin et 

ID CS) showed, the gas density profile steepens 



(/? increases) only slightly in the outer cluster parts and they concluded that mass estimates 
are not significantly affected. Nevertheless a possibly underestimated /? value for the low 
temperature clusters would rather lead to an overestimated value for /gas (Eq. 17.151 and 
Fig. 17.111) , which indicates that the low gas fractions for the low Tx systems found here are 
not an artefact. In any case deeper observations of low temperature clusters are needed to 
resolve this issue. Despite the problems of a high background level (Chap.El), XMM-Newton 
with its large effective area seems currently to be best suited for this purpose. 

The importance of more and especially deeper observations becom es even more obvious 
when the results from other works are compared. iDavid et al. I (jl 99,5h showed for a sample 



of 11 clusters that /gas increases continuously between ellipticals, groups, and clusters. The 
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anal ysis by Evra.rdl (Il997l ) , who combined /gas measurements from White and FabianI ( 1995h 
and Id avid et al ■I (|199,^ ^ (1 < Tx < 14keV), showed that a homogeneously analyzed cluster 
sample is vital if trends in the gas fraction are to be studied. First a trend of decreasing 
/gas with increasing Tx was indicated, but after the c luster gravitationa l masses had been 
homogeneously redetermined this trend disappeared. Allen and FabianI |l998l ) found that 
/gas is decreasing with increasing Tx for a sample of 30 high luminosity clusters after they 
had corrected the temperatures for cooling flows. Howeve r, their sample contains rather hot 
clusters mainly (5.5-26.4 keV). iRttori and FahianI (|l99 9'') found no statistically significant 
dependence of /gas on Mtot or Tx for a sample of 36 high luminosity clusters once the sample 
was corrected for a redshift dep endence^. When the reds hift dependence w as rieglecte d they 
found results in agreement with Allen and FabianI ()l998l ). lArnaud and EvrardI (jl999h found 
indications for an increase of f^as with incre asing Mtot using a sample of 24 clusters in the 



'gas 



with increasin: 



<i mill I I [ ^ 

range 2.2 < Tx < 14.6 keV. lMohr et al.l ( 19991 ) found a trend of increasing /j 
Tx fo r a sample of 45 galaxy clusters ranging in Tx from 2.4 to 10.1 keV. iRoussel et al.l 
()200nl ) analyzed 33 clusters in the temperature range 1-14 keV and found no clear trend 



a 



of /gas with Tx- iGrego et al.l (|200lh . who determined the shape of the gas density profile 
utilizing interferometric measurements of the Sunyaev-Zeldovich effect, found no indication 
for a dependence of /gas on Tx for a sample of 18 galaxy clusters. Their sample includes only 
rather hot (5.7-13.2 keV) clusters, however, whereas the effect appears to be strongest for 
cool clusters. 

It is worth noting that while David et al. ( 1995l l found that the ga s fraction comprises only 
a few percent of the total mass in ellipticals and small groups. iDavid (.1997. ). inspired by the 
early findings that a larg e fraction of the gal actic mass may be in the form of massive compact 
halo objects fMACHOs: Ulcock et ai]ll997l 'l. argued that the total baryon fraction may actu- 
ally stay constant at about 50 % independent of mass scale (fro m ellipti cals to rich clusters ) 
if MACHOs mainly consisted of stellar remnants (e.g.. Fields et a l. 1997i: Alcock et al. 19971 ). 
e.g., white dwarfs, which would also explain the abundance of heavy elements in the inter- 
galactic medium. However, t he galactic halo ma ss fraction con t ained in MACHOs has been 
significantly reduced recently JAIrork et allboonl V IPavid et all (jl 99nl 'l and lPavidl (|l997l ) also 
showed that the rati o gas mass over stellar mass increases with gas tem perature, confirmed b y 
Arnaud et al. (|l992h who found the same trend with cluster richness. iRoussel et al.l ()200nl ). 
however, did not find strong evidence of an increase of this ratio w i th Tv . The main part 
of the discrepancy seems to arise from the fact that lEous sel et al. did not find an 

increase of /gas with Tx- As mentioned later in more detail iBrvan (200Q ) compiled evidence 
for not only an increase of /gas with Tx but also for a decrease of the mass fraction contained 
in stars with increasing Tx- 

In summary it is concluded that the large spanned temperature range available in this 
work has allowed to find that the gas mass fraction A) does not vary systematically for 
clusters with Tx > 2 keV and B) drops abruptly for groups and clusters with lower gas 
temperatures in agreement with some of the previous works. The fact that the X-ray extent 
relative to the virial radius is usually smaller for groups than for clusters may be interpreted 
as a consequence of the lower gas fraction. Since proposed scenarios to explain a systematic 



■^The gas fraction depends on the diameter distance. Therefore choosing a 'wrong' cosmological model 
introduces an artificial redshift dependence for cluster samples spanning a large redshift range. Thi s effect 
may be used to constrain model parameters under the assumption of a constant gas fraction (e.g.. ISasakil 
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variation of the gas fraction with temperature or mass, e.g., pre-heating, should be considered 
also in the light of relations between other cluster parameters, they are further discussed in 
Sect. O 

In Fig. 16.51 it has been shown that for most clusters the gas mass fraction appears to vary 
with radius. For the majority of those clusters /gas seems to increase with radius, i.e. between 
rsoo and ^200- There appears to be no significant trend of this variation with system mass 
for the low to medium mass clusters. However, the gas fraction apparently becomes more 
likely to increase less fast or decrease with radius for the massive clusters. How may this 
be understood? Since /3 is an important parameter both for the gas and gravitational mass 
determination let us first estimate the behavior of /gas as a function of (3. Dividing (|2.3|1 by 
(|2.1(i)) one gets an approximate expression for the gas fraction 



/gas(< r) W /(rgas,/3gas(0),/3,rc) — 



-3/3+2 



— > 1 A/? < 1 
re 



(7.15) 



where the intracluster gas has been assumed in hydrostatic equilibrium and isothermal. 
Therefore one has by construction an increasing gas fraction with increasing radius for /3 < 2/3 
and a decreasing gas fraction for P > 2/3. This is illustrated more precisely (keeping only 
the latter two assumptions for the intracluster gas) in Fig. 17.111 for some typical parameter 
values. It is worth noting that from this Figure it is also clear that the gas fraction rises 
with radius out to > 5rc as long as /? < 0.8, i.e. for most of the clusters. Plotting the gas 
fraction ratios as determined for the clusters in this work as a function of the observed /? 
values (Fig. I7.12() confirms the decrease of /gas with increasing radius in the outer cluster 
parts for (3 > 2/3. As expected therefore /3 is the relevant parameter here. The variation 
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Figure 7.11 Radial variation of the gas fraction. Solid line: /? = 2/3, dashed line: /? = 0.8, 
dot-dashed line: (3 = 0.5. Vc = 150 kpc, fceTgas = 5keV, ne(0) = 0.01 cm~^. 



of /? as a function of mass therefore determines the variation of /gas, 200 //gas, 500 with mass. 
A quick look to Fig. 16.91 confirms this. So also trends of the radial variation of the gas mass 



88 



7. DISCUSSION 



Figure 7.12 Variation of the gas fraction as a function of /3. Error bars have been omitted. 



fraction need to be discussed in the context of other relations, which is done in the next 
Section. 



7.5 Relations between Shape Parameters, Temperature, Lu- 
minosity, and Mass 

In this Section the relations found in Sect. 16.31 are discussed. Suggested models to explain 
the observed behavior - including the variation of the gas mass fraction - are described. It 
is discussed whether the predictions of a single scenario fit all relations. 



It has long been known that a correlation exists between Lx and Tx fe.g.. iMitchell et al 

1973). Most authors find best fit relations similar to LboI ^X' ^ been found here 



too. This is significantly steeper than the simple self similar prediction LboI oc T^. Note 
that inclusion of line emission - more important for cooler systems- would further flatten 
this relation. 

In principle a dependence of the gas fraction on temperature could account for a steeper 
-^^Boi^^x relation (Eg. 17. And really it has been found in Sect. E21 that the gas fraction 
increases with increasing temperature. But this is true only for the lowest temperature 
clusters (Fig. 16. 4|) . So clearly this trend of the gas fraction cannot account for the steepness 
of the Lx-Tx relation over the entire Tx range. 

In addition to this deviation from self similarity over the entire Tx range, interestingly 
some authors have found a bend in the Lx- Tx relation in the sense that low Tx clusters give 
rise to an even steeper Lx-Tx relat ion (e.g., Helsdon and Ponma'n| [2000\ Other authors find 



no evidence for such a feature (e.g.. iMulchaev and Zabhidofjll998l ). In this work there is also 



no clear indication of any bend ranging from 0.7-13.3 keV. 

The Mtot^Tx relation has in general been found observationally to be steeper than the self 
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similar prediction Mtot oc Tgls (Eq. 17. 8|) and the normalization has been found by p revious 



authors to be lower than the normalization expected from simulations (e.g., Finoe uenov et al 



200ll . and references therein). Both trends are seen here, too, indicating that the interpreta- 
tion of temperature functions in a cosmological context depends systematically on the used 
method for the conversion between Tx and Mtot, in agreement with the findings of, e.g., 
Ikebe et al. (^001). Furthermore no clear deviation from a power law is found here. Note 
that Tgas is needed for the calculation of Aftot (Eq. I2.16() . These two quantities are therefore 
not determined independently of one another. 



A possible dependence of the gas density profile, i.e. of the shape par ameter 0, on tem^ 



perature has been indicated by observations of previous authors fe.g.. iDavid et al 



White II]|ll99ll: iMohr and Evrardlll997l: lArnaud and Evrardlll999l : Ijones and Formanlll99£ : 
T.lovd-Davies et a,l.l l2nnnl: iHelsdon and Ponmanl bnOft . H^ever. e.g.. iMohr et a,1.l (|l999lV 
Vikhlinin et al.l (1999|), and iRoussel et al.l (|2000|) have found no clear or only a weak trend. 




Fig. I6.1fll shows that for the 88 clusters with measured temperature analyzed in this work 
there is no trend found, especially not at the low temperature side. There may be indications 
of increasing /3 values with increasing Tx at the high temperature side, but the evidence is 
only weak. Stronger indications indeed are found if f3 is plotted as a function of gravitational 
mass (at the high mass side). Interpretations of this apparent trend, however, must be taken 
cautiously since /3 enters the calculation of Mtot (Eq. I2.16|) . There are two effects which 
are likely to enhance the very weak trend seen in Fig. I6.1UI as compared to Fig. 16.91 First 
of all M200 oc r^-65±o.05 i^'pg^j-,^ I6.11|) and therefore a possible trend with M200 will appear 
weaker in Tx. Secondly, the temperatures of the clusters with higher (] values at around 4~ 
6 keV (causing the slight bumb in the spline fit) are translated into relatively higher masses 

(Ea. ?nm . 

The trend of /3 at high masses may be understood by considering that the core radius 
clearly correlates with gravitational mass (Fig. 16.81 even though a large mass range is needed 
to detect this correlation due to substantial scatter). For large core radii /3 increases with 
increasing Tc (Fig. 16. 6|) . It is therefore not surprising that for large masses /? increases with 
increasing mass. 

A tight correlation exists between the gas mass and the gas temperature fFig. 16.1^ . The 
lowest temperature clusters have gas masses below the prediction by the best fit relation. 
This is consistent with the finding that the gas mass fraction drops suddenly for the lowest 
Tx clusters. A similar behavior, though less obvious, is seen in Fig. 16.121 low luminosity 
clusters have a smaller gas mass then expected from the best fit relation. 

As outlined above and in the previous Section measurements by different authors show 
different strengths of deviations from expected relations. This is obviously due to the fact 
that the measurements are difficult and several effects may lead to slight distortions of the 
results. Some of the main problems are summarized here. For low temperature clusters 
the gas is traced out least far relative to the virial radius, which affects the uncertainty of 
almost all measured quantities. Combination of heterogeneous samples, e.g., poor groups 
from author X combined with rich clusters analyzed by author Y, may introduce artificial 
trends. Purely X-ray fiux-limited samples may be less (or more) sensitive to deviations 
in relations compared to volume-limited samples (which are currently difficult to assemble 
with the necessary quality and size), especially when the X-ray luminosity is involved^. The 



3 



Consider for instance the case of a true Lx-Igas relation which has A) a larger scatter for lower Tg; 
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presence of central excess emission may bias the determination of the fit parameter values in 
the s ense that a too sm all core radius and therefore possibly a too small (3 value is measured 



(e.g., iMohr et al.lll999l 'l. This effect, if present, seems not to bias the determination of the 



gas mass, however, as has been shown by comparison of gas masses determined with single 
and double /? models (e.g., ReiprichI 1998h . Since the excess emission is often connected with 
a drop in temperature towards the very center, overall temperatures may be biased low if 
this effect is not taken into account (e.g.. Chap. El for more subtle effects on temper ature 
estimates due to accretion of small cool clumps of gas see Mathiesen and EvrardI 2001 ) . On 
the other hand the presence of substructu re may i ncreas e the temperatu re, the cor e radius and 
therefore possibly also /? (e.g., Jones an d FormanI 19991 . but seelNeuman n and A rnaud 1999). 
Mass estimates based on the hydrostatic assumption are also uncertain in major mergers and 
additionally the X-ray lu minosity may be strongly affected, at least for a short time (e.g., 
Eicker and SarazinI 200 ll ). A possible dependence of the surface brightness profile and gas 



temperature determination on the covered energy range may introduce artificial systematic 
trends, e.g., between /3 and Tx- If the gas temperature depends on radius in the outer 
parts, the assumption of isothermality leads to a biased mass estimate (Sect. 17.2.11 see also 
Chap. El). The contribution of line emission becomes non negligible for low Tx clusters (e.g., 
[Bohringcr and Hcnslcr 1982i), therefore deviations from the simplest self similar relations 
(assuming pure bremsstrahlung emission) are expected even if clusters were self similar. As 



a con sequence also possible metallicity variations between groups and clusters ([Davis et al 



I1999I for instance find lower metallicities for low Tx groups) may affect the group luminosities 
if typical cluster metallicities are assumed. Note that the deviations from the self similar 
relations are at least partly rather small, e.g., in the case of the Mtot^^'gas relation, therefore 
it is necessary to be aware even of the partially probably small effects that have just been 
outlined. However, they are not likely to strongly distort the results obtained from the 
mass function and some of the realistic possibilities have explicitly been shown to cause 
only deviations smaller than the statistical uncertainty, e.g., the possible presence of strong 
gradients in the gas temperature (Sect. [7A2|) . 

Accordingly a wide variety of models has been suggested to explain the observations. 
The currently most popular models concern heating. This heating may have oc cured uni- 



versally prior to the assemblv of groups, iust before, during or after collapse (e.g.. IWu et al 



1998a: 



Valaeeas and Si li] |l999l: lLoewensteinl '2000': Kravtsov and Yep cs '2000'; 



Tozzi and Normanf 200ll : 



Brighenti and Mathewai2001 



Wu et al 



2000; 



112001 h: .Rower et a,].H2001.V 



Since the required energy input depends on the gas density at the time of injection and the 
assumed fraction of the released energy that goes into gas heating, the cited amounts vary 
by an order of magnitude from 0.1-0.3 to 1-3 keV/particle. Various heat sources have been 

suggested. Among them SNe, AGN, and population III stars. 

But also other mechanisms have been proposed. For instance iMuanwong et al. I (|200lh find 
in simulations that cooling of the intracluster gas alone may account for obse rved deviations 
due to removal of low ent ropv gas. Based on some evidence (compiled from iMulchaev et al 



imd iHwang et all Il999l . and Icirimele et aD Il997l 'l that while the gas fraction in clusters 
decreases with dec reasing temperature the mass frac tion in stars increases (but see, e.g., 
Roussel et al. 2000, who do not find a similar trend) , iBrvanI (|2000h shows that a higher ef- 



clusters and B) a steeper slope for lower Tgas clusters. Both deviations may appear weaker then they are if a 
flux-limited sample is used because the low Lx clusters for a given Tgas have a smaller chance to be included 
in the sample than the high Lx clusters. 



7.5. RELATIONS BETWEEN SHAPE PARAMETERS, TEMPERATURE, 
LUMINOSITY, AND MASS 



91 



ficiency of galaxy formation in groups compared to clusters (e.g., David et al. 199Ci ) suffices 
to explain observations without the need of additional heating. A comparison of the predic- 
tions of this mod el to indications for a metallicity decrease in low temperature groups (e.g., 
DMi^£^^ 1999h may be interesting. 



Mathiesen et al. investigated the effects of clumping and substructure on measure- 



ments of Mgas- Using an ensemble of 48 simulated clusters they measured a mean mass 
weighted clumping factor C = (Pgas)/(/'gas)^ ~ 1.3-1.4 within rsoo- They showed that 
(-^gas'^^"'^*^/-^gas'^) ~ V^, indicating that gas masses estimated assuming a uniform den- 
sity profile at a given radius overestimate the true gas mass by about 16%, which is of the 
same order as the statistical uncertainties in this work ((AMgas,5oo/-^gas,50o) = 15%). Fur- 
thermore they found no trend with cluster bulk properties like temperature or mass. The 
possibility of scale dependent clumping therefore seems to be excluded as causing deviations 
from self similar relations. They showed further that the presence of moderate substructure 
does not bias the gas mass measurements significantly. 

A synopsis of the main trends indicated from the relations presented in this work is now 
given. The gas fraction seems to be constant for Tx ^ 2keV. For lower temperatures a steep 
decrease of /gas is indicated. A straightforward conclusion is that smaller systems contain 
comparatively less gas. At the same time there are no indications for a flattening of the gas 
density profile for these groups. Neither /3 becomes smaller, that is no flattening of the surface 
brightness profile in the outer part is observed (at least out to rx), nor does Tc become larger, 
that is no flattening in the very central regions is observed. The observed Lx-Tx and Mtot- 
Tx relations are steeper than the self similar expectation over the entire Tx range, whereas 
variations in /gas seem to occur only at low temperatures. It is therefore not necessarily clear 
that these deviations (observed at different scales) are caused by one and the same physical 
process. 

Heating of the ICG must have occured at some point as evidenced by the observed metal- 
licities. It is possible that the energy provided by SNe is sufficient to account for the decrease 
in the gas mass fraction for the lowest temperature systems, but possibl y does not suffic e to 
steepen the Lx-Tx and Mtot-Tx relations over the entire Tx range (e.g., Wu et al. 20001 ^ . A 
flatter gas density profile for low temperature clusters as one might expect in this scenario 
is not observed here. The steeper slope and lower normalization of the Mtot-Tx relation 
mav be caused by an in crease in average formation redshift for decreasing system mass (e.g.. 



Finoguenov et al.ll200lh . i.e. by hierarchical growth of structure. Since in this case low mass 



systems would have formed in a higher density environment this implies comparatively higher 
temperatures (and smaller core radii, in agreement with the steep relation between Tc and 
Tx observed here). However, neglecting cooling this effect would also increase the luminos- 
ity on average for lower temperature systems thereby tending to cause a flattening of the 
Lx-Tx relation. Since the emissivity depends on the density squared cooling may be more 
important for systems formed in a higher density environment resulting in a more efficient 
removal of gas and therefore a larger decrease in luminosity in line with observations. In this 
work no attempt has been made to model or simulate these effects in detail, but the aim has 
rather been to determine observationally relevant relations in a homogeneous way for a clus- 
ter sample not dominated by subjective selection criteria for a large temperature interval and 
describe possible scenarios accounting for them. Even though the overall physical processes 
determining the main observational appearance of galaxy clusters are well understood, at the 
moment it seems unclear which effect contributes most to the observed deviations from self 
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similarity. It appears unlikely that a single suggested effect can account for all of them. 



7.6 Mass Function 

This Section starts out with a comparison of the mass function determined here with previous 
determinations. In Sect. 17.6.21 constraints on cosmological parameters are derived. And in 
Sect. 17.6.31 the Lx-Mtot relation is applied to convert a luminosity function into a 'mass' 
function as an example. 



7.6.1 Comparison to Previous Estimates 
Bahcall and CenI |l99 



give a mass function constructed a) from optically selected clusters 
with masses dete rmined from the galaxv ric hness and b) from the cluster X-ray temperature 
function given by Henrv and ArnaudI lll991 ). Very good agreement is found for masses deter- 
mined within r\ between the Bahcall and CeiJ and the HI FLU GCS mass function (Fig. 17.13]) . 
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Figure 7.13 Cumulative gravitational mass function for HIFLUGCS using three different 
outer radii. Vertical error bars give the Poissonian errors. Horizontal bars indicate the 
individual bin sizes. Each bin contains 10 clusters, apart from the highest mass bin, which 
contains 13 clusters. The abundances from previous works are determined for Ma. 



White et al.l (ll993aD constrain the cluster abundance by using published values for the 
abundance and median velocity dispersion of richness class R > 1 Abell clusters. It is not 
surprising that their density is significantly higher than the HIFLUGCS density since they 
intentionally use conservative mass estimates, which are overestimates of the true cluster 
mass. 
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Biviano et al. and lCirardi et all (|l998al ) have determined the cluster mass function 



using optically selected cluster samples with masses determined from publishe d line-of-sight 
velocity dispersions of cluster galaxies. At the completen ess limit given by 



(|199^ 



Girardi et al 



is about a 
The latter authors 



( 1998al . triangle in Fig. I7.13|) the cluster densit y given bv Biviano e 
factor of two higher than the density given by iGirardi et al.l (jl998. 
explain this difference in density by their on average 40% smaller mass estimates due to 
an improved technique for removi ng interlopers and t he use of a surface-co rrection term 



in the virial theorem compared to iBiviano et al. I dials'). The value given by IGirardi et al 



itself lies significantly higher than the comparable HIFLUGCS density. The reason could 
lie in the fact that their optically estimated masses are on average slightly larger than the 
X-ray masses or that the external normalization for R > 1 {Ng^ > 50) clusters which they 
applied to their mass function is intrinsically higher than the normalization obtained through 
HIFLUGCS directly, or both. In Sect. 1 7. 2T^ it h as been found that for a common subsample 
of 42 clusters the virial masses determined by Girardi et al.l |l998a) are on average 25 % 
larger than the masses determined in this work. This difference might be smaller if masses 
out to the Abell radius were compared, since common radii would be used in this case. 
Increasing the X-ray masses artificially by 25 %, the diamonds in Fig. 17.131 shift towards 
higher masses, but the shift is too small to account for the difference to the triangle. The 
large scatter in the A^gx^Aftot diagram (Fig. 17. 9j) makes a reliable estimate of a best fit 
relation between these two quantities very difficult. Nevertheless, in order to get a rough 
idea of the mass within ta that corresponds to Agx = 50, I have performed fits using the 
minimization methods outlined earlier and find 5.1 < Mp^(Ngx = 50) < 8.8 x Mq. 
Note tha t this range is in agreeni ent w ith the ranges 5-8 x 10 ^^ h'^Q Mq and 5-7 x 10^^ h'^Q Mq 
given by Bahcall and GenI (|l993l ) and Girardi et al. ( 1998al ). respectively, for Ng^ = 50. This 
mass range corresponds to a cumulative number density obtained through HIFLUGCS in 
the range 1.7 < n(> Ma) < 8 .7 x 10"'^ /15Q Mpc~^. The external density estimate applied 
to normalize the IGirardi et al.l mass function therefore is a factor 1.2-6.2 higher than the 
corresponding estimate obtained here. It is therefore concluded that both effects (masses 
and normalization) are important but the latter factor is responsible for a larger fraction 
of the discrepancy. Assuming both normalizations to be determined from samples that are 
highly complete and representative of the local universe this may indicate that either X-ray 
and optical clusters are drawn from different populations or that projection effects (e.g., 
line of sight galaxy overdensities, which do not form a bound structure in three dimensions) 
po ssibly bias optically determi ned normalizations high. 

Girardi and Giurici (|200nh ently extended the mass function to loose galaxy groups, 



finding n(> 1.8 x 10^^ h^Q Mq) = 1.6-2.4 x 10"^ /ifg Mpc'^, which is outside the mass range 
that can be tested here. They find that t he group mass function can be described by a smooth 
extension of the cluster mass function hy Girardi et al.t (1998a) . Consistently this abundance 
is higher than the abund ance given by Bahcall and GenI (jl993l ) at that mass scale. 

Carlberg et al, Id 97l^ e compiled and partially reestimated abundances for local clus - 

" ^Il996h 



rg I . 

ter samples ()Henrv and Arnaudll99ll : iMazure et al.ll996l : Iwhite et ai]ll993al : like 

for comparison with higher redshift samples (the 'x' shows the density for a sample with 
higher mean redshift and therefore it cannot be compared directly). In general the compar- 



iso n to the HIFLUGCS mas s function shows be tter agreement than the abundances given 
bv IGirardi et all (|l998ah andlWhiteet all (|l993ah. 

The obvious importance of the definition of the cluster outer radius for the cluster mass 
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function can be directly appreciated by noting the large differences between the mass functions 
determined for HIFLUGCS for M500, M2Q0, and Ma in Fiar. I7.i;-il especially towards lower 
masses. Since for self similar clusters the mass scales with the third power of the characteristic 
radius (Sect. 17^ . determination of the mass within a characteristic overdensity is the natural 
choice. The formally determined Ma mass function is mainly given for the comparison with 
previous mass functions and it is recalled again that especially for the low mass systems the 
assumption of hydrostatic equilibrium may not be justified out to ta, and therefore the mass 
estimates for M500 and M200 should be considered much more precise than the estimates for 
Ma. Note that it also seems likely that especially small clusters and groups are not in virial 
equilibrium out to ta. 



7.6.2 Comparison to Predicted Mass Functions 

In this Section one of the main results of this work is derived and discussed. The mean 
density of the universe is determined by a comparison of the observed and predicted mass 
function. The formalism outlined in Sect. 12. 2. 3. 21 is used for the calculation of mass functions 
for given cosmological models. As mentioned earlier due to the low redshift range spanned 
by HIFLUGCS, the effect of a redshift correction is very small and therefore z = is set for 
all calculations, unless noted otherwise. 



Similarly to the work of llkebe et al.l (|200lf ) the statistical uncertainty in the mass deter- 



mination is incorporated in the model mass function as 

dh{M) _ 1 dn{M') 



dM ~ VmUM) y„oo dM' 
( -(logM'-logM)2 



X (2vrai,^„,,iog)-^/' exp ^ " ) dlogM' , (7.16) 




-A /tot ,log 

where ffMtot.log = 0-12 represents the logarithmic mean mass measurement uncertainty. Note 
that since HIFLUGCS is flux-limited and not volume-limited the weighting has to be per- 
formed on the mass distribution, N{M)/dM , rather than on the mass function itself. The 
effect of this weighting on the model mass function is a slight amplitude increase at the high 
mass end. 

For the modeling to be independent of the precise knowledge of the Lx"Mtot relation the 
quantitative comparison has been performed using a standard procedure on the differential 
binned mass function given in Fig. 17.141 (rather than using a maximum likelihood approach 
on the mass distribution). The values have been calculated in a fine grid of Qra and ag 
values assuming a flat cosmic geometry. The cosmological constant enters the calculation only 
through 6c and therefore has a negligible influence here. The minimum and statistical error 
ellipses for some standard confidence levels (c.l.) are given in Fig. 17.151 The tight constraints 
obtained show that with HIFLUGCS one can go beyond determining an ilm-crg relation and 
put limits on and dg individually. It is found that 

= 0.121°:°^ and as = O-QQIqH (7-17) 

(90% c.l. statistical errors for two interesting parameters), indicating a relatively low value 
for the density parameter. The large covered mass range and the specific region in Om/cs 
parameter space allow to derive these tight constraints from a local cluster sample. For 
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comparison for a given ag value the value which minimizes is calculated. These 
pairs can then roughly be described by a straight line in log space given by (using the usual 
notation) 

ffg = 0.43^1-°■3^ (7.18) 

In Fig. 17.141 also the best fit model mass functions for given = 0.5 and = 1.0 are 
plotted and one notes immediately that these value pairs give a poorer description of the 
shape of the mass function. Pre vious estimates have generally yielded a combination of 



slightly higher Q^/ag values (e.g., Peacocd 19991. Sect. 17.2). It has to be noted, however, 
that for instance in the important work of lwhite et al. ( 1993al ). who find ag « 0.57 O"'''^^, 



the authors explicitly state that their estimates of ag are probably biased high due to their 
conservative mass estimates. 
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Figure 7.14 HIFLUGCS mass function compared to the best fit model mass function with 
^^m = 0.12 and ag = 0.96 (solid line). Also shown are the best fit model mass functions for 
fixed = 0.5 (=^> (Tg = 0.56, dashed line) and Jim = 1.0 ag = 0.43, dotted line). 

Due to the large given ranges of several orders of magnitude in mass and especially 
density the values have been determined from comparison model/data naturally in log 
space. However, it has been verified that the same calculation in linear space yields best fit 
values lying within the 68 % error ellipse. 

In Sect. 16.41 arguments have been given why 14nax(-^^x) has been used for the determination 
of the mass function. Nevertheless if T4iax(-^tot) is used instead (see Fig. I6.15() . consistent 
results are obtained. Sinc e in this case one w ants to estimate L from M the relation (L | M) is 
the appropriate one re.g.. llsobe et al.lll99Ch . Performing a fit to a mass function constructed 



with 14nax(-^tot) results in best fit values Jim = 0.14 and erg = 0.86, which is consistent with 
the error range given in 1)7. 17(1 . 

As shown in Sect. 17.2.11 it is possible that the assumption of isothermality leads to an 
overestimate of the cluster masses on average. Therefore the robustness of the results has 
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been tested by multiplying the isothermal cluster masses by 0.799 and recalculating the 
minimum. As expected values for both fim and as are found to be lower. But the new 
minimum is contained well within the error ranges given in (|7.17() . On the other hand in 
Sect. 17.2.21 it has been shown by comparison with optical mass measurements that masses 
could possibly be underestimated. After increasing all cluster masses by 25 % a fit shows 
that both Sim and as become slightly larger, but again they lie well within the range (|7.17j) . 
Therefore these tests indicate that the constraints obtained here are fairly insensitive against 
systematic uncertainties in the mass measurements. 

As mentioned in Sect. 14.21 it has been shown that the small range of low redshifts covered 
here ensures that no redshift corrections need to be applied. Nevertheless it has been tested 
whether or not the best fit parameter values change if M200 is calculated using pc = Pc{z) for 
the extreme (strong evolution) case = 1, = 0), i.e. pc = 4.6975 x 10"^'^ {z + Vf gcm~^ 
for each cluster redshift (Eos. 1?!^ and 0'.4U() . The model mass function is then calculated for 
the mean redshift of HIFLUGCS, {z) = 0.05, using the formulae outlined in Sect. 12.2.3.21 It 
has been found that within the grid used here the best fit values do not change at all and also 
the error ellipses are almost not affected, thereby confirming that the application of redshift 
corrections does not affect the results. 

The value Hq = 71kms~^Mpc~^ has been adopted for the calculation of model mass 
functions b ased on the recent combination of constraints obtained using the Hubble Space 
Telescope (jlVTould et al .ll2000>). Setting the Hubble constant to their lower limit, Hq = 
65 kms~^ Mpc~^, does not affect the best fit parameter values significantly. Using even 
Hq = 60kms~^Mpc~^ changes the results for Q^i and as only well within the 68% error 
ellipse. Therefore the constraints obtained here on these cosmological parameters do not 
depend significantly on the specific choice of Hq. 

The data point in Fig. 17.141 that may be affected most by cosmic variance is the one at 
the lowest mass, since the maximum search volume is smallest for the clusters in this bin. 
Therefore the sensitivity of the best fit results on this last point has been tested by ignoring 
it. The decrease in the covered mass range of course increases the resulting error ellipse, but 
the best fit values vary only within the (smaller) 68 % error ellipse. It may be worth noting 
that leaving out the highest mass bin or leaving out the highest and lowest mass bin changes 
the best fit values only within the 90 % error ellipse. 

Since it has been found that for the estimate of the statistical errors one needs to explore 
ranges < 0.1 for Ob ~ 0.04 I regarded it necessary to check if th e approximation to 
the tr ansfer function as given in Sect. [T^.3.2l is still applicable. Recently Eisenstein and Hul 
(|l998l 'l derived a fitting function, that includes, e.g., also the oscillations induced by the 
baryons, which gives a bette r description of transfer functions computed with CMBFAST 
( Seliak and Zaldarriagal 1996 1) than fitting functions for zero or small baryon contribution to 
the total matter density derived previously. Therefore this improved version of an analytic 
transfer function has been incorporated in the procedure. It h as been found that wit hin the 
used grid the minimum does not change at all, the choi ce of the iBardeen et al. I (|l986l l fitting 
function combined with the shape parameter given by Sueivamal (|l995l l therefore seems to 
be accurate enough for the purpo ses here. However, the con fidence contours towards low Q,^ 
are getting compressed when the Eisenstein and Hul ( 1998h fitting function is used, thereby 
slightly decreasing the area of the error ellipse for a given confidence level. Since the latter 
statistical error ellipse can be regarded as more realistic, this one is shown in Fig. 17.151 

It has also been tested whether or not the recently found deviations of the PS formalism 
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Figure 7.15 Statistical confidence contours for the procedure. The cross indicates the 
minimum. Ellipses indicate 68 %, 90 %, 95 %, and 99 % confidence levels for two interesting 
parameters. 



compared to large A^-body simulations (e.g., Governato et alJ 19991 : Jenkins et ahlbnOll ) have 
a significant infiuence on the results obtained here. The best fit PS model (fim = 0.12, 
CTg = 0.96) has been coi npared to the mode l obtained using the 'universal' mass function (fit 
to A^-body simulations, Jenkins et ahlEoOll l for the same parameter values. These two models 
agree weh for M < 10^^ h^^ Mq. The differences become larger than the s ize of the Poiss onian 



error bars (Fig. [71^ for M > 2 x 10^^ h^^ Mr. 



0) 



in the sense that the iJenkins et al 



mass 



function predicts higher cluster abundances than PS. For larger values of the differences 
become comparable to the size of the error bars at lower masses, e.g., for Jim = fTg = 0.5 
around M ~ 5 x 10^^ h'^Q Mp). To estimate the influence of these differences on the best fit 



values derived using the PS mass function, the parameter values of the IJenkins et al.l model 
have been adjusted to reproduce the PS mass function, where = 0.15 and erg = 0.86 
have been found. The value for Vt^ becomes slightly larger but the combination of both 
values is still contained within the 90 % error ellipse. It is therefore concluded that the 
differences between the model mass functions do not significantly affect the interpretation of 
the HIFLUGCS mass function. Moreover this test can be regarded as confirmation of the 
validity of the PS mass function for the accuracy needed here. 

Also in this work the normalized baryon density may be estimated, using the mean gas 
mass fraction determined in Sect. 16.21 One may set the gas mass fraction in clusters, being 
the largest collapsed objects, equal to the baryon fraction in the universe, i.e. /gas = ^h/^m- 
It has been tested if consistent results are obtained if the value determined for fit, here is 
used for the calculation of model mass functions instead the one given in Sect. 12. 2. ^-i. 21 Using 

—3/2 

= 0.19 ± 0.08 /i^Q ilm results in model mass functions ve rv similar to the ones calculated 



using the baryon fraction given bv lBurles and Tvtlen ((1993). It is therefore not surprising 
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that the best fit values for 0,^ and fig vary only well within the 68 % error ellipse if the former 
rib determination is used. 

It is worth noting that combining these two measurements of J7b; i-e. 0.19ib0.08 h^Q n^n = 
0.0772 h7^, gives further evidence for a low value for by yielding an estimate = 



0.341°-^^ using Hq = 71 kms'^ Mpc-^ ([Mould et al 



120^ . 

termined from the standard deviation of Jib given above. This value for is an upper 
limit since baryons contained in, e.g., the cluster galaxies have been neglected. A low value 
for has b een indicated by this rnethod for sinaller cluster saniples by several works pre- 



vio 



usly (e.g., 'white and FrenlJ 



1991 



Briel et al 



aiier cmster saniples p y 
19921: lB5hringeil (1 9931: 



White and Fabian 1995: D avii 



d et a] 111 99,4 lEvraxdl 119971 : lEtto 



'White et al 



1993bl : 



ri and Fabian 1999; but see 



Sadat and Blanchardll200lh . One has to keep in mind, however, that this estimate extrapo- 



lates the gas fraction from cluster scale to cosmic scales. For the clusters in the sample used 
here it has been found that the gas fraction is not constant but varies with radius and cluster 
mass (Sect. ini21)i therefore further observational tests of this assumption may be useful. 

In summary, even though a conservative estimate has been made by neglecting the possible 
presence of gas temperature gradients, previous es timates obtained from cluster abundance s 
mostly yielded higher values for and ag (e.g..lWhite et al1ll993al : Girardi et al. 1998al ). 
For these works one could expect this already from Fig. l7.1^ and possible reasons for this have 
been discussed in Sect. 17.6.11 The results are, however, in good agreement with th e results 
from the power spectral analysis of the REFLEX clusters. Schuecker et al. ( 2001lj ) find for 
a given ACDM model (Om = 0.3) that erg = 0.7 represents the data well, which is very close 
to (Tg = 0.68 expected using relation (|7.18|) found here. Moreover the (1-cr) range 0.17 < 
rim < 0.37 (using h = 0.71 in their Eq. 18 ) quoted for d irectly is also consistent with the 
90 % range determined here. Furthermore llkebe et al. ( 200ll ^. who analyzed the HIFLUGCS 
temperature function using temperatures from homogeneously reanalyzed ASCA data, find 
Jim = 0.19^'lQ'Qg and ug = 0.96tQ }Q (90% c.l. statistical uncertainty; assuming an open 
cosmology) by comparison with Press-Schechter models, which is in good agreement with 
the results presented here. 



7.6.3 Mass Function estimated using Lx—Mtot Relation 

To show consistency and the power of the Lx-Mtot relation, fits to 'mass' functions, where 
masses have been estimated from the measured X-ray luminosity have also been performed. 
Relations for the flux-limited sample have been used to get the best mass estimate for the 
cluster luminosities included in HIFLUGCS. Mass functions for the two extreme relations 
(M I L) with a = 1.538 and {L \ M) with a = 1.310 are shown in Fig. 17.161 First of all one 
notes the fairly good agreement between the three mass functions. In detail the differences 
between the two mass functions estimated from different luminosity-mass relations can be 
understood by considering that at the low luminosity end the steeper relation predicts a 
higher mass for a given luminosity than the shallower relation, resulting in a shift towards 
higher masses of the mass function. At the high mass side the effect is opposite, resulting in 
a shift towards lower masses for the steeper relation. On average the points for the steeper 
relation lie higher which is caused by the fact that a steeper relation results in a smaller dM 
on average, which gives rise to an increased dn/dM. The differences to the mass function 
calculated using the measured masses are again understood by a similar comparison and are 
partially caused by a possible deviation of the shape of the Lx^Aftot relation from a pure 
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power law. 

Despite these small differences performing an actual fit^ results in the {0,^, ag) values 
(0.14, 0.85) for (L | M) and (0.22, 0.74) for (M | L). The first case is consistent with the 
error range given in ()7.17|) and in the second case the 90 % statistical error ellipse overlaps 
with the 90 % ellipse in Fig. 17.151 From the above and Fig. I7.l41 it is clear that using steeper 
luminosity-mass relations results in higher values for and lower values for erg. Here one 
wants to estimate M from L and therefore (M | L) is the appropriate relation to use. This test 
demonstrates that with in general easy to obtain X-ray luminosities of a statistical cluster 
sample and with the knowledge of the Lx-Aftot relation (even if approximated as simple 
power law) and its scatter realistic constraints on cosmological parameters can be set. 
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Figure 7.16 HIFLUGCS mass function (filled circles) compared to 'mass' functions estimated 
using measured luminosities and luminosity-mass relations (open symbols). Squares have 
been calculated using the (M | L) relation and triangles using the {L \ M) relation for 
H/FLL/GCS clusters. 



7.7 Total Gas and Gravitating Mass in Clusters 

To estimate the fraction of the gravitational mass density relative to the critical density 
contained in galaxy clusters, 

^^Cluster(> Mtot,min) = — / Mtot 'P{Mtot) dMtot , (7.19) 

^For the fit the corresponding scatter in logAf for the two relations (L | M) and (Af | L), criogMtot = 0-22 
and 0.21, respectively (Tab. IS^ . has replaced the mass measurement error, tJAftot.iog = 0.12, in 17.161 . 
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the individual cluster masses divided by the corresponding maximum search volumes have 
been summed up for HIFLUGCS, i.e. 

f^Clustcr = — ^ '200,1 ^ (7.20) 

Note that the determination of riciustor is independent of the Hubble constant. The cumu- 
lative diagram for riciustcr(> ^^20o) is shown in Fig. 17.171 In order to perform a conser- 
vative error estimate, HIFLUGCS has been split into two parts with bu > -|-20deg and 
bii < — 20deg, and the results for these subsamples are also shown in the Figure. This esti- 
mate is conservative because HIFLUGCS is about twice as large as each subsample. Taking 
the second and third lowest mass clusters together with the maximum mass range given by 
their individual uncertainties, one obtains 

^Cluster = 0-012lg;QQ4 (7.21) 

for masses larger than 6.4^Qg x 10^'^ h^Q Mq, i.e. the total gravitating mass contained within 
the virial radius of clusters amounts only to 1.2]']q^ percent of the total mass in a critical 
density universe. Combined with the best estimate from Sect. I7.6~2l = 0.12, this implies 
that about 90 % of the total mass in the universe resides outside virialized cluster regions 
above the given minimum mass. If galaxies trace mass it also follows that by far most of the 
galaxies do not sit in clusters. This result is consistent with the general presumption that 
clusters are rare objects, rare peaks in the density distribution field. 
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Figure 7.17 Gravitational mass density contained in galaxy clusters as a function of minimum 
mass. Filled circles indicate the complete HIFLUGCS, open triangles indicate the 34 clusters 
north of the galactic plane, and open diamonds the 29 clusters at southern galactic latitudes 
included in HIFLUGCS. 
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Comparing the diagram to the mass fraction Ociuste 



0.028^gQQg in clusters with masses 



larger than 2 x 10^^ h7n Mc^ given by iRikugita et alJ |l998) , based on the mass function 
determined by Bahcall and CenI ( 19931 ^. one finds that their estimate is a factor 4-5 higher. 
However, the Bahcall and CeiJ mass function is given for Ma and one gets a consistent result 
if r^ciuster is calculated using the cluster masses formally determined within here. It needs 
to be pointed out that at Ma ~ 2 x 10^*^ U^q Mq the typical virial radius is ~ 1 h'^^ Mpc and a 



mass determination at 3 K^q Mpc based on the assumption of virial equilibrium may therefore 
be rather uncertain and possibly leading to overestimates of J^ciuster- This becomes more 
crucial if mass functions for Ma are extrapolated even down to galaxy masses. For instance 
for the lowest mass group contained in HIFLUGCS one finds that M'^ nn = 3.5 x 10^^ h7a Mg 
while Ma = 
find r^ciust e 



1.2 X 10^^ h^Q Mq, i.e. Ma is a factor 3.4 larger. This wav lRikugita et al.l (|l998l ) 



0.12 lb 0.02 within the mass range 2 x 10 -2 x 10 /i^q Mq, which, compared 
to the results from Sect. I7.6.2( would account already for almost all mass in the universe. 

Replacing M200 with Mgas,200 in (|7.20j) and performing an analogous calculation yields 
the fraction of the total gravitating mass in the universe contained in the intracluster gas of 
galaxy clusters, fib,ciuster- The result is shown in Fig. 17.181 One notes that the cumulative 
curve flattens out with decreasing gas mass more strongly than the analogous plot in Fig. 17.171 
which is caused by the decreasing gas mass fraction with decreasing mass (Fig. 16.2]) . The two 
highest gas mass triangles lie a factor of two higher than the two highest gas mass solid circles 
by construction, since the same cluster masses (originating in one hemisphere) are divided 
by volumes which differ by about a factor of two. Combining the results for the second and 




Figure 7.18 Gas mass density contained in galaxy clusters as a function of minimum gas mass. 
Symbols have the same meaning as in Fig. 17.171 

third lowest gas mass clusters in the same way as for (|7.21() one obtains 

^^b,Cluster = O-OOl^iaoooi ^50^'^ (7.22) 
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for gas masses larger than 6.9"^^'^ 
fraction J7b,ciuste: 



1.5 X 10^^ hJ^^MQ. iFiikugita et al 



n 0044+0-0028 , -1.5 



give a gas mass 



in clusters with gravitational masses larger than 2 x 
10^* h^Q Mq. Using the value they quote for /gas this corresponds to gas masses larger than 
3.2 X 10^^ ^50^'^ again a factor of 4-5 higher than the value determined at this 

minimum gas mass from Fig. 17.181 Their result is still higher but consistent if gas masses 
determined within ta are compared. If f^b, cluster is taken as an estimate of the baryon content 
in clusters it has to be regarded as a lower limit since especially for the lo w mass clusters th e 



cluster galaxies may contribute a non negligible amount of baryons (e.g., David et al. IQQOl ) 



8 Conclusions 



In this work tight constraints on the mean mass density in the universe - one of the funda- 
mental cosmological parameters - have been obtained. Use has been made of high quahty 
X-ray observations of a weh defined sample of galaxy clusters. The strongest constraints ob- 
tained from the cluster mass function = 0.12+°;°^ are well consistent with the contraint 
from the cluster gas mass fraction ^ 0.34 and the estimate from cluster mass to light 
ratios ~ 0.15. Previous measurements of cluster abundances indicated somewhat higher 
values for than found here. Possible reasons for this have been discussed in Sects. 17.6.11 
and l7.6.21 However, the results obtained here are consistent with the results from the power 
spectral analysis based on the 462 galaxy clusters of the REFLEX survey. 

How do the results of this work fit into the general picture? How do they compare to 
constraints derived from completely independent measurements und what conclusions can 
be drawn fro m a comparison? Among the various methods that have been applied (e.g., 
Peebles! ^9^, currently two mor e methods se em especially encouraging for constraining the 



relevant parameters (e.g., Bahca ll et al. 1999f ). As mentioned in the Introduction one ap- 
proach uses measurements of temperature fluctuations of the CMB, another distant SNe as 
standard candles. The constraints (95% c.l.) a chieved recently by these two methods (e.g., 
de Bernardis et al. 20001 : Perlmutter et al. 199fl[ ) are compared in an Qm-^^A diagram to the 
results obtained here in Fig. 18.11 The first noteworthy aspect is that these independent 
approaches nicely complement each other in the way they constrain different regions in pa- 
rameter space. Secondly it is found that all three error ranges overlap; there is an area for 
which all measurements are consistent. This consistency is very encouraging. 

Note that Fig. l8.1l shows the statistical (90%) error range of the measurements performed 
here. Possible systematic uncertainties arise from the cluster mass measurements, from the 
sample construction, and from the formalism used for the comparison to predicted mass 
functions. However, several sources of possible systematic errors have been checked and the 
introduced uncertainties have been found to be smaller than the statistical ones. This gives 
strong confidence to the results. 

What are the implications of the fact that the concordance is confined to a specific region? 
In Fig. 18.11 some additional lines are drawn that separate special regions. The region of 
consistency is very close to the line -|- = 1, i-e. = (Eg. 12.30]) . This indicates that 
the universe may have a flat geometry. Furthermore the region is located clearly within a 
parameter region where the universe will never collapse again, but expand infinitely (Jim ^ 
1 A Q\ > 0). Moreover the region lies well above the line that separates the states in which 
the expansion of the universe is decelerating or accelerating (2J7a = ^^m)- This implies that 
not only will the universe expand forever but it will do so ever faster. 

Further it is assuring to note that the concordance region indicates that the universe is 
older than H^^ = 13.8 Gyr (for Hq = Tlkms"^ Mpc"^)^ since the age inferred for the oldest 
globular clusters ranges between 11.5-15.8 Gyr (e.g.. iBolte and HoganI Il99.4 ISalaris et al 
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19971 : ICarretta et alJEonol 'l 



The results on the mass to hght ratios and gas mass fractions indicate that the nature of 
most of the matter in the universe is as yet unknown. Moreover most recent measurements 

j 1, i.e. a geometry close to bei i ig fiat (e.g 



de Bernardis et al.l 20nnl 


Balbi et al.l 20001: 


2001 :lde Bernardis et al. 


20011. Combining 



in this work one finds i}\ ~ 0.88. This implies that the observed matter accounts only for a 
tiny fraction of the essence that determines the geometry of the universe. 

Observational progress continues in all wavelengths. Improvements on the accuracy of the 
derived parameters from the X-ray side are achievable in the near future, e.g., by a comparison 
of the local mass function presented here with mass functions of well defined higher redshift 
samples of massive clusters. A large spanned redshift interval allows to further weaken the 
degeneracy bet ween and ag, enabling a, more ind.GpGnd.Gnt niGcisurGniGnt of (^-S*? 
|PerreTiodlfl98nl v The reason lies in the fact that the evolution of clusters depends on f^m. 
For this purpose observations of the most luminous REFLEX clusters by XMM-Newton are 
already scheduled. 



CJ 




Figure 8.1 Observational constraints on two fundamental cosmological parameters. The area 
labeled galaxy clusters shows the constraints obtained in this work. 



9 Summary 



An X-ray selected and X-ray flux-limited sample comprising the 63 X-ray brightest galaxy 
clusters in the sky (excluding the galactic band, called HIFLUGCS) has been constructed 
based on the ROSAT All-Sky Survey. The flux limit has been set at 2 x 10"-^''" ergss"^ cm~^ 
in the energy band 0.1 — 2.4 keV. It has been shown that a high completeness is indicated by 
several tests. Due to the high flux limit this sample can be used for a variety of applications 
requiring a statistical cluster sample without any corrections to the effective survey volume. 

Mainly high quality pointed observations have been used to determine fluxes and physical 
cluster parameters. It has been shown that a tight correlation exists between the X-ray 
luminosity and the gravitational mass using HIFLUGCS and an extended sample of 106 
galaxy clusters. The relation and its scatter have been quantified using different fitting 
methods. A comparison to theoretical and numerical predictions shows an overall agreement. 
This relation may be directly applied in large X-ray cluster surveys or dark matter simulations 
for conversions between X-ray luminosity and gravitating mass. 

Data from the performance verification phase of the recently launched X-ray satellite 
observatory XMM-Newton on the galaxy cluster Abell 1835 has been analyzed, in order to 
test the assumption of isothermality of the cluster gas in the outer parts applied throughout 
the work. It has been found that the measured outer temperature profile is consistent with 
being isothermal. In the inner regions a clear drop of the temperature by a factor of two has 
been found. 

Physical properties of the cluster sample have been studied by analyzing relations be- 
tween different cluster parameters. The overall properties are well understood but in detail 
deviations from simple expectations have been found. It has been found that the gas mass 
fraction does not vary as a function of intracluster gas temperature. For galaxy groups 
(^B^x ^ 2keV), however, a steep drop of /gas has been observed. No clear trend of a varia- 
tion of the shape of the surface brightness profile, i.e. /3, has been observed as a function of 
temperature. The Lx-Tx relation has been found to be steeper than expected from simple 
self similar models, as has been found by previous authors. But no clear deviations from 
a power law shape down to /cb^x = 0.7 keV have been found. The Mtot-Tx relation found 
here is steeper than expected from self similar models and its normalization is lower com- 
pared to hydrodynamic simulations, in agreement with previous findings. Suggested scenarios 
to account for these deviations, including heating and cooling processes, and observational 
difficulties have been described. It appears that a blend of different effects, possibly includ- 
ing a variation of mean formation redshift with system mass, is needed to account for the 
observations presented here. 

Using HIFLUGCS the gravitational mass function has been determined for the mass 
interval 3.5 x 10^^ < M200 < 5.2 x 10^^ h^^ Mq. Comparison with Press-Schechter mass 
functions has yielded tight constraints on the mean matter density in the universe and the 
amplitude of density fluctuations. The large covered mass range has allowed to put constraints 
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on the parameters individually. Specifically it has been found that fim = 0.12_qq4 and 
(78 = 0.96^q}2 (90% c.l. statistical uncertainty). This result is consistent with two more 
estimates of obtained in this work using different methods. The mean intracluster gas 
fraction of the 106 clusters in the extended sample combined with predictions from the theory 
of nucleosynthesis indicates ^ 0.34. The cluster mass to light ratio multiplied by the 
mean luminosity density implies ~ 0.15. Various tests for systematic uncertainties have 
been performed, including comparison of the Press-Schechter mass function with the most 
recent results from large A^-body simulations, yielding deviations smaller than the statistical 
uncertainties. For comparison the best fit values for fixed as values have been determined 
yielding the relation erg = 0.43 

The mass function has been integrated to obtain the fraction of the total gravitating mass 
in the universe contained in galaxy clusters. Normalized to the critical density it has been 
found that f^ciuster = 0.012+o;oo4 for cluster masses larger than 6.4;o;g x 10^^ h^^ Mq. With 
the value for Q,^ determined here this implies that about 90 % of the mass in the universe re- 
sides outside virialized cluster regions. Similarly it has been found that the fraction of the to- 
tal gravitating mass which is contained in the intracluster gas, J^b,Ciuster = O-OOISI^qqqq^ ^50^^ 
for gas masses larger than 6.9lj^ 5 x 10^^ h^o^^ -^0' ^^^^ small. 
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